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ABSTRACT 


A numerical implementation of design sensitivity analysis of built- 
up structures is presented, using the versatility and convenience of an 
existing finite element structural analysis code and its database 
management system* The finite element code used in the implementation 
presented is the Engineering Analysis Language (EAL), which is based on 
a hybrid method of analysis. It has been shown that design sensitivity 
computations can be carried out using the database management system of 
EAL, without writing a separate program and a separate database. 

Conventional (sizing) design parameters such as cross-sectional 
area of beams or thickness of plates and plane elastic solid components 
are considered in this report. Compliance, displacement, and stress 
functionals are considered as performance criteria. The method 
presented in this paper is being extended to Implement shape design 
sensitivity analysis using a domain method and a design component 
method. Results of shape design sensitivity analysis will be reported 
in Part II. 
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CHAPTER I 
INTRODUCTION 

1 . 1 » Purpose 

To date there exist a wide variety of finite element structural 
analysis programs that are used as reliable tools for structural 
analysis. They give the designer pertinent information such as 
stresses, strains, and displacements of the structure being modeled. 
However, if this information reveals that the structure does not meet 
specified constraint requirements, the designer must make intuitive 
estimates on how to improve the design. If the structure is complex, it 
becomes very difficult to decide what step must be taken to improve the 
design. There is however, substantial literature [1] on the theory of 
design sensitivity analysis, which predicts the effect that structural 
design changes have on the performance of a structure. 

The purpose of this work is to develop and implement structural 
design sensitivity analysis, using the adjoint variable method that 
takes advantage of the versatility and convenience of an existing finite 
element structural analysis code and its database management system and 
the theoretical foundation of structural design sensitivity analysis 
that is found in Ref. 1. The finite element code that will be used is 
the Engineering Analysis Language EAL [2]. 
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Using the full capabilities of the EAL system, design sensitivities 
can be calculated within the program, without knowing the source code of 
the program. This has the advantage that the user deals with only one 
program, with only one database and no interfaces between different 
programs [3,4,5]. 

1.2. Continuum Approach of Design Sensitivity Analysis 

A number of methods could be used to implement structural design 
sensitivity analysis, but the most powerful is the continuum approach 
[1]. This method can be implemented outside an existing finite element 
code [3,4,5], using only postprocessing data. This approach is 
convenient, because design sensitivity analysis software does not have 
to be embedded in an existing finite element code. The continuum 
approach to design sensitivity analysis calculation can also be 
implemented using a powerful database management system such as the 
Engineering Analysis Language (EAL). Using the database management 
system of EAL, only one database is necessary for computation of design 
sensitivity information. That is, it is not necessary to create 
interfaces and other datafiles to compute sensitivity information. 
Information on element shape functions used in the finite element model 
is, however, necessary for design sensitivity calculation. 

The continuum approach to design sensitivity analysis can easily be 
extended to complex structural systems that have more than one 
structural component [6]. The design sensitivity vector is the 
derivative of a constraint functional with respect to design 
parameters. The magnitude of each component reflects how sensitive the 
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constraint functional is to a change in the associated design 
parameter. If a component of the vector is negative, the corresponding 
design parameter should be decreased to increase the value of the 
constraint functional. If a component of the vector is positive, the 
corresponding design parameter should be increased to increase the value 
of the constraint functional. In addition, if the magnitude of a 
component of the vector is large, then the corresponding design 
parameter will have a more substantial effect on design Improvement. 

When a designer uses a finite element structural analysis code in 
design of a structure, it is most likely that a number of program runs 
are necessary before a substantially Improved design is obtained. With 
the aid of a design sensitivity vector, the designer will know what 
direction to take to improve the design most efficiently. 
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CHAPTER II 

DESIGN SENSITIVITY ANALYSIS METHOD 

A detailed treatment of methods and procedures for calculating 
design sensitivities are given in Ref. 1, for constraint functionals 
such as compliance, displacement, stress, and eigenvalues. For 
compliance and eigenvalue functionals, the adjoint equation is the same 
as the state equation, thus no adjoint equation needs to be solved. 

Each displacement and stress functional requires an adjoint load 
computation and an adjoint equation must be solved. Note that the state 
equation and the adjoint equation differ only in their load terms. 

Using the reload or multiload option of an existing finite element code, 
the adjoint equation can be solved efficiently [3]. For the 
displacement functional, the adjoint load is a unit point load acting at 
the point where the displacement constraint is imposed. For the stress 
functional the shape function of the finite element code is necessary to 
calculate equivalent nodal forces of the adjoint load. 

The equivalent nodal force computation of the adjoint load for the 
stress functional can be based on different principles. To be 
consistent with EAL, a hybrid formulation should be used, which requires 
computation of the applied loads in terms of stress coefficients. This 
is impractical, because EAL is formulated in terms of nodal displacement 
coefficients. Another method, which is consistent with the hybrid 
formulation of EAL, is based on the modified Hellinger-Reissner 
principle [10,11]. Here the equivalent nodal forces are expressed in 
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terms of nodal displacement coefficients. To alleviate this difficulty, 
an equivalent nodal force computation Is based on the formulation used 
In the displacement method of finite element analysis. That Is, once 
the degrees of freedom (nodal displacements) of the element In EAL are 
known, a compatible displacement shape function that Is defined on the 
same finite element and has the same degrees of freedom can be used for 
computation of equivalent nodal forces. 

The flow chart of Fig. 2.1 shows the overall procedure. First, the 
model Is defined by Identifying the design variables, constraint 
functionals, finite element model, and loadings. In the next step, EAL 
Is used to obtain structural response. With the structural response 
obtained, an adjoint load is calculated, external to EAL, using assumed 
displacement shape functions. The adjoint load is input to EAL, to 
obtain an adjoint response for each constraint using reload option. 

Using the original structural response and the adjoint response, design 
sensitivity Information is computed for each constraint, by numerically 
Integrating the design sensitivity expressions. The process is 
convenient, since it uses data that are available or easily computable 
In the database. 

To give the basic idea of Implementation of design sensitivity 
analysis and computation procedures, simple prototype structural 
components are Investigated. Once design sensitivity analysis of 
structural components Is completed, the design component method of 
Ref. 6 can be used for design sensitivity analysis of built-up 
structures. The procedures and equations necessary for analysis of 



F£ Model Definition 
Preprocessing 



Figure 2.1. Flow Chart of Design Sensitivity Computation 
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structural components as membranes, plates, and beams are shown in the 
following sections • 


2mlm Membranes 

Consider a thin elastic clamped solid shown in Fig. 2.2, with 
thickness u = h(x) of the membrane as a design parameter where x * (x^, 
X 2 )* The energy bilinear form and load linear form are [1], 

__ 2 .. 

a (z,z) = // h(x) I ^ ^(z) £ ^(z)d^ (2.1.1) 

^ « i,j=l 

and 

2 2 

^ G) - a I F^z^dQ + / I T^z^dF (2.1.2) 

n 1=1 r 1=1 

where z [z^, z^j"^ is the displacement field, F = [F^, F^]^ is the 
applied body force, T = [T^, is the boundary traction, and o^^(z) 

and are the stress and strain fields associated with the 

displacement z and the virtual displacement z, respectively. The 
variational state equation is [1] 

a (z,z) = ^ (z) (2.1.3) 

u u 

for all kinematically admissible virtual displacements z. 

First consider the functional that represents compliance of the 
structure, 

’(',=// I F^z^dQ +11 T^z^dF (2.1.4) 

^ a 1=1 F i=i 

The first variation of Eq. 2.1.4 is 

I F^z^'dfi + / I T^z^ dF 

^ a i=i F i=i 


(2.1.5) 
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Figure 2.2. Clamped Plane Elastic Solid of Variable Thickness 

In order to eliminate the dependence on variation of the state variable 
in Eq. 2.1.5, it is necessary to define the adjoint equation as [1] 

2 2 

a (^,'^) = // I + I I T^l^dT (2.1.6) 

" i»l r i=i 

for all kinematically admissible virtual displacements Since the 
right side of Eqs. 2.1 ,3 and 2*1. 6 are identical if ^ = z, the adjoint 
equation does not need to be solved. Using the adjoint variable method, 
the design sensitivity is [1] 

2 

= -//[ I e^j(2)]6hdn (2.1.7) 

n i,j=l 

since z = ^ for the compliance functional. 
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To numerically integrate Eqs. 2.1.4 and 2.1.7, a two-by-two Gauss- 
point integration procedure is used. The equations then become 


N 

= I 

k=l 


1 


I [ 1 (F^)(z^)]w^J^ 

k 

I 

1 

N 

H 

1 

^1 h= 1 

k 

^=1 h=l J 


( 2 . 1 . 8 ) 


and 




N 

I 

k=l 


I [- I 

M i,j=l 


6h' 


(2.1.9) 


respectively, where J is the Jacobian, N is the total number of 
elements, superscript ^ is the counter for the number of Gauss points, 
superscript k is the counter for the element number, W is the weighting 
constant for the ^th Gauss point, and superscript h is the direction of 
the force and the displacement. 

Since EAL gives only the boundary displacement shape function for 
the 4-noded membrane element E41, a bilinear shape function for 
displacement is adopted for integration in Eqs. 2.1.4 and 2.1.7. Using 
stress computation of membrane element E41 in EAL, stresses and strains 
can be expressed in matrix form as 


{<^^^(z), O 22(2), <7^^(z)F = [p]{3} 

{e^^(z), e^^(z), e^^(z)r = [e]"M<^} 


where 


[p] 


0 0 1 X 2 0 

0 1 0 0 

1 0 0 0 0 


( 2 . 1 . 10 ) 

( 2 . 1 . 11 ) 


( 2 . 1 . 12 ) 


10 


{3} « fbj, b2, b3. 




V 0 



(2.1.13) 


(2.1.14) 


and [P] is the position coordinate matrix, which determines points where 
the stresses are obtained, {3} is the stress coefficient vector, and [E] 
is the elasticity matrix for a plane stress problem [8,9]. 

Next consider the functional that represents the displacement z at 

A 

a discrete point x. 


A A A 

= z(x) = // fi(x - x)z(x) d^ 


(2.1.15) 


A 

where 3(x) is the Dirac delta. The first variation of Eq. 2.1.15 is 


’t'l = // «(x - x)z'(x)dO (2.1.16) 


The adjoint equation in this case is [1,3] 


a (^,^) = // ^(x - x)^(x)d^ (2.1.17) 


for all kinematically admissible displacements This equation has a 
unique solution which is the displacement field due to a unit point 

A 

load acting at a point x. Using the adjoint variable method, design 


sensitivity of the displacement functional is 


n = -// I a^J(z) e^^(X) 6hd« (2.1.18) 

^ i,j=l 

where ^ is the solution of Eq. 2.1.17. 
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Using the same numerical techniques applied in the compliance 
constraint case, Eq. 2.1.18 can be evaluated as 


N 

= I 

k=l 


where 


I [- I 

M i,j=l 


1^^)} = [eTM^ca)} 


6h‘^ 


[E]-'[p]{a} 


(2.1.19) 


( 2 . 1 . 20 ) 


Finally consider the general functional that represents a locally 
averaged stress on an element as 


K - ff g(®(z)]m_d£2 


( 2 . 1 . 21 ) 




where mp is a characteristic function, defined on a finite element as 
1 


”P 




If dii 
a 

p 

0 




( 2 . 1 . 22 ) 




and g is the stress function. The first variation of Eq. 2.1.21 is 


= III 1 -^z) o^^(z’)]m d« 

3 a i,j=i 30^-J ^ 


(2.1.23) 


Replacing the variation in state z* by a virtual displacement the 
adjoint equation is obtained as [1,3] 


== //[ I ^f^ (z) o^^(l)]mpdJ2 


(2.1.24) 


Q i,j=i aoiJ' 

for all kinematically admissible displacements Equation 2.1.24 has a 



12 


unique solution for a displacement field Using the adjoint variable 
method, design sensitivity of the stress functional is 
2 

'(';=-//[ I e^^(^)]'Shd!^ (2.1.25) 

a i,j=l 

where ^ is the solution of Eq. 2.1.24, 

Using the same numerical techniques applied in the compliance 
constraint case, Eqs. 2 •1.21 and 2.1.25 become 

K = I lg^(Az))m W^]j^ . (2.1.26) 

^ JUi ^ 


and 


3 


N 

I 

k=l 


I [- I ei^(A)]V^^ 

M l,j-l 




(2.1.27) 


The right side of Eq. 2.1.24 can be evaluated as 


I [ 1 -^z) <^^^(^)] (2.1.28) 

i,j»l P 


After the adjoint load of Eq. 2.1.17 is calculated, the adjoint 
displacement field and resulting adjoint stress and strain fields are 
evaluated. The adjoint strains are then used in evaluating Eq. 2.1.27 
for sensitivity of the constraint functional. 

If von Mises’ stress criteria is selected in the constraint 
functional , then 

g = + 3(0^^} - ^ (2.1.29) 


and 
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3all 

3o22 


(2«j“ - <J^^)/2g 


(2g 22 _ a^l)/2g 


(2.1.30) 


3a 


12 


3<7^^/g 


which can be written in vector form as 


. {Js- _is_ _k_}T 

3aiJ S^ll’ aa22* 9^12^ 


(2.1.31) 


The equivalent nodal force is computed, based on the modified 
Hellinger-Reissner principle. The stress can be written as 


- [P]{3} 


(2.1.32) 


(see also Eq. 2.1.10). The stress coefficients {3} can be expressed in 
terms of the nodal displacement coefficients {q} as [7] 


Is} - [ h -‘]( t ](,1 


(2.1.33) 


if {7a} a (o|, where 


3a^^ 3a22 3a^2j.j. 


- 1 ^- 

^1 ^2 i 

fVa} 

can be written as 


i=l,2 


(2.1.34) 


M = [Vp]{6} 


(2.1.35) 


where 


[7p] 


3 , 
*1 

^i 


0 

1 

0 


1 y-2 0) 


0 0 x^) 


0 0 0) 


i=l,2 
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Using Eq. 2.1.34, [^p] in Eq. 2.1.35 becomes zero. Note that, if [^p] 
is not zero, Eq. 2.1.33 is not valid. Using Eq. 2.1.33, Eq. 2.1.32 can 
be written as 

fo} - [p][H~M[T]{q} (2.1.36) 

Using Eqs. 2.1.27, 2.1.31, and 2.1.36, the equivalent nodal forces {p} 
can be written as 

{f} - I {-^F[p][h”M[t]wS^ (2.1.37) 

JUl 

which is valid only if [^p] is zero. 

2.2. Bending of Plates 

Consider the clamped plate in Fig. 2.3 of variable thickness 
u - t(x), with a distributed load f(x). For this design independent 



Figure 2.3. Clamped Plate of Variable Thickness 
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loading, the energy bilinear form and the load linear form for the plate 
are given as [1] 

A 

a^(z,7) = // D(u)[zjjlj^ + Z 22^22 + ^(^22^11 ==11^22^ 

+ 2(1 - ( 2 . 2 . 1 ) 

and 

^(z) =* // fzd^ (2.2.2) 

a 

A o 

where D(u) * Et/[l2(l - '^ )] is flexural rigidity, E is Young ^s 
modulus, is the Poisson ^s ratio, and f is externally applied 
pressure. The governing state equation is [1] 

a (z,z) = ^(z) (2.2.3) 

u 

for all kinematically admissible displacements z. 

First consider the functional that represents compliance of the 
structure, 

= // fzd« (2.2.4) 

The first variation of Eq. 2.2.4 is 

’I'! = // fz'dS (2.2.5) 

The adjoint equation is defined as [1] 

a (X,"^) = // f^d^ (2.2.6) 


for all kinematically admissible displacements As in the membrane 


case, the right sides of Eqs. 2.2.3 and 2.2.6 is identical, so 1 = z. 
Using the adjoint variable method, design sensitivity of the compliance 
functional is [1,3,4] 

= // {-Et^[z^j + Z 22 + 2^z^Z 22 + 2(1 - v)zj2]/^(l " ^td^ 

2 

= - // [ I e^^(z)] «td« (2.2.7) 

“ l.j=l 

where ^^^(z) and e^^(z) are stress and strain of the extreme fiber, 
given as 


eij = - 



i.j “ 1, 2 


( 2 . 2 . 8 ) 


and 


,11 

Et 

<'11 

+ VZ 22 

2(1 

- v^> 

f22 

Et 

<"22 

^ "^11 

2(1 

- 

,12 

Et 



2(1 

- v2) 

^12 



(2.2.9) 


To evaluate Eqs. 2.2.4 and 2.2.7 numerically, a two-by-two Gauss- 
point integration procedure is used. Equations 2.2.4 and 2.2.7 become 


^ - hi 

^ k=i 


. ( 2 . 2 . 10 ) 


and 


’<'] = - I { I[ I (<J^^(z)]^[e^^(z))^]w^J^}^«t’^ (2.2.11) 

k>l A=1 i,j=l 

For the Gauss integration, and are obtained at each Gauss 


point, using the same moment formulation as in EAL* The moment vector 
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[M] is defined as 


{m^^(z), m^^(z), m^"^(z)P = [p]{e} 


-22 


12, 


( 2 . 2 . 12 ) 


where 


[p] 

{0} 


10 0 
0 1 0 
0 0 1 


X 2 0 


Xj 0 0 


x^X2 0 


Xj^ 0 X2 0 0 0 

X2 Xj^ 0 


X1X2 


(2.2.13) 


0 0 0 0 X, 0 0 

{bi, b2> b^j b^, bj, b^, by, bg, bg, ^10* (2.2.14) 


Thus, 


{o^^(z), <r^{z), {m^\z), M^^(z), -M^^(z)P (2.2.15) 

{e“(z), e^^(z), e^2(z)F= [eI'M^^Vz), (2.2.16) 


where [E] ^ is given in Eq. 2.1.14. 

Next consider the functional that represents displacement z at a 

A 

discrete point x. 


= z(x) = // fi(x - x)z(x)d^ 
^ Q 


The first variation of Eq. 2.2.17 is 


= // ^(x - x)z*(x)dQ 
^ Q 


The adjoint equation is defined as [1] 


(2.2.17) 


(2.2.18) 


a (^,^) = // ^(x - x)^(x)d^ 
“ 0 


(2.2.19) 


for all kinematically admissible displacements This equation has a 
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unique solution which is the plate displacement field due to a unit 

A 

vertical load at point x. Using the adjoint variable method [1], design 
sensitivity of the displacement functional is 


'll’ = // {— Et^[z ^ + z ^ + ^(z ^ + z ^ ) 

”c l ^1- 11 <10 '^i 1 11 n't' 


'11 11 22 22 


11 11 22 22 ' 


+ 2(1 - 


( 2 . 2 . 20 ) 


Using Eqs. 2.2.8 and 2.2.9, Eq. 2.2.20 can be rewritten as 


K ‘ - ■ff ^ 'J^^(z) e^J(X)6tdJ2 

n i,j=l 


( 2 . 2 . 21 ) 


Using the same numerical integration as in the compliance case, 
Eq. 2.2.21 becomes 


n- Il-I[ i (2.2.22) 

k=l ^=1 i,j=l 

where are strains obtained by applying the adjoint load. The 

procedure is the same as in Eqs. 2,2.12 thru 2.2.16, 

Finally, consider the functional that represents a locally averaged 
stress at the extreme fiber of the plate. 


'f'c “ // g(®(z)]m d^i (2.2.23) 

0 Q ^ 


where g[®(z)j is chosen as the von Mises' stress criteria and mp is a 
characteristic function on finite element ^p, defined as 


m^ 



X C % 


(2.2.24) 


0 
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The first variation of Eq. 2«2«23 is 


r.- a [ I -|f dn 

- i,j-l 3a J(z) 




The adjoint equation is defined as [1] 


(2.2.25) 


2 

a (X,X) = // [ ! ^f7(0 a^^(X)]m dJl (2.2.26) 

" ft i,j-l 3a^J 

for all kinematically admissible displacement X. Using the adjoint 
variable method, design sensitivity of the stress functional is 
2 

= // - [ I a^^(z) e^^(X)]6tdft + // -|^(a(z))mpdft (2.2.27) 
ft l,j-l ft 

where X is the solution of Eq. 2.2.26. Using two-by-two Gauss-point 
integration, Eqs 2.2.23 and 2.2.27 can be expressed as 


>^6 


J [g(a\z)Wlw^J^ 

A-1 


(2.2.28) 


and 


’I'l - I { I[- l (a^^(z))^(e^J(X))^]w^J^ + y 

k=l £=1 l,j=l i=l 

(2.2.29) 


(2.2.30) 


The right side of Eq. 2.2.26 can be written as 

I [ I -^(2) a^^(X)]Vw^J^ 

11=1 i,j=l 3o J 

For the von Mises’ stress criteria, the last term on side right side of 
Eq* 2*2.29 becomes 

11=1 


(2.2.31) 
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For the displacement method, the equivalent nodal force for the 
adjoint load can be computed in a consistent way, by using the 
displacement shape functions of the code and then using the same 
procedure as in Refs. 3,4 and 5. For the hybrid method, Reissner's 
principle can not be used, because [^p] is not necessarily zero. To 
alleviate this difficulty, it is proposed to select an acceptable 
displacement shape function for the adjoint load calculation. That is, 
once degrees of freedom of the element in EAL are known, a compatible 
shape function that has the same degrees of freedom can be selected. 

The term — f-r is the same as for membrane elements (Eq. 2.1.31). 
Stresses in a displacement formulation are given as [8,9] 


l"l- UlUlUl 


(2.2.32) 


Using Eqs. 2.1.31, 2.2.29, and 2.2.32, the equivalent nodal forces {f} 
can be evaluated as 


{f} = I [E][B]mpwS^ 


(2*2,33) 


2 « 3 e Beams 

Consider a beam with variable width and height, as shown in Fig. 2.4. 

T 

Width and height are the design parameters; l.e., u = [b(x 3 ),h(x 2 )l . 

The energy bilinear form and the load linear form of the beam are 


a (z,z) = / E[bh{w-w.j + 
m Q 3 3 


LV33V33)) 


3 ^ X 


'“ 33 ^ 33 ^^^*3 


(2.3.1) 


and 

_ i _ 

Mz) = J F"^zdx» 
0 


(2.3.2> 
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where F - [f^, f^] ig the distributed load, E is Young's modulus, 

3 3 

hb /12 is the moment of inertia with respect to Xj^— axis [2], bh /12 is 

the moment of inertia with respect to X 2 ~axis, z “ [u,v,w]'^ is the 





h 




Figure 2.4. Beam of Variable Thickness and Width 

displacement field, in Xj“, X 2 ~, and Xj -directions , respectively, 

are curvatures of the displacement field, and 

^ iT 

u^2 curvatures of the virtual displacement field z = Lu,v,wJ • 

In EAL [2], the beam element E21 is based on Timoshenko beam theory 
and includes torsion. The derivation here is based on technical beam 
theory without torsion. Thus, if the energy bilinear form of Eq. 2.3.1 
is used, results of sensitivity analysis are applicable only if shear 
and torsion are negligible. However, it can be easily extended to 
include torsional effects. 

Note that the load is applied in the positive direction of the 
coordinate axis. The state equation is [1] 
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a (z,z) = Hz) 
u 


(2.3.3) 


for all kinematically admissible displacements z. 

First consider the functional that represents compliance of the 
structure, 

I 

= / F^zdx_ (2.3.4) 

' 0 

The first variation of Eq. 2.3.4 is 
I 

'I'i - / F^z'dx- (2.3.5) 

0 

By replacing the variation of z' in Eq. 2.3.5 by a virtual 
displacement an adjoint equation is defined as [1] 

I ■ 

a (^,"^) ” / F"^^dx (2.3.6) 

“ 0 

for all kinematically admissible displacements Since Eq. 2.3.2 is 
identical to Eq. 2.3.6, ^ = z. Using the adjoint variable method of 
design sensitivity analysis. 




2 . 3Eb^h 2 . Eh^ ~2 

- J lEhw^ + -jx — V,, + u,,|4bdx, 

0 


33 12 


133J 


f 2 ^ Eb^ 2 ^ 3Ebh^ ~2 1 r. , 

- MEbw^ + j2 ^33 + 12 U33l^hdx3 


(2.3.7) 


A cubic displacement shape function is assumed for beam bending and a 
linear displacement shape function is assumed for axial displacement. 

To numerically integrate Eqs. 2.3.4 and 2.3.7, a two-point Gaussian 
integration procedure is used. Equations 2.3.4 and 2.3.7 can be 


evaluated as 
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7 


^ r ? At i 4 

I [ I F z J W j 
k=l A=i 


(2.3.8) 


and 


- II I E[h(w‘)2 , ‘ )2 . 

' k-1 A-1 ^ 4 33 12 33 

- I{ I E[b(w3>^ + |^V33)^ + (2.3.9) 

k= 1 1 


where N is the total number of elements, ^ is the Gauss point counter, W 
is the weighting constant for the Gauss point, and J is the Jacobian. 
Next, consider the functional that represents displacement z at a 

A 

discrete point X 3 , 

A ^ A A 

'I'g = z(x3) = / - X3)z(x3)dXg (2.3.10) 

A 

where ^ is the Dirac delta. The first variation of Eq. 2.3.10 is 

JIa 

’I'g = / '^(Xg - Xg)z'(x3)dx3 (2.3.11) 

The adjoint equation is defined as [1] 

a(X,T) s / 6 (x^ - Xg)^(x 3 )dx 3 (2.3.12) 

for all kinematically admissible displacements Equation 2.3.12 has a 

A# 

unique solution where ^ is beam displacement due to a 

A 

unit load acting at a point x^. Using the adjoint variable method the 
design sensitivity analysis is 
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= - / {Ehw3X^ + ^ V33X^3 + ^ u33A^3}6bdx3 


■ ff- ^33^33 ■" ^ ^33^33^^^^’'3 (2.3.13) 


Using the two-point Gaussian integration procedure, Eq. 2.3.13 


V ( Y itFu J. ^ J. Ill X^^1 it^t^1^6v>^ 

" , V? 1 ^^^3^3 4 ’^33^33 12 '^33^33^'^ ^ ^ 

k=l ^*1 


V 1 Y wfu ^x'^^ J. ^ X^^ J. t)h^ fo -i 1/X 

' iM , 1 ^^"3^3 12 ^^33^33 12 ^^33^33^^ ^ ^ (2.3.14) 

k=l *=1 


Finally, consider the functional that represents extreme fiber 
stresses in the beam. 


’f'o • / {w, - sign(b) y V 23 - sign(h) u^^ lEm^dx^ (2.3.15) 


rq J l"0 
0 


2 “33'"‘“P“"3 


where h /2 is the half-depth of the beam, b /2 is the half width of the 
beam, sign(b) and sign(h) are +1 or -1 and indicate at what extreme 
fiber the sensitivity is computed, and mp is a characteristic function 


that is defined on a finite element dXp as 


/ , dx 

dxp 


xCdx, 


X Gdx^ 


(2.3.16) 


The first variation of Eq 2.3.15 is 
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n, 

i>g = I Eniplwj - slgn(b) j - sign(h) y }dx. 


- I Em (sign(b) — ^ + slgn(h) 


u 
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5h)dx- 


(2.3.17) 


Replacing the first variation z' in Eq. 2.3.17 by a virtual 
displacement *= [^, the adjoint equation is defined as 


a(^,^) = / Enipl^ - slgn(b) J ^3 “ sign(h) (2.3.18) 

for all kinematically admissible displacements Equation 2.3.18 has a 
unique solution for a displacement field 1. Using the adjoint variable 
method of design sensitivity analysis, 

*ffi-:33X"3|6bdX3 • 


Z Z 3 

- f EnipSign(b) 5bdx^ - / {shw^l^ ^33^33 


0 


2 ~ 


Ebh 
12 “33 33 


}5hdx, - J Em sign(h) 


u 


0 


1^ ShdXg 


(2.3.19) 


where ^ is the solution of Eq. 2 .3. 18. With the two-point Gauss 
integration procedure, integrals in Eqs. 2.3.15 and 2.3.19 are evaluated 
as 


2 

’•'q “ ^ “ Sign(b) Vj 3 - sign(h) ^ (2.3. 


i:i ‘'"P 


2 “33 


20) 


I V J T r hb^ i h^ f k 

'^9 “ ■ uV ^ ^^33^3 J2 “33 ^3^^^ ^ 

k=l *=1 


bh^ “ 


~ . 0^,^ ^^^3^3 12 ^^33^33 ''' 4 “33^33^^ ^ ^ 

k=*l *=1 


^Z 


- ! I I Em„ ^ slgn(b) - ! 1 I Em^ ^ slgn(h)wb 


k=l i=l 


P 2 


k=l i=l 


(2.3.21) 
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The right hand side of Eq. 2.3.18 can be written as 

N - 2 Wn V V, £ £.^ 

Z tnp{ I E (^23 “ sign(b) Y ^33 - sign(h) y A^^)w J } (2.3.22) 

k=l 

2 >4. Built-Up Structures 

A general structure is a collection of structural components that 
are interconnected by kinematic constraints at their boundaries. 

Results stated in this section are from Refs. 1 and 6. The energy 
bilinear form of a general system, consisting of beams, membranes, and 
plates can be written as 

a^(z,z) = I [a^(z,z)]^ + I [a^(z,z)]^ + I [a^(z,z)j® (2.4,1) 

where [a (z,z)]^, [a (z,z)]^, and [a (z,z)]^ are given in Eqs. 2.1.1, 
u u u 

2.2.1, and 2.3.1, respectively. The load linear form of a general 
system can be written as 

^ (I) = I G)]^ + I (T)F + I (z)]® (2.4.2) 

u u u u 

where [^^(z)]^, [^^(z)]^, and [^^(z)]^ are given in Eqs. 2.1.2, 2.2.2, 
and 2.3.2, respectively. The state equation is [1] 

a (z,z) = ^ (z) (2.4.3.) 

u u 

for all kinematically admissible virtual displacements z. Since the 
energy bilinear and load linear forms of the state equation are just the 
sum of energy bilinear and load linear forms of each structural 
component, the design sensitivity equation of the system is a simple 
additive process [1,6]. The generalized design sensitivity of a built- 


up structure is 



27 


r = K^b + rsh + rst (2. 4,4) 

D n t 

which is the sum of the sensitivities of each structural component, 
given in Sections 2*1, 2.2, and 2.3, respectively. 

2.5. Coupling of Bending and Membrane Effect 
2.5.1. Formulation of Membrane Plus Bending 
A clamped plate combining bending and stretching is shown in Fig. 
2.5 with the laterally distributed load f(x) and in-plane traction load 
T = 



Assuming that bending and stretching are decoupled, one can obtain the 
energy bilinear form by adding the plate and membrane energy bilinear 
forms as in the case of the built-up structure. 

a^(z,7) - [a^(z,z)]^ + [a^(z,z)]^ (2.5.1) 

““ M ^ P 

where [a^(z,z)] and [a^(z,z)] are given in Eqs. 2.1.1 and 2.2.1, 
respectively. The membrane thickness h(x) has to be replaced by t(x) in 
this case. Likewise, the load linear form can be expressed as 


(2.5.2) 
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— M — P 

where and [£.^(z)] are given in Eqs. 2.1.2 and 2.2.2, 

respectively. The state equation is 


a (z,z) = £ (z) (2.5.3) 

for all kinematically admissible virtual displacements z# 

Consider the functional representing the allowable stresses in the 
middle plane due to stretching and on the surface due to bending as 

^10 ^3 ^ ^6 (2.5.4) 


where and are given in Eqs. 2.1.21 and 2.2.23, respectively. 
Forces and moments of the component are [8] 


{n N N M M M y 

Xi XjX^ Xj x^ x^x^^ 



(2.5.5) 


where {a^} = [E] {e^^} are the membrane stresses and {<7g} = x^[E]{k} are 
the bending stresses. The curvatures {<} are 


{•c} 




(2.5.6) 


Membrane and bending stress resultants can be decoupled if the plate is 
symmetric with respect to the Xj^-X 2 plane. In EAL, moments and forces 
are related as 



(2.5.7) 


where [C] is the coupling coefficient matrix [2] 
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Design sensitivity analysis of the stress functional is chosen to 
illustrate the procedure. The sensitivity is the sunnnation of membrane 
and plate sensitivity as 





2 . . 2 

* // [- T o^^(z)e^^(X)]6t dfi + // [- y a^^(z)e^^(X)]6t dfi 

^ i,j = l n i,j=l 

+ // If- 5t dfi (2.5.8) 

which are given in Eqs. 2.1.25 and 2.2.27. Membrane stresses and 
strains are expressed in Eqs. 2.1.10 and 2.1.11, respectively and 
bending stresses and strains on the surface are expressed in Eqs. 2.2.15 
and 2.2.16, respectively. 

For the design sensitivity of the displacement functional, Eq. 

2.1.18 can be used for the in-plane displacements due to stretching and 

3 

Eq. 2.2.22 for the displacement z due to bending. 


2.5.2. Numerical Examples 

To demonstrate the numerical accuracy of the approach, a numerical 
example is tested. The finite element plate model is given in Fig. 2.6 
which is restrained at one side and loaded with a distributed tensile 
load in the positive Xj direction and nodal loads of 501b at nodes 21 , 
42, 63, 84 and 105 in the negative Xj direction as shown in Fig. 2.6. 

It contains 80 E43 elements (EAL bending plus membrane element type) , 
105 nodal points, and 500 degrees of freedom with the design variable 


►>< 
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thickness u = t(xj,X 2 )* The material properties are the same as the one 
in Section 4.1* 

Design sensitivity results of the von Mlses^ stress functional are 
given in Table 1 for Membrane, Bending and Total sensitivities with the 
perturbation of O.Olt. Since the stress resultants can be decoupled, 
the membrane sensitivity is expected to be the same as that of the 
original membrane model given in Table 7* However slightly different 
numerical values are obtained because numerical difference can occur in 
the decomposition process of membrane plus bending stiffness matrix when 
the process is compared to the membrane stiffness matrix alone. 



Figure 2. 6* Membrane plus Bending Plate Finite Element Model 

The same nodal points as in Table 6 are selected to check the 
accuracy of the design sensitivity of the displacement functional in the 
Xj, X 2 and negative x^ directions. Design sensitivity predictions and 
finite differences, with 6t = 0.05t are given in Table 2. 


Table !• Design Sensitivity Check for von Mises' Stress 
of Membrane plus Bending Plate with 6t = O.Olt 


El. 

No. 



iJ»(t+At) 


’I*' 

Ratio 

% 


M 

10053.27 

9953.71 

-99.56 

-100.53 

101.0 

1 

B 

12459.37 

12194.25 

-265.12 

-249.86 

94.3 


T 

22512.63 

22147.96 

-364.67 

-350.93 

96.2 


M 

9995.58 

9896.59 

-98.99 

-99.96 

101.0 

10 

B 

7897.13 

7732.11 

-165.01 

-158.86 

96.3 


T 

17892.71 

17628.71 

-264.00 

-258.81 

98.0 


M 

9999.86 

9900.97 

-98.89 

-100.00 

101.1 

20 

B 

748.51 

731.58 

-16.93 

-16.72 

98.9 


X 

10748.38 

10632.55 

-115.82 

-116.72 

100.8 


M 

8570.37 

8485.51 

-84.86 

-85.7 

101.0 

21 

B 

14344.51 

14040.40 

-304.05 

-289.13 

95.1 


T 

22914.89 

22525.91 

-388.89 

-374.84 

96.4 


M 

10019.71 

9920.49 

-99.21 

-100.19 

101.0 

30 

B 

7727.56 

7565.95 

-161.61 

-155.44 

96.2 


T 

17747.27 

17486.45 

-260.82 

-255.64 

98.0 


M 

9999.20 

9900.91 

-98.30 

-99.99 

101.7 

40 

B 

675.18 

642.77 

-14.41 

-13.94 

96.7 


T 

10657.08 

10543.67 

-113.41 

-113.93 

100.5 


M: membrane, B: bending, and T: total 
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Table 2. Design Sensitivity Check for Displacement of 
Membrane plus Bending Plate, 6 t = 0.05t 


Node 

No. 

Dir. 

Mt) 

i|;(t+At) 


r 

Ratio 

% 

63 

*1 

6.631E-3 

6.315E-3 

-3.157E-4 

-3.316E-4 

105.0 

63 

’^3 

-2.585E-1 

-2.214E-1 

3.710E-2 

3.895E-2 

105.0 

74 

*1 

3.297E-3 

3.141E-3 

-1.570E-4 

-1.649E-4 

105.0 

74 

X2 

-2.011E-4 

-1.915E-4 

0.961E-5 

1.004E-5 

105.0 

74 

X 3 

-7.995E-2 

-6.842E-2 

1.153E-2 

1.205E-2 

104.5 

105 

*1 

6.631E-3 

6.315E-3 

-3.157E-4 

-3.316E-4 

105.0 

105 

X2 

-3.996E-4 

-3.804E-4 

1.914E-5 

1.995E-5 

104.2 

105 

^3 

-2.582E-1 

-2.211E-1 

3.707E-2 

3.895E-2 

105.1 


Since the displacements can be decoupled in membrane plus bending 
elements, the sensitivities for Xj^ and X 2 directions are the same as 
that of Table 6 . The sensitivity for X 2 direction displacement at node 
63 is not considered since the displacement is zero due to the symmetry 
of loads and structure. 

2.6. Design Sensitivity Analysis of Pointwise Stress Functional 

2.6.1. Membranes 

Consider the general functional that represents a locally averaged 
stress on an element as in Eq. 2.1.21. 

= // g(^(z))m dn 


( 2 . 6 . 1 ) 
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where m is defined on a finite element in Eq. 2.1.22. If we have a 
P P 

smooth problem so the stress is continuous, we can consider polntwlse 
constraint. In this case, by letting the test area shrink to a point 

A 

X, m becomes a Dirac delta measure. Thus 
P 


’I'll * g(<J(z(x))) 

“ // ))<S(x - x)dft 

a 

The first variation of Eq. 2.6.2 is 


(2.6.2) 




2 

I 

i.j=l 


3o J 


// 

a 


2 

y 

i.j=l 


o^^(z')6(x - x)d£l 


(2.6.3) 


Replacing the variation in state z' by a virtual displacement X, the 
adjoint equation is obtained as 


aJX,X) 


If I M£l a^J(X)6(x - x)dO 
ft i,j=l 3a^J 


(2.6.4) 


for all kinematically admissible virtual displacements X. Using the 
adjoint variable method, the design sensitivity of the polntwlse stress 
functional is 




-// [ I a^^(z)e^j(X)]6h dft 

i,j=l 


(2.6.5) 


which can be computed using information at Gauss points as 


ni 


N 

I 


k=l 


y [- I (o^j(z)^(e^^(X))^lw^J^ 
£=1 i,j-l 


k 


6h 


k 


( 2 . 6 . 6 ) 
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where N is the total number of elements, I is the counter for the number 
of Gauss points, k is the counter for the element number, J is the 
Jacobian, and W is the weighting constant for the ilth Gauss point. 

If von Mises* stress functional is selected as the constraint 
functional, the formulations for stresses are given in Section 2.1. 


2.6.2 Plates 

Consider the functional that represents a locally averaged stress 
on the surface of the plate as in Eq. 2.2.23. If we consider the 
pointwlse stress, the stress functional and its first variation will 
become 


’*'12 “ g['^(z(x))) 


and 


f f g(o(z))6(x - x)dfi 

h 


'*'12“ ^ a^^(z»(x)) 

i,j=l 


* // I a^^(z’)6(x - x)d£2 

n i,j-l 

The adjoint equation can be defined as 


(2.6.7) 


( 2 . 6 . 8 ) 


,(x.« - I 

i,j*l 9a ^ 

= // I a^j(I)6(x - x)dn 

Q i,j=l 3a ^ 


(2.6.9) 


Then the design sensitivity of the pointwise stress functional can be 
obtained using the adjoint variable method as 
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r. 


12 


2 . . 

// “ [ I a^^(z)e^J(X) ] 6 t dft + // -^ (a(z)) 5 (x - x)dfl ... 
£1 l.j®! 

( 2 . 6 . 10 ) 


which can be computed as 


A* 

U2 


N 

I 

k=l 


I [- I + I 

i»l i,j=l i,j=l 


3gfg^^(z)) 


3t 

( 2 . 6 . 11 ) 


6(x - xM 


6t' 


where N, i, and k are explained in Section 2.6.1. In case of von Mises' 
stress constraint, the formulations of stresses for plate are given in 
Section 2.2. 


2.6.3. Numerical Examples 

The examples given in Sections 4.1 and 4.2 are used to test the 
accuracy of design sensitivity of the pointwise stress constraint. For 
the pointwise stress sensitivity, the first Gauss point shown in 

A 

Fig. 2.7 is selected as x. 


i 


1 

GP4 

• 

GP3 

• 


GPI 

• 

GP2 

• 



\*- 1 1 “H 



Figure 2.7. 2-by-2 Gauss Points 


To check sensitivity of the stress constraint of Eq. 2.6.2, the 
equivalent nodal forces of the adjoint load are computed using right 
side of Eq. 2.6.4. Design sensitivity analysis results for the von 
Mises' pointwise stress are given in Table 3 for several finite elements 
with the perturbation of 6t = O.Olt. 
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Accuracy of the sensitivity analysis results of the pointwise 
stresses at the first Gauss point of elements are almost equal to that 
of averaged stresses except at the first element. The design 
sensitivity results compared to the finite difference approximations are 
excellent. Especially the sensitivity of pointwise stress of the first 
element is better than that of the averaged stress in Table 7. 

The design sensitivities of the pointwise stress of plate problem 
are computed using Eq. 2.6.11. The adjoint load of Eq. 2.6.9 is used to 
get the equivalent nodal forces. Design sensitivity results of the von 
Mises' functional at the first Gauss point are given in Table 4 for 
several finite elements. The design perturbation for the finite 
difference computation is 6t = O.OOlt. The sensitivity results of 
pointwise stress are not quite excellent when compared to Table 10. It 
is found that the bending stresses are not evenly distributed on the 
plate surface for a given finite element element. That is, values of 
stress at four Gauss points are rather different from the average value. 

Table 3. Membrane Design Sensitivity Check for Pointwise 
Stress at the first Gauss point, 6t = O.Olt 


El. 


No. 

Kt) 

i|»(t+At) 

1 

10496.95 

10392.99 

10 

9990.33 

9891.39 

20 

9999.99 

9900.96 

21 

8453.86 

8370.15 

30 

10015.01 

9915.84 

40 

9999.93 

9900.94 


Ratio 

Ai|) t|»' % 


-103.95 

-104.97 

101.0 

-98.94 

-99.90 

101.0 

-99.03 

-100.0 

101.0 

-83.71 

-84.54 

101.0 

-99.17 

-100.15 

101.0 

-98.99 

-99.99 

101.0 
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Table 4. Plate Design Sensitivity Check for Pointwise 
Stress at the first Gauss point, 6t = 0.00 It 


El. 

No. 


i(i(t+At) 


r 

Ratio 

% 

1 

35.64 

35.57 

-0.071 

-0.074 

104.7 

2 

830.29 

828.62 

-1.661 

-1.665 

100.2 

3 

1841.86 

1838.17 

-3.691 

-3.609 

97.8 

4 

2716.80 

2711.37 

-5.432 

-5.331 

98.1 

5 

3252.14 

3245.63 

-6.510 

-6.374 

97.9 

7 

1157.15 

1154.83 

-2.314 

-2.376 

102.6 

8 

1273.81 

1271.26 

-2.553 

-2.473 

95.7 

9 

1169.60 

1167.25 

-2.341 

-2.112 

90.2 

10 

1038.10 

1036.02 

-2.078 

-1.735 

83.5 

13 

1310.30 

1307.68 

-2.620 

-2.670 

101.9 

14 

1151.13 

1148.83 

-2.307 

-2.432 

105.4 

15 

901.55 

899.75 

-1.805 

-1.985 

109.9 

19 

1402.82 

1400.01 

-2.813 

-3.057 

108.7 

20 

1533.27 

1530.20 

-3.069 

-3.291 

107.0 

25 

1883.25 

1879.47 

-3.77 

-3.901 

103.5 
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CHAPTER III 
PROGRAMMING ASPECTS 

So far analytical results and numerical algorithms for design 
sensitivity analysis have been stated. This Chapter outlines the basic 
organization of the EAL database management system and a structural 
design sensitivity analysis program that has been implemented using EAL. 

3.1. The EAL Database Management System 

Design sensitivity analysis of structural components and built-up 
structures can be implemented using a database management system. A 
survey of database management systems can be found in Ref. 12. EAL can 
be considered as a database management system [12], as well as a higher 
order programming language with advanced programming concepts [13]. 

EAL is a set of independent processors that communicate with each 
other through a random access database (Fig. 3.1). The database is 
manipulated according to user commands. The commands can be combined in 
a runstream, which can be stored in the database as a runstream 
dataset. Runstream datasets are driven by the Execution Control System 
(ECS), which also allows branching and looping within a runstream 
dataset. The ECS allows the user to call a runstream dataset from the 
database within a called runstream dataset. The ECS also initiates the 


execution of a new processor 





ECS-Executlon Control 
System 


EAL Processor 

one in central memory at a time 


Working Storage 
Dynamically Allocated 



Data Handling Routines 

I 



Register Set 
(In central 
memory all 
the time) 


Data Base 





Library 1 


Library 2 





Data set 1 


Data set 1 


Data set 1 



Data set 2 


Data set 2 
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• 

• . ♦ 




• 


• 
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Runstream 1 



• 
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• 


• 


. 



• 


• 


. 




1 


Figure 3.1 Data Flow in EAL 
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, Besides the regular database, EAL has a set of registers, which are 
stored in core. The user can assign a name, a type code (real, integer, 
or alphanumeric), and a value to each register. The contents of a 
register can be manipulated by register action commands. The registers, 
together with the ECS, enable controlled branching to be performed in 
runstream datasets. Register action commands have a higher priority 
than other commands, so whenever a register action command is given, the 
regular command procedure will be interrupted and the register command 
will be executed first. 

To manage the database, a set of database utility functions are 
provided. The database consists of one or more libraries, where each 
library contains a set of named datasets. The contents of each library 
are stored in a separate table, which for each dataset stores the name, 
data type, block size, number of columns, total number of words in a 
dataset, and location in the library. If many datasets are stored in 
one library, database overhead disk I/O can be very large. Database 
utility functions allow the user to change information in the table of 
contents, copy datasets from one library to another, and erase 
information in libraries. 

To manipulate information within each library, EAL provides the 
Arithmetic Utility System (AUS). AUS allows a wide range of vector and 
matrix manipulations as well as storing new datasets that are specified 
by the user. Matrices can be stored as full matrices or in a sparse 


storage forum 
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Besides the database management system, EAL contains a variety of 
processors common to finite element analysis programs. Although EAL's 
element library Is limited. It contains one, two, and three dimensional 
elements. The program Is able to solve static as well as buckling 
problems. The dynamic analysis part can handle eigenvalue and 
eigenvector solution and forced dynamic response analysis. 

Substructuring and graphical pre- and post-processing are also 
available. Because the global system matrices (stiffness, mass, and 
geometrical stiffness matrices) are usually very large and sparse, a 
sparse storage technique for hypermatrices Is used for global 
matrices. Thus, the user does not have to worry about bandwidth (as In 
SAP IV [14]) or wave front length (as In ANSYS [5]). For the solution 
of a large system of equations, the global stiffness matrix Is factored 
according to Ref. 15. 

EAL also allows the user to write his own program and combine It 
with the EAL database system. Because the user can use Fortran callable 
data handling routines , It Is relatively easy to create new processors 
for the EAL database management system [16]. Writing a separate 
processor for a specific task Is time consuming, but the result will be 
a more efficient system than just using the registers and the arithmetic 
utility processor. No separate processor Is written In this work, 
because the goal Is to show the feasibility of combining the design 
sensitivity analysis method with the EAL database management system. 

3.2. Program Organization 

Using the Engineering Analysis Language EAL, a general purpose 
design sensitivity analysis program has been written. The program can 
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handle three types of constraint functionals; compliance, displacement, 
and stress. Three types of elements (element types E21, E41, and E42) 
can be used to model a structure and evaluate design sensitivities. A 
flow chart of the program Is given In Figure 3.2. To use the program, 
the user sets the system control parameter, gives Information about the 
design variables, specifies the constraints, sets up his finite element 
model. The finite element model Is described In the runstream dataset 
INIT MODL 0 0, where all other Information Is given In the runstream 
dataset PARA SET 0 0. After that, the program automatically computes 
sensitivities for the given constraints and design variables. System 
control parameters are as follows: 

NLST - Number of load cases 
LCAS - Actual load case 

NBTD - Number of Independent beam/ truss design variables 
NMDV - Number of Independent membrane design variables 
NPDV - Number of Independent plate design variables 
NDV - Total number of design variables 

DE21 - Number of beam/truss elements 
DE41 - Number of membrane elements 
DE42 - Number of plate elements 
DETO - Total number of elements 

For simplification of the checking process only, the program allows 
only one control parameter per finite element. For beam elements, the 
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user can chose the height, the width, or a combination of both as 
independent design parameter. Weighting values for beam design 
parameter are introduced for selecting the height or width or any 
combination as design variables. The weighting values for beam design 
parameters are given in the table DESV VALU 0 0. 

The relationship between elements and design parameters are defined 
in the lists given below. For each element, there is one entry in the 
appropriate table that gives the design variable group. 

ED21 REL 0 0 Relationship for beam/truss elements 

ED41 REL 0 0 Relationship for membrane elements 

ED42 REL 0 0 Relationship for plate elements 

The constraint control parameters are 

CCOM - Compliance constraint 

GDIS - Number of displacement constraints 

CS21 - Number of stress constraints for element typ E21 

CS41 - Number of stress constraints for element typ E41 

CS42 - Number of stress constraints for element typ E42 

CTOT - Total number of constraints 

For every constraint group, with the exception of the compliance 
constraint, a table is required to describe the location of the 
constraint. For displacement constraints, two entries are needed for 
each constraint. The first entry is the node number on which the 
displacement sensitivity is evaluated and the second entry is the 
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direction of the displacement constraint. For each stress constraint 
group (E21, E41, and E42) a separate table is needed that gives the 
element numbers on which the stress constraint functionals are 
evaluated. The tables are 

ST21 LIST 00 - Gives the constraints for beam/ truss elements 

ST41 LIST 00 - Gives the constraints for membrane elements 

ST43 LIST 00 - Gives the constraints for plate elements 

For the stress constraint in a beam element (E21), the maximal 
stress is at one of the four corners of the beam cross section. The 
program computes stresses at four corners, finds the corner with the 
largest absolute value of stress , and computes sensitivity of the 
maximum stress. 

Additional parameters are; 

IDGP=1 - Pointwise stress at Gauss point case 
IDGP=0 - Averaged stress case 

(1) For E43 elements, 

DE43 - Number of E43 plate elements 
ED43 REL 0 0 Relationship for E43 plate elements 
CS43 - Number of stress constraints for element type E43 
For results, 

CIND has an additional number between 60,000 and 70,000 
indicating the stress constraint of E43 element. 

(2) For pointwise stress constraint, 

GPNO LIST 0 0 Indicates the Gauss point at which pointwise 
stress will be evaluated 
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For results, 

DVAL E42 1 1 The value of pointwise stress constraint functional 
at Gauss point listed in input data set of E42 

DVAL E41 1 1 The value of pointwise stress constraint functional 
at Gauss point listed in input data set of E41 
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CHAPTER IV 
NUMERICAL EXAMPLES 

In order to check whether the design sensitivity Information 
obtained is accurate, a comparison is made with the finite 
difference An appropriate design perturbation 4u must be selected, 

in order to obtain a meaningful finite difference of the constraint 
functional. That is, if 4u is too small, 4t|» = t|<(u + 4u) - 'l'(u) may be 
inaccurate, due to loss of significant digits in the difference. On the 
other hand, if ^u is too large, ^’1' will contain nonlinear terms and the 
comparison with ’I'’ will be meaningless. The design sensitivity 'P* of a 
constraint functional is the scalar product of the design sensitivity 
vector and the design variable perturbation vector ou. That is, the 
perturbation of an element design parameter is multiplied by the 
corresponding sensitivity component and the sum of all products is the 
design sensitivity ’I'*. 


4.1. Membrane 

The finite element membrane model in Fig. 4.1 is a simple plane 
elastic solid that is restrained at one end and loaded with a 
distributed tensile load at the other end. It contains 80 isoparametric 
elements (EAL plane stress element, type E41), 105 nodal points, and 200 
degrees of freedom, with the design variable thickness u = h(x). 

Young's modulus and Poisson's ratio are given as E = 3x10^ psi and 
Y = 0.3, respectively. Each finite element has uniform thickness. 
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Figure 4.1 Plane Elastic Solid Finite Element Mode 
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so that a maximum possible number of design variables is 80. For 
simplicity, a uniform thickness of h * 0.5 in. is used for the 
sensitivity check. 

Compliance sensitivity results are shown in Table 5, where 
s t|;(h + Ah) - ^^(h) and is the predicted value computed from 
Eq. 2.1.9, with design perturbations of O.Olh and ^h =* 0.05h. The 

percent accuracy of the sensitivity prediction is computed using the 
ratio X 100/^4^. 


Table 5. Membrane Design Sensitivity 
Check for Compliance 


5h 

t(h) 

'l'(h+4h) 

Aijj 


Ratio % 

O.Olh 

265.24 

262.62 

-2.63 

-2.65 

101.0% 

0.05h 

265.24 

252.62 

-12.63 

-13.26 

105.0% 

Several 

discrete 

nodal points 

shown in 

Fig. 4.1 

are selected 


check accuracy of design sensitivity of the displacement functional of 
Eq. 2.1.19. In order to compute this equation, the adjoint 
strain due to the adjoint load is needed. For each direction of 

displacement on a node, there is a separate sensitivity calculation that 
produces an adjoint strain Design sensitivity predictions and 

differences, with ^h = 0.05h, are given in Table 6. 
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Table 6. Membrane Design Sensitivity 
Check for Displacement 


Node 

No. 

Dir. 

Kh) 

Kh+4h) 

4<|> 

r 

Ratio 

% 

74 

X 

3.297E-3 

3.141E-3 

-1.570E-4 

-1.649E-4 

105.0 

74 

Y 

-2.011E-4 

-1.915E-4 

0.961E-5 

1 .004E-5 

105.0 

105 

X 

6.631E-3 

6.315E-3 

-3.157E-4 

-3.315E-4 

105.0 

105 

Y 

“3.996E-4 

-3.804E-4 

1.914E-5 

1.995E-5 

104.2 

63 

X 

6.631E-3 

6.315E-3 

-3.157E-4 

-3.316E-4 

105.0 


To check the stress constraint sensitivity of Eq. 2.1.27, the 
equivalent nodal force of the adjoint load on the right of Eq. 2.1.37 is 
computed for the finite element adjoint analysis and is obtained 

for each constrained element. Design sensitivity results for von Mises’ 
stress functionals are given in Table 7, for several finite elements. 
Perturbations are ^h = O.Olh and = 0.05h, for the von Mlses’ stress 
criteria. 

With all three constraint functionals, design sensitivity results 
compared to the finite difference approximation are excellent. It is 
interesting to note that in Tables 5, 6, and 7, the finite difference 
approximation is nearly 1% of the constraint functional when = O.Olh 
and nearly 5% of the constraint functional when <Sh = 0.05h. These 
results also show that as approaches zero, the ratio approaches 


one. 



Table 7. Membrane Design Sensitivity Check for Stress 
(a) von Mises’ Stress with * O.Olh 


El. 

No. 

"Kb) 

Kh+4h) 

Ai|) 

ft 

Ratio 

% 

1 

10053.402 

9953.863 

-99.539 

-106.804 

107.3 

10 

9995.646 

9896.680 

-98.966 

-99.957 

101.0 

20 

10000.141 

9901.130 

-99.011 

-99.999 

101.0 

21 

8570.358 

8485.503 

-84.855 

-86.209 

101.0 

30 

10019.743 

9920.537 

-99.206 

-100.197 

101.0 

40 

10000.065 

9901.054 

-99.011 

-99.996 

100.9 



(b) von 

Mises’ Stress with 

6h = 0.05h 


El. 

<Kh) 

Kh+4h) 


Ratio 

No. 




% 

1 

10053.402 

9574.668 -478.734 

-534.021 

111.5 

10 

9995.646 

9519.663 -475.983 

-499.958 

105.0 

20 

10000.141 

9523.944 -476.147 

-499.996 

105.0 

21 

8570.358 

8162.246 -408.112 

-431.046 

105.6 

30 

10019.743 

9542.612 -477.131 

-500.983 

105.0 

40 

10000.065 

9523.871 -476.194 

-499.982 

105.0 
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4,2 Bending of Plates 

The clamped plate element model shown In Fig, 4.2 is uniformly 
loaded with a pressure f(x) ■“ -1.5 Ib/in in the z direction. Since the 
model is symmetric with respect to the center and symmetric boundary 
conditions are applied, only one quarter of the plate is analyzed. The 
quarter model contains 25 4-node quadrilateral thin plate elements of 
type E42. It has 36 nodal points and 85 degrees of freedom. 

The design variable is the plate thickness u = t(x). Young's 
Modulus and Poisson's ratio are E = 30.5x10^ psi and Y = 0.3, 
respectively. The constant plate thickness is t = 0.4 in. Self-weight 
of the plate is neglected. 

Compliance sensitivity results are shown in Table 8, where 
4*1' <|'(t + 4t) - ’l'(t) and ’I*' is the predicted value that is computed from 
Eq. 2.2.11, with design perturbations ^t = O.Olt and 5t = 0.05t. The 1% 
design perturbation for the compliance constraint functional gives good 
correlation between design sensitivity and the finite difference 
approximation. However the 5% perturbation shows nonlinearity in the 
compliance of the plate element. 

Table 8. Plate Design Sensitivity 
Check for Compliance 


4t 

l'(t) 

Kt+4t) 

Alp 

'I'' 

Ratio % 

O.Olt 

5.0016 

4.8545 

-0.1471 

-0.1501 

102.0 

0.05t 

5.0016 

4.3205 

-0.6811 

-0.7502 

110.2 
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Several discrete points in Fig. 4.2 are selected to check design 
sensitivity accuracy for the displacement functional in Eq. 2.2.22. 

In order to compute this equation, as in the membrane case, the adjoint 
strain due to the adjoint load is needed. Some displacement 

results are shown in Table 9 for a design perturbation of ^t = O.Olt. 
The design sensitivity results of Table 9 agree very well with the 
finite difference approximation. 

Table 9. Plate Design Sensitivity 
Check for Displacement 


Node 

Ht) 


Ai|) 

I'* 

Ratio 

No. 





% 


9 

1.9287E-3 

1.8720E-3 

-5.6720E-5 

-5.7861E-5 

102.0 

10 

3.0860E-3 

2.9952E-3 

-9.0760E-5 

-9.2579E-5 

102.0 

14 

1.9287E-3 

1.8720E-3 

-5.6720E-5 

-5.7861E-5 

102.0 

17 

1.1349E-2 

1.1015E-2 

-3.3380E-4 

-3.4046E-4 

102.0 

20 

3.0860E-3 

2.9952E-3 

-9.0760E-5 

-9.2579E-5 

102.0 

23 

1.8585E-2 

1.8038E-2 

-5.4660E-4 

-5.5753E-4 

102.0 

27 

1 . 1349E-2 

1.1015E-2 

-3.3380E-4 

-3.4046E-4 

102.0 

28 

1.8585E-2 

1.8038E-2 

-5.4660E-4 

-5.5753E-4 

102.0 

32 

4.1105E-3 

3.9896E-3 

-1.2088E-4 

-1.2331E-4 

102.0 

35 

2.5256E-2 

2.4514E-2 

-7.4280E-4 

-7.5767E-4 

102.0 

36 

2.7124E-2 

2.6327E-2 

-7.9780E-4 

-8.1372E-4 

102,0 
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To check sensitivity of the constraint functional of Eq, 2.2.29, 
the equivalent nodal forces of the adjoint load are computed with 
Eq. 2.2.33. Design sensitivity results for the von Mises’ stress 
functional are given in Table 10, for several different elements. The 
perturbation for the finite difference calculation is 5t = O.OOlt. Note 
that even though the equivalent nodal force calculation for the adjoint 
load is not consistent with the hybrid method, since a displacement 
shape function is used , design sensitivity accuracy is excellent . 

Table 10. Design Sensitivity Check for Stress 
Von Mises' Stress with 6t = O.OOlt 


El. 

No. 

’l'(t) 

i|»(t+At) 

Alp 

r 

Ratio 

% 

1 

451.13 

450.23 

-0.90 

-0.889 

98.7 

2 

1128.39 

1126.14 

-2.25 

-2.185 

97.1 

3 

1762.13 

1758.61 

-3.52 

-3.476 

98.7 

4 

2268.37 

2263.84 

-4.53 

-4.518 

99.7 

5 

2549.64 

2544.54 

-5.10 

-5.097 

99.9 

7 

1253.26 

1250.76 

-2.50 

-2.543 

101.7 

8 

1248.91 

1246.42 

-2.49 

-2.471 

99.2 

9 

1034.00 

1031.93 

-2.07 

-1.967 

95.0 

10 

809.66 

808.04 

-1.62 

-1.460 

90.5 

13 

1288.46 

1285.89 

-2.57 

-2.548 

99.1 

14 

1210.50 

1208.09 

-2.41 

-2.380 

98.8 

15 

1084.12 

1081.96 

-2.16 

-2.134 

98.8 
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Table 10 — continued 


19 

1534.12 

1531.06 

-3.06 

-3.057 

99.9 

20 

1680.28 

1676.92 

-3.36 

-3.358 

99.9 

25 

1963.52 

1959.61 

-3.91 

-3.927 

100.4 


4.3 Beams 

A cantilever beam finite element model shown in Fig. 4.3 is loaded 
with a force F = [0.0 10.0 -10.0] lb. at the tip. It contains 20 2-node 

beam elements of type E21 and 21 nodal points with six degrees of 

freedom each. The beam has a rectangular cross-section with constant 
width and height, b = 0.5 in. and h = 0.25 in., respectively. Young's 

modulus and Poisson's ratio are E = 30.5x10^ psi and Y = 0.3, 

respectively. Self weight is excluded in the analysis. 

Compliance sensitivity results are shown in Table 11, where 
Alp s ip(u + Au) -ip(u) and *P' is the predicted value calculated from 
Eq. 2.3.9, with design perturbations iSfe = 0.01b and 6h = O.Olh. 

Table 11. Beam Design Sensitivity 
Check for Compliance 


iSh iSb ip(u) ip(u+4u) A<P <P' Ratio % 


O.Olh 0.01b 22.524 21.646 


-0.8789 


-0.9010 102.5 
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Figure 4«3« Beam Finite Element Model 

Several discrete points along the beam are selected to check 
accuracy of design sensitivity of the displacement functional of 
Eq* 2.3 • 14# In order to compute Eq. 2 .3 .14, the beam curvature of the 
adjoint displacement field is needed* Displacement results are shown in 
Table 12 for design perturbations of » 0.01b and ^h = O.Olh. Results 
show that design sensitivity predictions are close to the finite 
difference approximation. 
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Table 12. Beam Design Sensitivity 
Check for Displacement 


Node 

No. 

Direc 

tion 

Ku) 

4^(u+^u) 


■P' 

Ratio 

% 

3 

3 

-0.0272 

-0.0261 

0.00106 

0.00109 

102.6 

6 

3 

-0.1613 

-0.1550 

0.00629 

0.00645 

102.6 

9 

3 

-0.3904 

-0.3752 

0.01523 

0.01562 

102.5 

12 

3 

-0.6955 

-0.6684 

0.02713 

0.02783 

102.5 

15 

3 

-1.0577 

-1.0164 

0.04130 

0.04231 

102.4 

18 

3 

-1.4578 

-1.4009 

0.05690 

0.05831 

102.5 

21 

3 

-1.8769 

-1.8037 

0.07320 

0.07508 

102.6 

3 

2 

0.00546 

0.00525 

-0.2126E-3 

-0.2192E-3 

103.1 

6 

2 

0.03229 

0.03104 

-0.001259 

-0.001294 

102.8 

9 

2 

0.07814 

0.07509 

-0.003046 

-0.003129 

102.7 

12 

2 

0.13918 

0.13375 

-0.005430 

-0.005572 

102.6 

15 

2 

0.21162 

0.20337 

-0.008250 

-0.008470 

102.6 

18 

2 

0.29167 

0.28030 

-0.011370 

-0.011670 

102.6 

21 

2 

0.37551 

0.36087 

-0.014640 

-0.015028 

102.6 


To check stress constraint sensitivity of Eq. 2.3*21, the 
equivalent nodal force for the adjoint load and the right side of 
Eq* 2*3*18 must be computed for finite element adjoint analysis and the 
adjoint displacement field ^ must be obtained for each constraint 
element. Stress results for several finite elements are shown in Table 
13 for design perturbations * O.Olh and = 0.01b* 
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The reduced accuracy in correlation between the finite difference 
calculation and the first variation for the element number 20 results 
from an inaccurate finite difference calculation. 

Table 13. Beam Design Sensitivity Check for Stress 


El. 

No. 

Fiber 

t(u) 

iPCu+Au) 

A4* 

,j,t 

Ratio 

% 

1 

4 

60613.5 

58831.9 

-1781.6 

-1818.0 

102.0 

5 

4 

48187.0 

46770.9 

-1416.1 

-1445.2 

102.0 

10 

4 

32654.9 

31695.6 

-959.3 

-979.1 

102.1 

15 

4 

17121.3 

16621.7 

-499.6 

-512.8 

102.6 

20 

4 

1592.9 

1541.7 

-51.2 

-62.3 

121.8 


4.4 Built-Up Structures 

A built-up structure that uses beams and plates is shown in Fig. 4.4. 
The structure is clamped on two edges, with symmetric boundary 
conditions applied along the other two edges. A uniform pressure 
f(x) =1.5 Ib/in is applied on the top surface of the plate. The model 
contains 20 beam elements (E21) and 25 plate elements (E42). Beam width 
and beam height are two dependent design variables of the first 
independent design parameter, and the plate thickness is the second 
independent design parameter. The initial values are b = 0.05 in., 
h = 0.40 in., and t = 0.1 in. Young's modulus and Poisson's ratio for 
the beams and plates are E = 30.5x10^ psi and 7 = 0.3, respectively. 




















61 


Self weight is neglected. An example of the design sensitivity 
vector is given in Appendix Al. 

Compliance sensitivity results are shown in Table 14, ^ere 
Ai|> = <l'(u+4u) - ^(u) and ’I'* is the predicted value computed from 
Eq. 2.4.4, with the design perturbations = 0.01b 
and ^h = O.Olh for the first control parameter and perturbation of 
fit =• 0.0 It for the second control parameter. 

Table 14. Built-Up Structure Design Sensitivity 
Check for Compliance 


Control 

Parameter 

'J'(u) 

4'(u+^u) 

Al); 

r 

Ratio 

% 

1 

110.084 

109.086 

-0.998 

-0.9893 

99.1 

2 

110.084 

107.560 

-2.524 

-2.5577 

101.3 


Several discrete points in Fig. 4.4 are selected to check design 
sensitivity accuracy for the displacement functional in Eq. 2.4.4. In 
order to compute this equation, just as in the case of single 
components, the adjoint strain £^^(1) is needed. Some displacement 
results are shown in Table 15 for the design perturbations fib = 0.01b 
and fih = O.Olh for the first control parameter and a perturbation of fit 
= 0.0 It for the second control parameter. 



Table 15. Built-Up Structure Design Sensitivity 
for Displacement 


Node 

No. 

Control 'l'(u) 

Parameter 

4>(u+4u) 



Ratio 

% 

8 

1 

-0.02115 

-0.02131 

0.000197 

0.000196 

99.6 

8 

2 

-0.02115 

-0.02102 

0.000492 

0.000498 

101.2 

9 

1 

-0.06217 

-0.06157 

0.000605 

0.000606 

100.1 

9 

2 

-0.06217 

-0.06078 

0.001394 

0.001420 

101.3 

10 

1 

-0.10086 

-0.09992 

0.000941 

0.000937 

99.6 

10 

2 

-0.10086 

-0.09857 

0.002295 

0.002324 

101.3 

11 

1 

-0.12491 

-0.12372 

0.001190 

0.001192 

100.2 

11 

2 

-0.12491 

-0.12209 

0.002820 

0.002855 

101.2 

12 

1 

-0.13473 

-0.13348 

0.001250 

0.001242 

99.3 

12 

2 

-0.13473 

-0.13166 

0.003070 

0.003111 

101.3 

15 

1 

-0.18204 

-0.18032 

0.001720 

0.001715 

99.7 

15 

2 

-0.18204 

-0.17792 

0.004120 

0.004174 

101.3 

16 

1 

-0.29650 

-0.29375 

0.002750 

0.002727 

99.1 

16 

2 

-0.29650 

-0.28974 

0.006760 

0.006845 

101.3 

17 

1 

-0.37298 

-0.36955 

0.003430 

0.003407 

99.3 

17 

2 

-0.37298 

-0.36447 

0.008510 

0.008626 

101.4 

18 

1 

-0.40049 

-0.39683 

0.003660 

0.003628 

99.1 

18 

2 

-0.40049 

-0.39133 

0.009160 

0.009285 

101.4 

22 

1 

-0.48626 

-0.48186 

0.004400 

0.004365 

99.2 

22 

2 

-0.48626 

-0.47511 

0.011150 

0.011300 

101.3 

23 

1 

-0.61524 

-0.60970 

0.005540 

0.005478 

98.9 
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Table 15 — continued 


23 

2 

-0.61524 

-0.60110 

0.014140 

0.014329 

101.3 

24 

1 

-0.66123 

-0.65533 

0.005900 

0.005840 

99.0 

24 

2 

-0.66123 

-0.64600 

0.015230 

0.015434 

101.3 

29 

1 

-0.78180 

-0.77486 

0.006940 

0.006861 

98.9 

29 

2 

-0.78180 

-0.76377 

0.018030 

0.018279 

101.4 

30 

1 

-0.84074 

-0.83332 

0.007420 

0.007334 

98.9 

30 

2 

-0.84074 

-0.82131 

0.019430 

0.019688 

101.3 

36 

1 

-0.90440 

-0.89647 

0.007930 

0.007836 

98.8 

36 

2 

-0.90440 

-0.88347 

0.020930 

0.021217 

101.3 


To check design sensitivity of the stress constraint functional for 
plate elements in the built-up structure, the equivalent nodal force for 
the adjoint load of each constraint must be computed. Design 
sensitivity results for the von Mlses' stress functional are given in 
Table 16, for several different elements. The perturbations *Sh » O.Olh 
and ^b = 0.01b are for finite difference calculation of the first 
control parameter and ^t = O.Olt is the perturbation of the second 
control parameter. Note that the equivalent nodal force calculation for 
the adjoint load is done with shape functions that are Inconsistent with 
the hybrid method, but this has no effect on accuracy of design 
sensitivity calculations. 



Table 16. Design Sensitivity Check for Plate Stress 
Von Mises’ Stress with 6t = O.Olt 


El. 

No. 

Control ^(u) 

Parameter 

<Ku+Au) 


’1" 

Ratio 

% 

1 

1 

3669.32 

3635.34 

-33.98 

-34.11 

100.4 

1 

2 

3669.32 

3621.67 

-47.65 

-47.95 

100.6 

2 

1 

9094.17 

9005.53 

-89.17 

-86.17 

96.6 

2 

2 

9094.17 

8979.13 

-115.04 

-111.05 

96.5 

3 

1 

14410.77 

14276.46 

-134.31 

-132.65 

98.8 

3 

2 

14410.77 

14223.77 

-187.00 

-184.91 

98.9 

4 

1 

18484.13 

18309.74 

-174.39 

-172.25 

98.8 

4 

2 

18484.13 

18245.95 

-238.18 

-238.23 

100.0 

5 

1 

20882.59 

20688.57 

-194.02 

-192.30 

99.1 

5 

2 

20882.59 

20611.15 

-271.44 

-272.81 

100.5 

7 

1 

10370.80 

10277.77 

-93.03 

-92.73 

99.4 

7 

2 

10370.80 

10233.58 

-137.22 

-140.06 

102.1 

8 

1 

10381.71 

10288.66 

-93.05 

-91.21 

98.0 

8 

2 

10381.71 

10244.29 

-137.42 

-137.10 

99.8 

9 

1 

8802.67 

8730.74 

-71.93 

-68.44 

. 95.1 

9 

2 

8802.67 

8681.11 

-121.56 

-115.74 

95.2 

10 

1 

6956.20 

6900.33 

-55.87 

-51.49 

92.2 

10 

2 

6956.20 

6859.40 

-96.80 

-87.96 

90.9 

13 

1 

10757.59 

10664.65 

-92.94 

-90.63 

97.5 

13 

2 

10757.59 

10612.64 

-144.95 

.-143.94 

99.3 
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Table 16 — continued 


14 

1 

10079.36 

9990.00 

-89.36 

-87.01 

97.4 

14 

2 

10079.36 

9945.17 

-134.19 

-131.87 

98.3 

15 

1 

9011.52 

8931.79 

-79.73 

-77.53 

97.2 

15 

2 

9011.52 

8891.51 

-120.01 

-117.59 

98.0 

19 

1 

12794.24 

12682.81 

-111.43 

-109.85 

98.6 

19 

2 

12794.24 

12622.52 

-171.71 

-172.01 

100.2 

20 

1 

14073.98 

13952.86 

-121.12 

-119.39 

98.6 

20 

2 

14073.98 

13884.00 

-189.98 

-190.94 

100.5 

25 

1 

16610.25 

16472.22 

-138.03 

-135.45 

98.1 

25 

2 

16610.25 

16382.40 

-227.85 

-229.33 

100.6 


To check stress constraint sensitivity for beam elements in the 
built-up structure, the equivalent nodal force for the adjoint load on 
beam elements must be computed, so that the adjoint displacement field 1 
can be calculated. Stress results for several beam elements are shovm 
in Table 17. 

The forward finite difference method has been used so far to check 
the accuracy of the design sensitivity prediction. The prediction of 
the gradient for the first design parameter is very small for the beam 
element, which implies that the function has a nearly zero slope. For 
better finite difference approximation, the central finite difference 
method is used, to compare the accuracy of the prediction. The central 


finite difference method is defined as 



^u) ^ "■ ^u ) 

2 
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The perturbations are = 0.01b and ^h = O.Olh for the first control 
parameter ^u, and = O.Olt is the perturbation for second control 
parameter. 


Table 17. Design Sensitivity Check for Beam Stress 


El. 

No. 

Fiber 

Control 'l'(u) 

Parameter 

4^(u+^u) 

'l'(u-Au) 


r 

Ratio 

% 

1 

4 

1 

40182.0 

40191.2 

40161.1 

15.1 

17.8 

118.7 

1 

4 

2 

40182.0 

39277.9 

41109.3 

-915.7 

-915.7 

100.0 

2 

4 

1 

7430.8 

7457.8 

7401.8 

28.0 

27.9 

99.9 

2 

4 

2 

7430.8 

7245.7 

7621.7 

-188.0 

-188.0 

100.0 

3 

4 

1 

-9945.1 

-9947.6 

-9939.7 

-3.9 

-4.6 

117.4 

3 

4 

2 

-9945.1 

-9721.0 

-10174.6 

226.8 

226.8 

100.0 

4 

4 

1 

-17328.3 

-17337.9 

-17313.5 

-12.2 

-13.1 

107.4 

4 

4 

2 

-17328.3 

-16934.5 

-17732.3 

398.9 

399.0 

100.0 

5 

4 

1 

-20339.1 

-20363.2 

-20309.4 

-26.9 

-28.1 

104.5 

5 

4 

2 

-20339.1 

-19867.9 

-20823.6 

477.9 

477.8 

100.0 

6 

4 

1 

82011.8 

82049.4 

81950.6 

49.4 

56.0 

113.4 

6 

4 

2 

82011.8 

80152.9 

83919.1 

-1883.1 

-1883.1 

100.0 

7 

4 

1 

17977.1 

18035.4 

17914.1 

60.6 

62.9 

103.8 

7 

4 

2 

17977.1 

17534.3 

18433.3 

-449.5 

-449.6 

100.0 

8 

4 

1 

-18646.6 

-18639.5 

-18648.1 

4.3 

1.7 

40.0 

8 

4 

2 

-18646.6 

-18234.8 

-19067.7 

416.5 

416.5 

100.0 



Table 17 — continued 


■36874.1 -36910.3 -36648.1 
•36874.1 -36024.9 -37746.3 
44467.5 -44534.1 -44388.5 
■44467.5 -43426.8 -45537.6 


-41.5 

-43.5 

104.7 

860.7 

860.8 

100.0 

-72.8 

CNI 

• 

1 

106.1 

1055.4 

1055.3 

100.0 
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CHAPTER V 
CONCLUSIONS 

Results of this study show that it is possible to combine the 
design sensitivity algorithms of Ref. 1 with the database management 
system of EAL. For stress constraints as performance criteria, it is 
necessary to compute equivalent nodal forces for the adjoint load. For 
plate elements, the EAL finite element analysis is based on a hybrid 
formulation, but a displacement finite element formulation is used for 
evaluating the equivalent adjoint nodal forces. Nevertheless, results 
of the design sensitivity analysis are very accurate, which indicates 
that it is not necessary to compute equivalent nodal forces for the 
adjoint load using exactly the same shape functions that are employed in 
finite element analysis. 

A database management system with a finite element capability and 
the adjoint variable method of design sensitivity analysis, permit 
implementation of a design sensitivity analysis method that does not 
require differentiation of element stiffness and mass matrices. It is 
shown that a database management system can be used to implement design 
sensitivity analysis, so only one program with one database is 
necessary. 

Work is progressing to extend the methods presented in this report 
to include shape (geometric) design parameters. A domain method [17] 
for shape design sensitivity analysis and a design component method [6] 
for sensitivity analysis of built-up structures are used for software 
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implementation. Numerical implementation and results of shape design 
sensitivity analysis will be reported in Part II: Shape Design 


Parameters 



70 


REFERENCES 

1. Haug, E.J., Choi, K.K., and Komkov, V., Structural Design 
Sensitivity Analysis , Academic Press, New York, N.Y., 1985. 

2. Whetstone, W.D., EISI-EAL Engineering Analysis Language. 

Reference Manual. EIS Inc., July 1983. 

3. Choi, K.K., Santos, L.T., Frederick, M.C., "Implementation of 
Design Sensitivity Analysis with Existing Finite Element Codes," 
Journal of Mechanics. Transmissions, and Automation in Design, 85- 
DET-77 . 

4. Frederick, M.C., and Choi, K.K., "Design Sensitivity Analysis with 
APPLICON IFAD Using the Adjoint Variable Method," Technical Report 
84-17, Center for Computer Aided Design, University of Iowa, 1984. 

5. Santos, J.L.T., and Choi, K.K., "Implementation of Design 

Sensitivity Analysis with ANSYS," Technical Report 85- , Center for 

Computer Aided Design, University of Iowa, 1985. 

6. Choi, K.K., and Seong, H.G., "Design Component Method for 
Sensitivity Analysis of Built-up Structures," Journal of Structural 
Mechanics , to appear, 1986. 

7. Plan, T.H.H., "Derivation of Element Stiffness Matrices by Assumed 
Stress Distribution," AIAA, July 1964, 

8. Cook, R.D., Concept and Application of the Finite Element Analysis , 
John Wiley and Sons Inc., New York, 1981. 

9. Zienswicz, 0.0., The Finite Element Method, McGraw-Hill, London 
1977. 

10. Washizu, K., Variational Methods in Elasticity and Plasticity , 
Pergamon Press, 1982 

11. Plan, T.H.H., "Finite Element Methods by Variational Principles with 
Relaxed Continuity Requirement,", in Variational Methods in 
Engineering, edited by Brebbia, C.A. , Tottenham, H. , Southampton 
University Press, 1973. 

12. Sreekanta Murthy, T., Arora, J.S., "A Survey of Database Management 
in Engineering," Advanced Engineering Software, Vol 7 (3), 1985. 

13. Prasad, B. , "An Integrated System for Optimal Structural Synthesis 
and Remodelling," Computers and Structures, Vol 20 (5), 1985. 



71 


14. Bathe, K.J., Wilson, E.L., Peterson, F.E., SAP IV; A Structural 
Analysis Program for Static and Dynamic Response of a System , 
College of Engineering, University of California, Berkley, 1974. 

15. Whetstone, W.D., "Computer Analysis of Large Linear Frames,", ASCE, 
Journal of the Structural Division. Vol 95, November 1969. 

16. Giles, G.L., Haftka, R.T., "SPAR Data Handling Routines," 

NASA Technical Memorandum 78701, September 1978. 

17. Choi, K.K. and Seong, H.G., "Domain Method For Shape Design 
Sensitivity Analysis of Built-Up Structures", Computer Method In 
Applied Mechanics and Engineering, to appear, 1986. 



72 


APPENDIX A1 - DESIGN SENSITIVITY VECTORS 

3ij/ 

This appendix lists the design sensitivity vectors for the 
built-up structure for the compliance constraint, for the displacement 
constraint at node 36, for the stress constraint in beam element 1, and 
for the stress constraint in plate element 25. 

Table 18. Sensitivity Vectors for the Compliance Constraint 


Beam element 




Plate element 



1 

2 

3 

4 

5 

6 
7 
,8 
9 

10 

11 

12 

13 


-11.62 

-30.97 

-0.62 

-1.65 

-0.73 

-1.94 

-1.97 

-5.26 

-2.72 

CM 

• 

1 

-47.55 

-126.80 

-3.13 

1 

00 

• 

u> 

-2.51 

-6.68 

-8.94 

-23.85 

-12.97 

-34.58 

-11.62 

-30.97 

-0.62 

-1.65 

-0.73 

-1.94 

-1.97 

-5.26 


1 

-9.83 

2 

-47.94 

3 

-127.21 

4 

-219.05 

5 

-282.67 

6 

-47.94 

7 

-50.04 

8 

-50.25 

9 

-38.58 

10 

-27.41 

11 

-127.20 

12 

-50.25 

13 

-59.09 

14 

-62.54 


14 
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Table 18 — continued 


15 

-2.72 

-7.26 

15 

-59.77 

16 

-47.55 

-126.80 

16 

-219.03 

17 

-3.13 

-8.34 

17 

-38.58 

18 

-2.51 

-6.68 

18 

-62.54 

19 

-8.94 

-23.85 

19 

-114.17 

20 

-12.97 

-34.58 

20 

-145.47 




21 

-282.65 




22 

-27.41 




23 

-59.77 




24 

-145.47 




25 

-202.79 

I 

-185.52 

-494.66 

I 

-2557.65 

«u 

0.004 

0.0005 

6u 

0.001 


-0.7421 

-0.2473 

ij,. 

-2.557 


In Table 18, the sum of the sensitivity components i®l,***,n 

is given in the third row from the bottom. When multiplying the sum of 
the sensitivity components with the perturbation ^ of the design 
variable u, one gets the first variation, which is given in the last row 
of Table 18. Note that this result coincides with results given in 
Table 14. Results in Tables 19, 20, and 21 coincide with results given 
in Tables 15, 16, and 17, respectively. 
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Table 19. Sensitivity Vectors for the Displacement Constraint 


1 element 

3ij; 

d\|) 

Plate element 


1 

0.0604 

0.1610 

1 

0.040 

2 

0.0074 

0.0197 

2 

0.251 

2 

0.0010 

0.0028 

3 

0.755 

4 

0.0144 

0.0386 

4 

1.426 

5 

0.0286 

0.0763 

5 

1.934 

6 

0.3205 

0.8547 

6 

0.251 

7 

0.0405 

0.1079 

7 

0.355 

8 

-0.0068 

-0.0181 

8 

0.429 

9 

0.0587 

0.1564 

9 

0.394 

10 

0.2098 

0.5595 

10 

0.311 

11 

0.0604 

0.1610 

11 

0.755 

12 

0.0074 

0,0197 

12 

0.429 

13 

0.0010 

0.0028 

13 

0.448 

14 

0.0144 

0.0386 

14 

0.457 

15 

0.0286 

0.0763 

15 

0.389 

16 

0.3205 

0.8547 

16 

1.426 

17 

0.0405 

0.1079 

17 

0.394 

18 

-0.0068 

-0.0181 

18 

0.457 

19 

0.0589 

0.1564 

19 

0.958 

20 

0.2098 

0.5594 

20 

1.490 


Table 19 — continued 





21 

1.934 




22 

0.311 




23 

0.389 




24 

1.490 




25 

3.744 

I 

1.469 

3.9176 

I 

21.217 

«u 

0.004 

0.0005 

6u 

0.001 

I'* 

0.005876 

0.0001958 

<!>• 

0.021217 
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Table 20. Sensitivity Vectors for the Beam Stress Constraint 


1 element 

3i|i 

■5b 

3\|^ 

Plate element 

dip 

■5T 

1 

77168 

-62100 

1 

-15002 

2 

9230 

24614 

2 

-52185 

2 

-652 

-1739 

3 

-80766 

4 

603 

1610 

4 

-87753 

5 

1108 

2955 

5 

-83260 

6 

-18728 

-49942 

6 

-88775 

7 

-97 

-258 

7 

-13790 

8 

-1728 

-4608 

8 

4890 

9 

-2162 

-5765 

9 

3402 

10 

-1163 

-3102 

10 

-1870 

11 

-13565 

-36174 

11 

-111076 

12 

618 

1648 

12 

23378 

13 

-2072 

-5524 

13 

-45076 

14 

587 

1565 

14 

-21198 

15 

1405 

3746 

15 

-4547 

16 

-15697 

-41858 

16 

-104176 

17 

-709 

-1889 

17 

1248 

18 

-1428 

-3807 

18 

-29999 

19 

-2440 

-6506 

19 

-30033 

20 

-1815 

-4841 

20 

-20950 
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Table 21. Sensitivity Vectors for the Plate Stress Constraint 


element 

9i|; 


Plate element 

9i|) 

X . 

1 

-641 

-1711 

1 

-361 

2 

-117 

-312 

2 

-2660 

2 

19 

51 

3 

-8478 

. 4 

-201 

-538 

4 

-16560 

5 

-444 

-1185 

5 

-22898 

6 

-3784 

-10090 

6 

-2660 

7 

-865 

-2307 

7 

-4650 

8 

1316 

3509 

8 

-6543 

9 

2632 

7018 

9 

-6805 

10 

-10611 

-28298 

10 

-5715 

11 

-641 

-1711 

11 

-8477 

12 

-117 

-312 

12 

-6543 

13 

19 

51 

13 

-7439 

14 

-201 

-538 

14 

-6114 

15 

-444 

-1185 

15 

4866 

16 

-3784 

-10090 

16 

-16560 

17 

-865 

-2307 

17 

-6805 

18 

1316 

3509 

18 

-6115 

19 

2632 

7018 

19 

-803 

20 

-10612 

-28299 

20 

2472 
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Table 21 — continued 





21 

-22897 




22 

• 

-5714 




23 

4864 



- 

24 

2472 




25 

-79195 

I 

■ -25392 

-67726 

I 

-229328 


0.004 

0.0005 

«u 

0.001 

'1^' 

-101.57 

-33.86 


-229.328 




APPENDIX A2 - SOURCE PROGRAM 


vJaUT Ul. 


^ 

rp 

^ 4 * Program far computing the first variation of stress# displacement 
and compliance constraints 

Written bg B. DOPKEF^ 


iii29 DR IV GLOB) 

‘i 

^^OCALL# (29 PAHA SET 0 O) 
-^OCAI L# (29 INIT hODL O O) 
^K)CAL.L# (29 DSA COM 0 0) 

^^KETURW 

it END DR IV 


FNDDHIV 




4.' 

><29 PARA 
ir4LST^l 
! L.CAS = 1 


SET 0 O) ENDPARA 

$ NUMBER OF LOAD CASES 
% NORMAL LOAD CASE 


•NBTD-=1 $ NUMBER OF 
•NMDV-=0 % NUMBER OF 
*NPDV-=1 % NUMBER OF 
’4» 

I NDV-^-NB TD+NMDV+NPDV 

,« 

•DE21=20 % NUMBER OF 

• 0E41---0 -4» NL'MBER OF 

*DE42-=25 % NUMBER OF 


INDEPENDENT BEAM/TRUSS DESIGN VARIABLES 
INDEPENDENT MEMBRANE DESIGN VARIABLES 
INDEPENDENT PLATE DESIGN VARIABLES 


BEAM/TRUSS ELEMENTS 
MEMBRANE ELEMENTS 
PLATE ELEMENTS 


' DET0“-^DE2H-DE4MDE42 
% 

^ DEFINE INITIAL VALUE FOR ALL DESIGN VARIABLES 




«XGT AUS 

TABLE(NI=:2# NJ^^-NBTD" ) : DESV VALU O 0 

J=l: 10 1.0 


^ DEFINE RELATION BETWEEN ELEMENTS AND DESIGN VARIABLES (SIZING) 

% BEAM/TRUSS 
>JZ(DE2i# 10) 

TABLE<NJ.-=n)E21"i TYPE^-^0): ED21 REL 0 0 
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I 

I 

i 

1 

> ORIGINAi: PAGE IS 

'I OF POOR QUAUTY 

i 


1 
i 
i 
I. 

1 
1 
1 
1 
i 
1 

M\EL U) 

JZ(D£*l 1, :30 ) 

TABI-t type; -^0) : 1 ED41 REL O O 

: 1 
J-k : I 

: i 

: X 

: 1 
J--6 ; 1 

; I 
J-=-S ; 1 

v;-9 : X 

J=10: 1 

J-11: 1 

J=12: 1 

J-13: 1 

J“- 14: 1 

1 

J^i6: 1 

t 

J^lB: I 

J---19: 1 

J-30; 1 

J-31: 1 

vR---c‘ai: i 

1 

J=24: 1 

I 1 

LAB EL 30 
PlAlL- 

-fi-jz ( or:4:>; 40> 

TABLE(M0-^“DE42", TYPE=0> : ED42 REL 0 0 


i 

■ .» 

J--3 

- » 

J -- 3 

;J 

J^4 






J^/ 

; » 


J -- 1 
j 

Jr.:j 

w '--4 

J^-6 
j- / 
J-L 

1 0 
1 

J-0 2 
J - i 3 
J- i 4 
1 ‘3 
.J -16 
J 3 / 
J - 1 13 
I 9 
J -30 
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J^8 : 

: :> 
vJ^=10: 2 

J. 

1.2: 2 
J*-13; 2 

J-=i4. 2 

2 

J-=^16: 2 

v^-'17: P. 

2 

J - 1 9 : 2 

J'-P.O: 2 

2 

J-^22; 2 

J-=^23: 2 

J-'=24: 2 

J=-25: 2 

:H, ABEL 40 


qt: pooh 


PAGS is 

quaint Y 


■i CON r K A I N r dsf i t' 4 1 r i on 


! CClJM =0 
' CDIS^O 
I CS21=^0 
? CS41 ^0 
‘ CS4 2'^1 


i COHPL I ANCE CONTRAINTS 

NUMBER OF DISPLACEMENT CONSTRAINTS 
'4> NUMBER OF STRESS CONSTRAINTS IN BEAM/TRUSS ELEMFZNTS 

^ NUMBER OF STRESS CQf4STRAINTS IN MEMBRANE ELEMENTS 

* NUMBER OF STRESS CONSTRAINTS IN PLATE ELEMEN fS 


? CGTR==CS2i+C34i‘*‘CS4P 


if 

‘ (riOr-'CCGM^ COIGfCSTR 


^TABLES FOR CONSTRAINT NODES/ELEMENTS 
i^xJZ<CDIS» 100) 

TAEL E(Nl:-2. NJ~"CDIS‘L TYPE=0) : DICO LIST 1 

: 36 3 ^ FIRST NUMBER 

X- J-2 : 9 3 

^ ; 10 3 

T J^4 : 11 3 

T : 12 3 

H LABEL 100 
>JZ<CS2l. 110) 

TaQL.E(NJ=“CS21 “4 TYPE==0> : ST21 LIST 
J=^a : 1 

^a..ADEL no 

% 


njlAQS^li 130) 

TABLE (NJ^^“CS41”/ TYPE==0): ST41 LIST 


J-^1 

1 

J^-2 


, 

10 


127 

y)=^^ 

1 


31 

A-- -^ 7 

61 

J-S 

64 


^a.AoeL 130 


«JZ(CS^^2. 140) 

TAGL.E(NJ-'*CS42n TYPE^O) : ST42 L.IST 


1 

IS THE NODE, 


SECOND IS DIRECTION 



- i : 

^H..ABF.L U!0 
»iUFTURW 

'A 

4. 


ENDPARA 


it 

-^^^29 lUlJ HOOL p O) END I ^4 IT 

■i 

aXfir TAB 
START 3A 

TITIC. ' res TING OF MEMBRANE Oh 80 ELEMENTES 
MATERIAL CONSTANTS 
:l 30500000. 0. 3 300. 0 

JOINT LOCArU'JNS 

I -20. O -20. 0 0. 0 0.0 —20. 0 O. 0 6 i i 

■i, -20 O O. O O. 0 O. O O. O O. O 

HREF 

FORMAl “2 

1,1, 30. O, 30. 0, 0. O 
H.A 

1 0 , 10 
B.T 

RLCT 1,0.40 0.05 

CONSTRAINT OEFINIT ION 1 
/E«D 1 2 6 


1 , 36 

Symmetry plane- i 

BVl'iMETRY PI .ANE-L’ 

ZERO 1 rz 3 4 5 6 
1. 6. 1 
/, 31,6 
iiXQ T EL.l) 

I.'; 2 1 

13 14 1 5 
25 26 1 5 
3 9 15 
5 11 1 5 
L42 

1 2 0 7 1 5 5 

aXOT E 
aXOT EKS 
aXOT IAN 
*XOT K 
aXOT DRSI 
ax or ,AOf) 

ALPH.A: CA.5E TITLE 

1 'UNIFORM ELEMENT PRESSURE 
TABLE; NODAL PRESSURE 
J--=l,36: -1.0 

-aXQT ENNF 
RESET SET--="LCAS'' 

»XOT SSOL 
RESET SET---“LCAS‘‘ 
aXOT VPRI 

PRINT 1 GTAT DISP "LCAS" 1 


ORK 


PAGE B 




0^7 


rjTY 


! TFMP^-OF.1 0-0E21 
iKiZn'FHP/ 200 > 
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ii-JZ(DE:41, 1^36) 

*;<QT ES 
PrlODE: -- 2 
QUTLIH^l 
OUTPUT TYPE^-ESU 
U- STAT DISP “LCAS" i 
l:4 1 1 

QUrPUT TYPr=^ESl) 

U=^ STAT DISP "LCAS" 1 
E4 1 1 

-.a ADEL i'Oh 
*JZ(DE42, 157) 

*XQT ES 
PhfJDE^2 
GUILIfi=^l 
OUTPUT TYPE-ESD 
U-== BFAT DISP “LCAS" 1 
L4:2 1 

OUTPUT TYPE-ESI) 

U= STAT DISP “LCAS" 1 
1 

*l-ADEl 15/ 
iiLADLI... 200 
! TEMP=FRLE( ) 

4> 

^ PhTURN 

^ i:md I N I r 
i> 

^(29 DSA COM O 0) ENDDSAC 

T 

% PREPARE DATA FOR LATER USAGE IN SENSITIVITY CALCULATION 
*JZ(0L2i, 40) 

j* DC ALL, (29 PREP £21 0 0) 

I ADEL 40 
ifJ2(DE41, 50) 

^^nCAI_L, (29 PREP E41 0 0> 

^l.ABLL 50 
->J2U)E42, 60) 

J^OCALL, (29 PREP E42 O 0) 

^t-I.ABEL. 60 

T LOOP OVER ALL CONSTRAINTS 

•4. 

!CLC-0 ^ REGISTER FOR CONSTRAINT COUNTING 

M. ABEL 500 
!CLC^CLC-» 1 

! adlc=ioook:lc 

■PSI=0.0 ^ REGISTER FOR VARIATION SUMMATION 

' VMST^O. 0 
U'UIFLU 
‘ CLC 

i^DCALL, (29 ADJQ COMP 0 0) 

">DCAIL, (29 AD JO SOLU O O) 

DC ALL, (29 SENS CALC O O) 


GRIcmAi: PAGE IS 
OF POOR QUALITY 
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ICLCC^^CLC- C TOT 
■JK JNZ ( CLCC. r>00) 

$ 

-rtKtVrURN 

‘4i t 

ii'C.NDDSAC 

^ 

‘4^ 

ii-(29 ADJa COMP O O) ENDAOCC) 

‘ V£hP-Ul .C-1 

' jr-:fxn =^0 

;^JGZ( TEhP/ B70) 

JZ (CCOM/ B/0) 

$ COMPUTE COMPLIANCE 
% 

i*-XGr DCU 

TOCC 1 ST AT 01 BP "LCAS" 1: N3^“ADLC" 

t 

>*JUMP 1000 

^LAUEL 870 
■t 

' TEPIP^Cl.C-CLUM- GDIS 
.»vJGZn‘El iP. 080) 

^ COMPUTE ADJOINT FOR DISPLACEMENTS 

i' 

‘ 1 EMP-^^CL.C-CCOn 

^GET TEMP'S ENTRY OF DICO LIST 


! xJENT^^DS 

1, "TEMP", 

1 ( 1 

DICO 

LIST i 

1 ) 

' IDIR^DS 

a; 

'-) • 
c. , 

•TEMP", 

1 { 1 

DICO 

LIST i 

1 ) 

(iXQT AUS 
SYSVeC ; 

APPl 

.. FORC: 

"ADLC " 

1 




IDIR** 

Jt'MT“ 

1 . 0 
T 

! TDIR^FREE C ) 

^sMJMP 1000 
-riLAUPL 8U0 

^ COMPUTE STRESS ADv-IOINTS 


% 

* FOR E21 STRSS CONSTRAINTS 
! VE MP- CLC COM- CD I S-C S21 
! Tl"MP 

’tJOZ (TEMP, 890) 

•f 

! vemp-^clc-ccom-cdis 

I VLMP 
‘ DE21 

:JiGEr TH:MP'S entry of ST21 LIST 
! JEN r=--OS ’'TEMP’L 1 ,1 (1 ST21 LIST 1 1) 


OmcmAL PAGE IS 
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* COMPUTE STRESS CONSTR FOR ELEMENT: JENT £21 

* 

a-DCALL, (29 ADJQ £21 0 0) 

» 

a JUMP 1000 

t- 

•iiL.MlEL 090 

V 

P FOR E4 1 STf^SS CONSTRAINTS 

! n:;MP^= ct.oocoM-CDis- cs?i-cs4i 

• VENP 

iiJOl ( TChP, 73 0) 

? f EMP - c Lc *-c c on - c D I s-c sr> 1 

! TEMP 

3»QET TEHP'S ENTRY OF ST41 LIST 
! Dii41 

‘JENT^DS "TEHP“/ 1 ,1 (X ST41 LIST 1 1> 

^ COMPUTE STRESS CONSTR FOR ELEMENT: JENT £41 

i|> . 

^DCALL. (29 AO JO E4i 0 O) 

« ,‘UMP 1000 

ia.ABEL 910 
I 1 t.MP^FRFF( ) 


t FOR E42 STRSS CONSTRAINTS 
? TEhP = CLC-CC0M-CDIS-CS21-CS41-CS42 
‘ TEMP 

^tJGZ(TENP, 1000) 

! TEhP=CLC-’CCOM-CDIS- C32i~CS41 
• TEMP 
! DE42 

$GET TEMP'S ENTRY OF ST42 LIST 
iJENT^DS ••TEMP'S! ,1 (1 ST42 LIST 1 1) 

% 

i COMPUTE STRESS CONSTR FOR ELEMENT: JENT E42 

% 

*DCALL, (29 ADJO E42 O 0) 

>^l_ABEL 1000 
! TEMP=-FREF( ) 

! JENT=^FREF( ) 

a RETURN 

iiENDADCO 

^ : 

'i 

h^(29 AO jo SOLU 0 0) 

$ 

! rEMP==(:LC--i 
JOZ< TEMP, 095) 


ENDAOSO 
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(ccoM, 

$ THIS JUMP IS ONLY FOR COMPLIANCE CONSTRAINTS 

>JUMP B86 

M-ABEL 800 
*:<GT SSQL 
RESEr SET=“ADL(:“ 

:|i 

i^XGT VPRT 

* 4 : 

PRINT 1 STAT OISP "ADLC“ 

^lY-^BEL 886 
•i 

! TFMP^DFTO 
( TEMP. J 08) 

•:LJZ(DL41, 106) 

>^XQT LS 

PHOrE^- 2 

OUTLIlL-^^1 

OUTPUT TYPE-ESH 

U-- STAT DISP "ADLC" 1 

L41 I 

^4 LABEL 106 
>JZ VDE42, 107) 

«XGT ES 
PhaDE->;? 

OUTLIB^l 
GUTPUl TYPE^ESH 
U= STAT DISP "ADLC“ 1 
E42 1 

*LABFL 107 
^LABEL 108 
' TEMP ^ FREE ( ) 

‘t 

■^RETURN 

* ENDAOSO 


^ 

^(2V SENS CALC O O) .ENDSECA 

‘ ICDV^ 1 

! i i::mf=-o 

^l,AUEL 300 
IPSI^O O 


•«JZ(D£21, 140) 

»DCALL, (29 SEGA E2t O 0) 
*LABEL 140 

>JZ(DE41. 100) 

^DCALL, (29 SEGA E41 O 0) 
•:t-LABEL 100 

^JZ(DE42, 160) 

^DCALL. (29 SEGA E42 O 0) 
« LABEL 160 
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i u:dv 
! 

‘ i;CDV^ICOVi 3. 

! IODV+ j. 

f ICDV 
‘ I t-MP 

aJGZ< TEMPi :.)00) 

! TEMP-FPErC ) 

? ICDV=^f-RL' P< ) 

‘■t 

iUiETURN 

I 

:iENDSECA 

'I 


ORIGINAE PAGE IS 
OE POOR QUALITY 


*(29 PREP fc;2l, O 0) ENDPE21 

•4 

S THIS SUBROUTINE IS FOR ELEMENT TYP E21 
» THE PRINCIPAL LIRRARY IS 13 

$ T'EMPORV LIBRARY IS 17. WILL BE DELETED AT THE END OF THIS ROUTINE 
* 


ij^XQT Eli 
EXTRACT 

S0URCE=E21 : 


CONTENT 

SPEC: 

GEOM 

1 % LENGTH OF ELEMENTS 

CREATE 13 

ii. 

LEWG E21 0 

0 





EXTRACT: 

B0URCE=E21: 


CONTENT 

SPEC: 

GEOM 

5, 13 DIRECTION COSINES 

CREATE 13 

DIRC £21 0 

0 





EXTRACT: 
CREATE 13 

EXTRACT: 

S0URCE=^E21: 
ELNQ E21 0 

0 

CONTENT 

SPEC: 

GEOM 

14,31 ^ ELEMENT-NODAL RELATION 

5GURCE-E21: 


CONTENT 

SPEC: 

MATE 

1,2 $ MATERIAL PROPERTIES 

CREATE 13 

MATP E2i 0 

0 





EXTRACT. 

S0URCE-E21: 


CONTENT 

SPEC: 

INTE 

1, 18 % CONNECTIVITY 

CREAi'E 13 

CONN E21 0 

0 





■p 

EXTRACT: 

S0URCE--=E21: 


CONTENT 

SPEC: 

SECTION 16,23 $ STRESS POINT LOCATION 

CREATE 13 

4: 

SPLO E21 0 

0 





J^XQT DCU 







TOCC 13 CONN E21 0 0 


TYPE^O 





* IF NOT ALL NODAL INFORMATION IS IN GLOBAL COORDINATES USE LTOG 
1» 

*XQT AUS 
! ICQU-0 
«LABEL 100 
! ICOU=ICOUi 1 

■P1=DS. 1.1, •• ICOU", 1 (13 CONN E21 O 0) 

•P2=DS, 14. "ICOU", 1 (13 CONN E21 O 0) 


UXl^OS. 1, “PI". 1 

(1 

5TAT 

DISP 

“LCAS“ 

1 ) 

UY1=DS, 2. “PI", 1 

( 1 

STAT 

DISP 

“LCAS“ 

1 ) 

UZ1=DS, 3, “Pl“, 1 

( 1 

STAT 

DISP 

“LCAS** 

1 ) 

RX1=^DS, 4, “Pl“, 1 

( 1 

STAT 

DISP 

“LCAS“ 

1 ) 

RY1=DS, 5, “PI", 1 

( 1 

STAT 

DISP 

“LCAS" 

1 ) 


I RZl^DS, 6i ‘ 

PI " / 1 

( 1 

STAT 

DISP 

“LCAS" 

1 > 

I UX2=--=DS, 1, ' 

P2“. 1 

< 1 

GTAT 

DISP 

“LCAS“ 

1 ) 

! UY2-^0G, 2i * 

P2‘S 1 

( 1 

STAT 

DISP 

•‘LCAS“ 

1 > 

MJZ2--DSi 3, ‘ 

P2‘S i 

( 1 

GTAT 

DISP 

“LCAS“ 

1 ) 

! RX2^^DS/ 4/ • 

P2“/ 1 

( 1 

STAT 

DISP 

“LCAS” 

1 ) 

!RY2^DS> 5/ * 

P2'*/ 1 

( 1 

GTAT 

DISP 

“LCAS" 

1 ) 

!RZ2=DS/ 6/ * 

P2"/ 1 

( 1 

STAT 

DISP 

“LCAS*’ 

1 ) 

! P2=FREE( ) 

4: 







T' 

TABLE(NI^3/ 

N J"-^^ > : 

17 

EL21 

DISP 

•• ICOU“ 

0 

J=l: “UXl" 

“UYl “ 

II 

UZl" 




J^2; “RXl" 

‘*RY1*‘ 

** 

RZt “ 




J=3: “UX2“ 

••l.JY2‘* 

II 

UZ2“ 




J^-4; •*RX2“ 

“RY2‘' 

II 

RZ2“ 





DEFINE 811^13 DIRC E21 O O 

‘Pl = IC0U-l-^9 

!P1 

TABLE(NI=3, NJ^3) : 13 EL21 ROTA “ICOU" O 

TRANSFER (BGURC£=B 11/ ILIN^9/ SBA3E=*'P1" ) 
!Pl==FRr.E( ) 

DEFINE EL21 RDTA “ICOU“ 0 


DEFINE 812^^17 


EL21 DISP 

**ICQU' 

0 


13 DI21 

ELRF 

It 

ICQU" 

0=^RPROD(811/ 812) 


‘UXl^DS/ 

1/1/ 

1 

(13 

DI21 

ELRF 

"ICOU" 

0> 

lUYl^^'DS. 

2 / 1 / 

1 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

f UZl^DS/ 

3/ 1, 

J 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

!RXl=-DSi 

1/ 2/ 

i 

( 13 

DI21 

ELRF 

"ICQU" 

0) 

! RYl^DSi 

2, 2i 

1 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

IRZl-DS. 

3/ 2/ 

i 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

!UX2-0S, 

1/ 3/ 

1 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

! UY2-DS, 

2/ 3/ 

1 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

»UZ2^DS. 

3/ 3. 

1 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

!RX2^DS, 

1. 4/ 

1 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

•RY2==DS/ 

2/ 4/ 

1 

( 13 

DI21 

ELRF 

"ICOU" 

0) 

!RZ2-^DS/ 

3, n, 

1 

( 13 

DI21 

ELRF 

" ICOU" 

0) 


'L=I)S. “ ICaiJ“/ 1/ 1 (13 LEND E21 O 0) 

! ICOU 

IXGPl-i. O'-O. t?7/35^LaO. 5 
! XGP2-1. 0 *0. 57735^L^<0. 5 
! Al= 1. 0/L^UZi 
‘ A2=l 0/L*tTZ2 
‘UZGl^Al *A2 
»UZG2=Al'*A2 
! Al=l 0/L 
! A2-=Al >A1 
I A3-‘A2^A1 
!El=-6. 0^^A2 
f E2-=12. 0^XGPl ttA3 
'F3^-4. 0^<A1 
!E4=6. 0^X0P1«A2 
!E5^E3-J*0. *7 
f WXG1^E1 h“E2^UY1 
• WXG1^-E3- E^^R:<1+WX(:)1 
! N X G I E 1 - E2 ^ U Y2h W X Q 1 
INXGl— •E5~E4^t-RX2+WXGl 
! WYG1=-£1 + £2^UX1 
! WYGl=E3 *E4itRYl-*-WYGl 
! WYG 1 =-E 1 -E2ttUX2+WYG 1 
! WY(n^ES+E4>RY2-HWY01 


! t:2=l2. 0*XGP2*A3 

• 0^XGP2^A2 

• WXG2^^lTH b*:;25»UYl 

! WXG2 — E3-E4 >^K X 1 +WX02 
! UXG2^-£j -E2 ifrUY2 1 WXG2 
' UIXG2^-£5-E4^RX2+WXG2 
f WYG2=£l>£2i«UXl 
! WYG2^£3 +^£4 ifR Y 1 +WYG2 
! WYG2^'- E 1 “E2-rt-UX2+WYG2 

• WYG2^^£5+E4'JtHY2+WVG2 
! XGPl^FREEt ) 

! XGP2 =FRE£( ) 

! A1 =FRt-E( ) 

’ A2=FRE£( ) 

! A3-FREE ( ) 

!L=^FREE( ) 
f £1=FR£E( > 

!£2=KREE( ) 

!E3^FREE< ) 

!E4-FRt;t ( ) 

!ES-FREE( ) 

TABLE ( NI ==^=2, NJ=--3) : 13 ELST E21 

J=l; *'UZG1'* *‘UZG2*‘ 

J-2: ••WXGl" '‘14X02“ 

j.^3- m^^yGI" “WYG2'' 

• TEriP-icrju-Dt:2i 
-^jr4Z(TEMP/ 100) 

!TEMP=FREE( ) 

I ICGU-P‘REE( ) 

!UZG1=--FREE< ) 

‘UZG2---FREE< ) 

*WXG1--FRCE( ) 
iWXG2=FREE< ) 

! WYG1=FREE< ) 

•WYG2=-FREE( ) 
f UX1=FR£E( ) 

‘UY1=FRFE( > 

•U21=:FRFE( ) 
fRXl-=FR£E( ) 

!RY1-^FREE( ) 

!RZ1=FREE( > 

!UX2=^FREE( ) 

*UY2-=FR£F ( ) 

! U22=FREE< ) 

»RX2-^FR£E< ) 

!RY2=FREE( ) 
f RZ2^FREE( ) 

«XQT DCU 
ERASE 17 
^RETURN 

EN0PE21 

^ 

«M29 SEGA E21 0 0) 

% 

! TEMP^O 

^XGT AUS 
ICOU^O 


©RIGIMAL PAGE IS 
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ICOU“ 0 


ENDSC21 



^l.ADEL iOO 
‘ ICGU^ICDU-» 1 

! rEf1P = OS (1 tD21 REL 0 O) 

! ICDV 
! TEMP 

! TEMP ==TEhP- ICDV 
JfrJN2(Te:PIP. 200) 


!P1=DS, 13, “ ICOU", 1 

(13 CONN E21 

0 0) 


IP2=nS, 14, “ICOU", 1 

(13 CONN E21 

0 0) 


lUX 1. "PI" 

. 1 

( 1 

ST AT 

DISP 

"ADLC" 

1 ) 

‘UYl-^DS. 2/ "PI" 

. 1 

( 1 

GTmT 

DISP 

"ADLC" 

1 ) 

MJZi^OS, 3, -PI'* 

i 1 

< 1 

STAT 

DISP 

"ADLC" 

1 ) 

!RX1^D3. 4, "PI" 

. 1 

( 1 

STAT 

DISP 

"ADLC " 

1 ) 

'NYl^DS. 5; "PI" 

/ 1 

( 1 

STAT 

DISP 

"ADLC" 

1 ) 

‘RZl==D3i 6. "Pi" 

. 1 

( 1 

STAT 

DISP 

"ADLC" 

1 ) 

'UX2=0G. ii "P2" 

. 1 

( 1 

STAT 

DISP 

"ADLC" 

1 ) 

! UY2==03i 2/ "P2" 

, 1 

( 1 

STAT 

DISP 

"ADLC" 

1 ) 

?UZ2-DS. 3, "P2“ 

, 1 

( 1 

STAT 

DISP 

"ADLC" 

1 ) 

!RX2^=DG. 4. "P2" 

, 1 

( 1 

STAT 

DISP 

"ADl.C" 

1 ) 

!RY2^DS. 5. "P2" 

, 1 

( 1 

STAT 

DISP 

"ADLC" 

1 ) 

f RZ2^0Si 6s "P2" 
!P1^FHEE( ) 
!P2=FREE( > 

* 

table (Ni ^3, NJ-: 

, 1 

( 1 

STAT 

DISP 

"ADLC" 

1 > 

4) : 

17 

EL21 

DISP 

"ICOU" 

"ADLC 

J--=l. "UXl" "UYi" 

II 

UZl" 




J^2: "RXl" "RYl" 

II 

RZ1“ 




"UX2" "UY2" 

II 

UZ2" 




J^4; "RX2" "RY2" 

•' 

RZ2" 




DEFINE B 11-13 

EL21 

ROTA " 

ICOU" 

0 


DEFINE B 12=^1 7 

EL21 

DISP " 

ICOU" 

"ADLC " 



17 DI21 ELRP~ "1000“ "ADLC'*-RPROD( Bll, B12) 


UXl =:DSi 

li li 

1 

( 17 

DI21 

ELRF 

"ICOU" 

"ADLC" ) 

UYl-=DSi 

2 / 1 / 

1 

( 17 

DI21 

ELRF 

"ICOU" 

"ADLC" ) 

UZl -DSi 

3i 1, 

1 

( 17 

DI21 

ELRF 

"ICOU" 

"ADLC" ) 

RXl^DSi 

li 2i 

1 

( 17 

DI21 

ELRF 

"ICOU" 

“ADLC" > 

RYl^DSi 

2i 2i 

1 

( 17 

D121 

ELRF 

"ICOU" 

"ADLC " ) 

RZ1-=DS, 

3i 2i 

1 

< 17 

DI21 

ELRF 

" ICOU" 

"ADLC" ) 

UX2^DSi 

li 3i 

1 

( 17 

DI21 

ELRF 

" ICOU" 

"ADLC" ) 

UY2^DSi 

2i 3i 

1 

( 17 

DI21 

ELRF 

"ICOU" 

"ADLC" ) 

UZ2=DSi 

3i 3i 

1 

( 17 

DI21 

ELRF 

"ICOU" 

"ADLC" > 

RX2-DSi 

li 4i 

1 

< 17 

DI21 

ELRF 

"ICOU" 

"ADLC" ) 

RY2-^DSi 

2i 4i 

1 

( 17 

DI21 

ELRF 

"ICOU" 

"ADLC" ) 

RZ2“DSi 

3> 4i 

1 

( 17 

DI21 

ELRF 

"ICOU" 

"ADLC " ) 


! ICGU 

!L-^D5, "ICOU". 1 , 1 (13 LENG E21 O O) 

!XGP1 = 1, 0-0. 57/35-i^L-i*0. 5 

fXGP2--=l. 0 + 0. 57735+^L-MO. 5 

* Al =-1. 0/L+4-UZ1 

■ A2=l. 0/L*UZ2 

!LUZ1=A1 ^A2 

!LUZ2=AH A2 

! Al^l. 0/L 

! A3^A1*A1 

! A3-A2>A1 

I El ^-6. 0i+A2 

! E2=12. 0++X0Pl+^A3 

lE3=-4 0++Ai 

!E4^6. 0-f+XGPl^A2 



92 


E5--fc:3*0. t'l 
WXGl-El+E2■i^UY1. 
WXGl-=-E3-E4it-RXl+WXGl 
WXG1--EX-E2’HJY2+WXG1 
LWXl --ES-E4*RX2+WXG1 
WYGl^El+ETliiiUXI. 
WYGl-E3+E4-»RYi+WYGl 
WYG1=-E1-E2 *UX?h WYO l 
LWYl -E5+E4«RY2 -kWYG 1 
E2=12. 0^»XGP2«A3 
E4-6. 0itXGF2-«A2 
WXG2-~E1+E2>HJV ), 
WXG2=-E3-E4»RXi+WXG2 
WXG2--E1-E2*UY2 h WX 02 
LWX2^--E5-E4*RX2-*WXG2 
WYG2 --Ei+E;2«UX J. 

W YG2-^E3-» E4»R Y X +W Y02 
WYG2— E t-E2»UX2H WY02 
L WY2^E5+E4*R Y2+WYG2 
XGPX=FRe.E( ) 
XGP2-FREE( ) 

A3 -FREE ( ) 

£3=FRF£( ) 

E4=FRC:E( > 
t£3-=F-Rb:E( ) 


UZG1 = DS. 1. 1. 1 

( 13 

ELST 

E21 

"ICOU" 

0) 

UZG2=DS/ 2, i, 1 

( 13 

ELST 

E21 

"ICOU" 

0) 

WXGl^DS, 1, 2. 1 

( 13 

ELST 

E21 

"ICOU" 

0) 

WXG2=0S; 2/ 2> 1 

( 13 

ELST 

E21 

"ICOU" 

0) 

WY01=DS; 1, 3, 1 

( 13 

ELST 

E21 

"ICOU" 

0) 

WYG2=^0S, 2; 3, 1 

( 13 

ELST 

E21 

"ICOU" 

0) 


it 

* COMPUTE DESIGNS SENSITIVITY 


!E=DS, X, “ICOU", X <X3 MATP E2X O 0) 

• D -DS, 2, ■* ICDU", X ( X3 SPLO E2X O 0) 

in-DS, X, "ICOU", X (X3 SPLO E2X 0 0) 

! D--U-M-2. 0 
0 

!WEIK=DS X,"ICDV",X (X DESV VALU 0 0) 
!NEIH=DS 2, "ICDV", X (X DESV VALU 0 0) 
! D£LB=X. 0«l^)EID^< B 
!DELH=X. 0^<WEIH*H 
! ICDV 

! AX =--UZGX-;iLUZX 
' A2-UZG2»LUZ2 
! E X --AX +A2*E*H»DELB*L/2 
! £2 =A X +A2*F>B*DELH*L/2 
!PST=PSI-EX-E2 
: £X •=AX+A2*E*H»WEIB*L/2 
!E2=AX+A2*E*B-»WEIII»L/2 
! BD£L=-EX-;i X. 0 
!HDEL=-E2^H. 0 
! AX -WXGX-aLWXX 
! A2-WXG^^^LWX2 

! EX-=AX-* A2*£»B^B*B*DELH«L/24. O 
!E2-AX+A2*3. 0*fe>H»B*B*DELB»L/24. O 
•PSI^PSI-EX-E2 

!EX=AX+A2*E»B^tB<H)*WEIH-»L/24. 0 
!E2-AX i-A2^»3. O^E<*H-i< B^»B*WE I B*L/24. O 
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! nDF-L-’-^^DOEL' E2 
!HDEL~HI)EL“-E1 
! Al^WYGl^LWYi 
! A2=WVG2iU,WY2 

! El=Al+A2*t:*H^HitH^l)ELB-*i’L/24. O 
?E2-A1*^A2-«3. 0^^E*B^H*H*DELH*L/24. O 
!PBI'--Pai-E1*E2 

! El -A l+A2i^E-^Hi^H*>H-« WEIB^^L/24. 0 
! E2^Al-»-A2->^3. Oi^E^B^H^H*WEIH*M'L/24, 0 
! DDEL=BUEL-E1 
•HDEL=HD£L-E2 
IPS I 
% 

I rEMP=CLX“CCOM-CDIS“CS21-CS4l-l 
! TEHP=^XLC-CCOM- CDIS“l 
itJL.Z(TEhP. 1^5) 

I TEhP=TEHPi l-CS2i 
#JGZ (TEMP, 155) 

I TEMP=-CLC-CCOM-Cf) I S 

!JENT=DS "TEMP", 1 ,1 (1 ST21 LIST 1 1) 

I JEN r 

! TEMP‘D I COU-JENT 
I JEN'l =f-.REF( ) 

*JNZ(TEMP, 155) 

ISlB-1. 0 
!SIH=1. 0 
! TEMP=PniN-2 
itJNZ (TEMP, 444) 

' SIB--1. 0 
I S I H= i. . 0 
*JUMP(448) 
ii-LABEL 444 
1 TEMP=PDIN-3 
^^JNZ<TEMP, 446) 

• SIB^^“1. 0 

0 

*JUMP(44a) 
it-LABEL 446 
I TEMP-PaiN-4 
^tJNZH EMP, 448) 

'SIB-1. 0 
> 311-1-=^- 1. O 
i^LABEL 448 
‘SIB 
! SIH 

it- 

‘ INTl=WXGl + WXG2^E-«-0. 25 h^DEUB-»'SII3 

• INT2==WY01+WYC2-JtE'^0. 25^DELH*»SIH 
IPSl 

I PSI=PSI-INT1-INT2 
IPSI 

‘ INT1=WXG1+WXG2^P:-»-0. 25*WEIB*SIB 
I INT2=WYGl+WYQ2^E^tO. 25*WEIH*SIH 
!HDEL=HDEL-INT2 

• BDEL==BDP.L“INn 
! INTl^FREEl ) 

! INT2^FRE£( ) 

!SIB=FREE( ) 

!SIH^FREE( > 

$ 


♦LABEL 15S 
! TEMP--FREE< > 

! DELB-FREE ( > 

! DELII--FREE( ) 

S' 

I ICOU 
IPS I 
I tlDEL 
! UDEL 

! A1=FREE( ) 

I A2==FREL ( ) 

!E1'='^FRI-:E( ) 

IE2-FREEC ) 

•L -:FREE( ) 

I SIGl -FREE( > 

ISIG2=FREE< ) 
i 

■iiLABEL 200 
I TEMP=ICQU-DE2i 
♦.JNZ^TF.MP. 100> 

' ICOU--FREE( ) 

I 1£MP~FREE< ) 

>VJZ01=-FREE( ) 

•UZ02-PUEE( ) 

IWXG1=-FREE( ) 

IWXG2=FRE£( ) 

•WYG1-^FREE( ) 

! WYG2-FREE< ) 

‘LUZ1-^-FRE£( ) 

!LUZ2=FREE( ) 

!LWX1=--FREE< > 

•LWX2-=^FREE( ) 

!LWY1-FRL£( ) 

•LWY2=FR£E< ) 

IE --FREE ( ) 

> i3=FREE( ) 

IH^FREE< ) 

!UX1-FREE( ) 

! UY 1 -FREE ( ) 
lUZl =FREE( ) 

!RX1^FREE( > 

1RY1-FREE( ) 

•RZ1^^FREE( ) 

!UX2=FREE( ) 

!UY2=FREE( ) 

IUZ2-FRE£< ) 

IRX2-FR£E( ) 

!RY2-FR£E( ) 

!RZ2=FREE( ) 

♦RETURN 

♦ 

S. 

•4- 

<H29 AD JO E2i O O) 

* 

♦XOT AUS 

t COMPUTE AVERAGE ELEMENT STRESS 
* 
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UZG1=DS, 1, 1 (13 ELST E21 "JENT‘* 0) 

UZG2=DS. 2. 1. i (13 ELST E21 “JENT" 0) 

WXGi=DS, li 2, 1 (13 ELST E21 “JENf 0) 

WXG2^DS, 2. 2, 1 (13 ELST E21 “JENT*' 0) 

WYG1=DS, 1. 3i 1 <13 ELST E21 '‘JENT“ 0) 

WYG2=DS, 2. 3, i (13 ELST E21 “JENT“ O) 

L=DS, “vJENT“/ LI (13 LENG E21 0 0) 

£ = DS, 1. “JENT". 1 (13 MATP E21 0 0) 

I3==DS, 2. “JENT", 1 (13 SPLO E21 0 0) 

H=D3, L “JENT". 1 (13 SPLQ E21 O 0) 

O 

0 

H 

U 

SIG1=^E*UZGI 

SIG1=-E*0. ^j-»B^i'WXGl+SIGl 
SIG1=-E*0. 5*H»WYG1+SIGI 
SIG2^E*UZG2 

SIG2=-E«-0. 5*B*WXG2+SIG2 
3IG2^'-E^0. 5i^H«WYG2+SIC2 
VMSl==SlGl+SIG2-itO. 5 
SIG1-=EMJZG1 

5IG1^+E>0. ^7>D^^WXG1-»SIG1 
SIG1=- E^O. t!)^Hit-WYGl-»-SIGl 
SlG2-^Ei^UZG2 

SI(52=+£*0. 5^^D*WXG2+SIC2 
SIG2^-E^tO. 5*H*WYG2-*“SIG2 
VMS2^SIG1+SIG2«0. 5 
SIG1=E*U2G1 

SIG1=+E*0. S-J^D^WXGT+SIGl 
SXG1=+E*0. 5-»H*WYGl+SIGl 
SIG2“--E^UZG2 

SIG2^+E*0. 5*B*WXG2+SIG2 
SIG2.-=+E*0. 5^^H^<*WYG2+SIG2 
VMS3=5IG1+SIG2*0. 5 
SIG1=E*UZG1 

SIGI=-E*0. 5-«-D^^WXGl+SIGl 
SIG1=+E^^0. fj^tHJ^WYGi+SlGl 
SIG2=E^UZG2 

SIG2^~E*0. 5»I5^K4XG2+SIG2 
SIG2^-H£«0. 5^iH^t-WYG2+SrG2 
VMS4^SIGH-SIG2-i^O. 5 
UZGl=-f-REF( ) 

UZG2=^FREE( ) 

WXQ1--=FREF< ) 

UXG2^FREE( ) 

WYGl^FREFC ) 

WYG2=‘-FREF( ) 

VHS 1 
VMS2 
V11S3 
VHG4 

YMCl---ABS(VhSn 

VMC2=ACS(VMS2) 

VMC3^ADS( VhS3) 

VMC4=ACS(VHS4) 

VMST^VMSl 
VnCS^VMCl 
P0IN--=1 
GIB-1 0 


0 

!C0MP=^VriCS-VMC2 
J^v^QZCCOflPi 444) 

• VM3T-VHS2 
!PniN=2 
!GIB=- 1. 0 
! 3 1 H= 1 0 
i^LABEL 444 
I C0MP^^VriCS“VMC3 
JGZiCUMPi 445) 
f VMST=VI'i33 
!PCHN^3 

0 

O 

*L.ABEL 445 

• C0I1P=:=VMCB VMC4 
^JGZ (C(3hP, 446) 

! VhST^VhS4 
!POIN=-4 

O 

*SIH='-1 0 
^l.ABEL 446 
! VPiST 
IPO IN 

ICOMP=P'REE( ) 

! VNS1*=PREE( ) 

I VNS2--=FR£E( ) 

' VHS3=P'REE ( ) 

I VhS4=FPx£E( ) 

! VhCl^F«FE-:( ) 

I VHC2=FR£E( ) 

! VMC3=FR£F.( ) 

I VhC4=FRE£( ) 

I VMCS=FR£E( ) 

!SIG1=^F‘REE< ) 

!SIG2==KREE( ) 

IXGP1 = 1. 0-0. 57735^L«0. 5 
!XGP2=1. 0+0. 57735«L^^0. 5 
I AD3=*-1. 0*E/L 
I AD9=1. O^E/L 
I Al^l. O/L 
I A2=A1^AI 
I A3=A2*Ai 
!El*-^~6. 0i^A2 
!E2=12. 0^XGP1*A3 
IE3=-4. O^tAl 
•E4---6 0^XGPH*A2 
IE5=E3-^0. 5 
lf'6==12. 0^XGP2»A3 
‘£7=6. 0^XGP2^A2 
•AD=2 0«El+E2+E6/2. 0 
! AD2=-AD*E4^SIB^B/2. 0 

• ADl=-Ar/^E^SIH^H/2. 0 
•AD=-2. 0i^El-E2«E6/2. 0 
‘ AD8=-AD+E^5IQ+B/2- 0 

• AD7=-AD^tE*SIH*H/2. 0 
‘AD=2. 0-J*E3+E4+E7/2. 0 
I AD4=AD*E*SIB^B/2. 0 

I AD5=-Ar)*E^*-SIH^H/2. O 
•AD=2. O^E5-»E4fE7/2. O 


ORTGINAE PAGE 

C •{' p-'>0 ; ' O" ’ ' ’ : ■ 



! AD10=AD^fc:^SIB^B/2. 0 
!A011=-AD*E^SIH^H/2. 0 
! AD=FRf:E( ) 

! XGPl=Fnt.E< ) 

! XGP2=FF^EE( ) 

!F-=FREE( ) 

IB-FREE ( ) 

!M=^FREE( ) 

!L:-FREE< ) 

! A1-=FREE( ) 

I A2=FREE( ) 

! A3-FREE( ) 
lEl-FREEC ) 

! E2-FREE( ) 

I E3-l- REE( ) 

!£4=FREE( ) 

IE5=FRF.£( ) 

!E6-FREE( > 

!E7=^FREE< ) 

!SIB-FREF( ) 

!SIH-FREE< ) 

TABUE(NI=3, NJ:^4) : 17 ROTA AD21 “AOL.C** 0 


J=1 ; 

’•ADI" 

‘‘AD2" 

"AD3" 

J=2; 

"AD4" 

"AOS" 

0. 0 

J=3: 

"AD7" 

"ADC" 

"AD9" 

xJ=4: 

"ADIO" 

"ADll" 

0. 0 


DEFINE Q 11-13 EL21 ROTA *\iENT*' 0 
17 EL21 RGTT "JENT" O-RTRAN(Bll) 
DEFINE B 12-1 7 EL21 RQTT “JENT“ 0 
DEFINE B 11=17 ROTA A021 “ADLC" 0 
17 CLAD E21 “ADLC" O-RPROD ( B 1 2. B 1 i > 


ADl=DirS, 

1, 1, 1 

( 17 

GLAD 

E21 

"ADLC" 

0) 

AD2=DB, 

2, 1, 1 

( 17 

GLAD 

E21 

"ADLC" 

0) 

AD3-DS; 

3, 1, 1 

(17 

GLAD 

£21 

"ADLC" 

0) 

AD4-DS, 

1, 2, 1 

( 17 

GLAD 

£21 

"ADLC" 

0) 

AD5-D3, 

2, 2, 1 

( 17 

GLAD 

E24 

"ADLC" 

0) 

AD6-DS, 

3, 2, 1 

(17 

Gl.AD 

£21 

"ADLC" 

0) 

AD7-DS, 

i, 3, 1 

(17 

GLAD 

£21 

"ADLC" 

0) 

AD8-DS, 

2, 3, 1 

( 17 

GLAD 

E21 

"ADLC" 

0) 

AD9=DS, 

3, 3, 1 

( 17 

GLAD 

E21 

"ADLC" 

0) 

AD10=1)S, 

1, 4, 1 

( 17 

' GLAD 

^ E21 

"ADLC 

* 0) 

ADI 1=DS, 

2, 4, 1 

< 17 

' GLAD 

> £21 

"ADLC 

• 0) 

AD12-DS, 

3, 4, 1 

( 17 

GLAD 

^ E21 

"ADLC 

• 0) 


IPl-DS, 13, "JENT", 1 <13 CONN E21 0 0) 
IP2-DS, 14, "JENT". 1 (13 CONN E21 0 0) 


SYSVEC: 

APPL FORC 

1 = 1 

j n 

; "ADI" 

J--P2" 

: "AD7" 

1=2 

J="P1" 

: "AD2" 


: "AD8" 

1=3 

J=“P1" 

"A03" 

J="P2" 

: "AD9" 

I-=4 

J=:"P1" 

: "AD4" 

J="P2" 

: "ADIO" 

1=1') 

j^npiM 

"AD!?" 


"ADll 
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J^“P^ “ 

1-6 

J=^“P5“ 
!P1=FRC 
!H2^FRt:E:( 

! ADI REE. 

! AD2-KREE 
• AD3-=FRE£ 

? AD4 =^FHEE 
! ADS^FREE 
! ADo^FREt". 

» AD7=FREE 
! AD8-=FREE 
! AD9^FREE 
? AD10^=^FRE:i 


: "AD6" 

; "AD12“ 

( ) 

) 

( ) 

) 

) 

) 

) 

) 

) 

) 

) 

( ) 


! ADI 1^FREE( ) 
? AD12--FREE( ) 
*i^XQT Ui 


^SHQW 
^XQT DCU 
ERASE 17 
^RETURN 
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^ E^4DAD21 

^ ( 29 PREP E41 O 0 ) ENDPE41 

‘F 

% THIS SUUROUTINE IS FOR ELEMENT TYP E41 
% 

* THE PRINCIPAL LIBRARY IS 14 

€ TEMPORY LIBRARY IS 15 CAN BE DELETED AT THE END OF THIS ROUTINE 
i^XQT Ell 
EXTRACT: 


EXTRACT: 
CREATE 1 

EXTRACT: 
CREATE 1 
% 

EXTRACT: 
CREATE 1 

a> 

EXTRACT; 
CREATE 1 
'4. 

EXTRACT: 


S0URCE-E41: 
AREA E41 0 

0 

CONTENT 

SPEC: 

GEOM 

S0URCE^E41: 
ELRF E41 0 

O 

CONTENT 

SPEC: 

oeoM 

S0URCE^F.41: 
ELNO £41 0 

0 

CONTENT 

SPEC: 

GEOM 

S0URCE-E41: 
MATP £41 0 

0 

CONTENT 

SPEC: 

MATE 

S0URCE-=E41: 
CONN E41 0 

0 

CONTENT 

SPEC: 

INTE 

SQURCE-E41: 
HIT E41 0 0 

CONTENT 

SPEC: 

S 1, ; 


1 ^ AREA OF ELEMENTS 


5. 12 % LOCAL ELEM. REF. FRAME 


22, 57 % ELEMENT-NODAL RELATION 


$ CONNECTIVITY 


HINV T 


^XGT DCU 

TGCC 14 CONN E41 0 O 
^XGT Ul 

'ENTl^O. 6220084 
•ENT2^-1. 0/6. O 
!ENT3=0. 0446582 
^^TI(14 SNAP FUNC 1 1) 


TYPE-0 



"ENTl " “EN'r2" 

"ENT3“ 

"ENT2" 

99 


"ENT2" "FN'l 1" 

"ENT2" 

"ENT3" 



"ENT3" "EN 1 2" 

"ENTl " 

"ENT2" 



"ENT2" "ENT3" 
!FNT1=FRFE( ) 
IENr2=KREE( ) 
!£NT3=FRFE( ) 

$ 

KXQT AUS 

"FNT2" 

"ENTl" 



I ICOU=l 
f TEMP=0 





^^LADEL 200 





! TEMP=FREE( > 
!B1 = IC0U-U<8 






DEFINE Al--^14 ELRF E41 0 0 
TABLE(NI=^2. NJ-=4> ; 15 ELEM REFT “ICOU” 99 
TRANSFER (SOURCE^Al, ILIM^S* SBASE=“Q 1 ‘* ) 
!B1=FREE( ) 

DEFINE BQi-15 El.EM REFT “ ICOU*' 99 

14 ELEM REFE “ICOU" 0-^RTRAN ( BQl ) 

* COMPUTE DETERMINATES FOR INTEGRATION 
! 1L=1 

!GPl=--0. 57/35 
*LABEL 55 
!PXI^-.- 1. OiiGPl 
!HFr-=~l. O'iiOPl 
! 1 FMP^IL-S 
^JNZ ( TEMP, 35) 

IPXI-+1 O'itCPl 
f PET'^-1. OaGPl 
^LABEL 35 
' TEMP^^^IL-*3 
JtJNZ(TEr-iP, 37) 

•PXI-=-i 1. 0»GPi 
‘PET^^-^ l. OGPl 
:n..AB£L 37 
! TEMP= IL- 4 
i^JNZ(TEMP, 39) 
f PXI=“i. 0>GP1 
•PET-^ t l. O^OPi 
^vLABEL 39 

‘DUM1==--1. 0^-PET>0 25 
!DUM2=^~1. O+PXr^O. 25 
!DUM3--“-» i. 0-PE1>C). 25 
IDUM4--1. O-FLXI-i^O. 25 
!L)UM5^+1. 0+PEI*0. 25 
! DUM6'=+1. O+PXI itO. 25 
! DUM7^~1 . O'-PET^O. 25 
!DUM8=+l. O-PXn^O, 25 

TABLE(NI==2, NJ^-4) : 15 HELP E41 “IL" 99 

J=1 : “DUMl" “DUM2" 

J-2 ; “DUM3" “DUM4" 

J=3 ; “DUM5" “DUI16" 

J-4 : "DUM7" “DUMB" 

DEFINE El =15 HELP £41 “ IL" 99 
DEFINE EB2-14 ELEM REFE “ICOU" 0 

15 JAC E41 "ICOU" "IL“=RPROD (E1,BB2) 
!DUM1=DS 1,1,1 (15 JAC E41 "ICOU" "IL") 
*DUM2=DS 1,2,1 (15 JAC E41 "ICOU" "IL") 
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»DUM3^DS 2,1/1 <15 JAC E41 •’ICOU’* '‘IL“> 
!DUr14-=DG 2/2,1 (15 JAC £41 " ICOU" “IL“) 
i DUH5-^DUIil*frI)Ur'14 
! DUri6^DUM2^tDUM3 
! DUH8=nUl15-DUM6 
! DUMB 

! TEMP^ IL-1 
*JNZ(TEHP/ 45) 

! OETl-^DUhB 
^JUflP 53 
LABEL 45 
! TEriP==IL “2 

( TEhP/ 47) 

» DE T 2^DUM8 
53 

it^LABEL, 47 
! TEMP=IL-3 
ifv)NZ(TEf4P/ 49) 

! DET3=^DUM8 
^JUMP 53 
< LABEL 49 
! TEMP^II .-4 
^JNZ(TEMP, 53) 

! DET4-=DUM8 
^ LABEL 53 
* IL^IL-t-1 
f TEMP-*-- IL-5 
frJNZ( TEMP/ 55) 

TABl-E (141 = 4/ NJ>=1 ): 14 DETE E41 "ICOU" 0 
J=l, "DETl” "DET2" "DET3" "DET4“ 

! DETl=FRt:E( ) 

I DET2=EREE( ) 

! DET3=FR£E( ) 

!DFr4=FREE< ) 

!DUHl=rRFE( ) 

'DUM2=FREE( ) 

! DU13=FREE( ) 

I DUM4=FREE( ) 

!DUM5=FREE( ) 

!DUM6=FRt;E( ) 
f DUM7 =FREE( ) 

!DUM8=FREE( ) 

! IL=FREE( ) 
f PXI=FREF( ) 

!PET=FREE( ) 

f gpi =fre:e( ) 


DEFINE BB1=14 


SNAP 

FUNC 

1 1 



DEFINE BD 

2=14 


EI.EM 

REFE 

“ icau 

*' 0 


15 GAUS POIN 

cC 

II 

ICOU 

* 0= RPROD 

(BBl/ BB2) 

! XGPl^DS 

1/ 1/ 

1 

( 15 

GAUS 

POIN 

“ICOU" 

0) 

I YGPl^DS 

1/ 2/ 

1 

( 15 

GAUS 

POIN 

»ICOU“ 

0) 

1 XGP2^DS 

2, 1, 

1 

(15 

GAUS 

POIN 

**ICOU" 

0) 

! S'GP2=DS 

2/ 2/ 

1 

(15 

GAUS 

POIN 

** ICOU" 

0) 

! XGP3=D3 

3/ 1, 

1 

(15 

GAUS 

POIN 

"ICOU" 

0) 

! YGP3-US 

3/ 2/ 

1 

(15 

GAUS 

POIN 

" ICOU" 

0) 

• XGP4=DS 

4/ 1, 

1 

(15 

GAUS 

POIN 

"ICOU" 

0) 

' YGP4=DS 

4, 2/ 

1 

( 15 

GAUS 

POIN 

"ICOU" 

0) 
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TABLE(NI=1 

2, NJ= 

5) : 14 

PMAT 

E41 "ICOU" 0 







J=1 

0. 0 

0. 0 

1. 0 

0. 0 

0. 0 1.0 

0. 0 

0. 0 

1. 0 

0. 0 

0. 0 

1. 0 

J=2 

0. 0 

1. 0 

0. 0 

0. 0 

1.0 0. 0 

0- 0 

1. 0 

0. 0 

0. 0 

1. 0 

0. 0 

J=-'=3 

1. 0 

0. 0 

0. 0 

1. 0 

0. 0, 0. 0 

1. 0 

0. 0 

0. 0 

1. 0 

0. 0 

0. 0 

J=4 

"YGPl 

“ 0. 0 

0 0 

•'YGP2 

" 0. 0 0. 0 

"YGP3** 

0. 0 

0. 0 

"YGP4" 

0. 0 

0. 0 

J=5 

0. 0 " 

XGPl" 

0. 0 

0. 0 

"XGP2" 0. 0 

0. 0 “ 

XGP3'* 

0. 0 

0. 0 '* 

XGP4“ 

0. 0 


! XGPl=FR6lr"( ) 

! YGPi^FRet: ( ) 

! XGP2=FREE( ) 

! YGP2^KRF.E( ) 

! XGP3^FR6:F< ) 

! YGP3=FREE( ) 

» XGP4^FREE( ) 

‘ YGP4^FREE( ) 

% 

!NU=OS 2, "ICaU".l (14 MATP E41 O 0) 
!r4UPl^NUH 1. O 
(NU=-1. O^tNU 

TABLE(^4I = 12. h4J^12) : 14 SIEP E41 “ICOU‘* O 


J- 

= 1 

1 

0 


"NU 

» 

0. 

0 

0 

. 0 0. 0 0. 

0 


0 

0 

0 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

J~- 

=2 

“NU* 


1. 0 

0. 

0 

0 

. 0 0. 0 0. 

0 


0 

0 

0 

0 

0. 

0 

0. 

o 

0. 

0 

0. 

0 

J= 


0 

0 

0. 0 ‘ 

'NUPl" 

0 

. 0 0. 0 0. 

0 


0 

0 

0 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 


^4 

0 

0 

0 

0 

0. 

0 

1. 

0 

"NU" 0. 

0 


0 

0 

0 

0 

0. 

0 

0, 

0 

0. 

0 

0. 

0 


= 5 

0 

0 

0 

0 

0. 

0 

"NU" 

1.0 0. 

0 


0 

0 

0 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

J 

=6: 

0 

0 

0 

o 

0. 

0 

0. 

0 

0. 0 "NUPl 

II 


o, 

. 0 

0. 

0 

0. 

0 

0 

0 

0. 

0 

0. 

0 

J 


0 

0 

0 

0 

0. 

0 

0. 

0 

0 0 0. 0 

1. 

0 

"NU 

II 


0. 

0 

0. 

0 

0. 

0 

0. 

0 

J 

-8. 

0 

0 

0 

0 

0. 

0 

0. 

o 

0. 0 0. 0 

"NU‘ 

1 

1. 

0 


0. 

o 

0. 

o 

o. 

0 

0. 

o 

J 

-9 

0 

0 

0 

. 0 

0. 

0 

0 

0 

0. 0 0. 0 

0 

0 

0 

0 

"NUPl ' 

1 

0. 

0 

o, 

0 

0. 

0 


J- 

10; 

0. 

0 

0. 

0 

0. 

0 

0 

0 

o. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

1 

. 0 

**NU" 0. 

0 

‘J- 

1 1 

0. 

0 

o. 

0 

0. 

0 

0. 

0 

0. 

0 

0, 

0 

0. 

0 

0- 

0 

0. 

0 

II 

NU* 

• 1,0 0. 

0 


12; 

o. 

0 

0. 

0 

0 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0 

. 0 

0,0 "NUPl 


!^4U==FREE( ) 

•NlJPF^FREEC ) 

!TEMP=FR£E( ) 

!B1 =ICGU“*i^5 
DEFINE A2^ESB E41 1 1 

TABLE(NI=tj, NJ^1):15 BBB E41 "ICOU** 0 
TRANSFERR(SaURCE--=A2i ILIM=5, SBASE^"B1" ) 

! Dl-FREEC ) 

% 

DEFINE BB1^14 PMAT £41 “ ICOU“ 0 
DEFINE BB2=15 BBB E41 "ICQU" O 
•SAGR^DS S. “ICOU'S 1 (14 CONN E41 O 0) 

!OOTH=DS 26» "SAGR"; 1 (1 SA BTAB 2 13) 

14 GPST E4i “ICOU** RPRQD ( "OQTH" BB1,BB2) 
!OOTH=FREE( ) 

!SAGR=FREE(> 

! ICOU:-ICQU'i 1 
! TEMP=DE41- ICOU+1 
>JCZ(TEhP. 200) 

! in=FREE( ) 

! ICOU-FREE( ) 

‘TEMP=FREE( ) 

-«XQr DON 
ERASE 15 
RE TURN 


^‘NDPE41 
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«-(29 SEC A E41 0 O) 
‘4> 

«XQT AUS 


ENDSC41 


ip 

!iCOU^l 
• Tt- l'IP-^0 
^<LaIjI"L 200 

‘ELGR^OS l,"lCQU'\i (1 ED41 REL 0 0 ) 
! TEMP-ELGR- ICDV 
*JNZ< TEMP, 195) 


ORIGINAL PAGE IS 
OF POOR QUAUTY 


IBl^^ICOU- 1-afi 

DEFIME A2=l ESB E41 "ADLC“ 1 
TADI.E(NI:-^5; 14J=1):15 BBB E41 “ICOU” “CLC‘ 

I RANSFERR (S0URCE^A2, ILIM=5, SBASE=‘*Dl '• ) 

‘ 01 ^-FRL“E( ) 
it 

nEFlNE BBl:-^t4 PMAT E41 “ ICOU“ 0 
DEFINE 032^15 BBB E41 "ICOU" "CLC" 

‘SACK''=DS S, '‘ICaU“/ 1 (14 CONN E4 1 0 0) 
lOOTH^OS 2&, *'SAGR“, 1 (1 SA BTAB 2 13) 

15 GPST E41 *'ICrjLr» “CLC“ = RPROD < ''DOTH*' BB1/BB2) 
inorH-FREEC ) 
iSAOR-FREE( ) 
it 

DEFINE Bil-14 SIEP E41 “ICOU'* 0 
DEFINE B12==15 GPST E4 1 "ICOU" "CLC" 

14 GFEP E41 "ICOU" "CLC"= RPROD <B11,B12) 

S 


!SIl -03 

1,1,1 

< 14 

GPST 

E41 


ICOU" 

0) 


!BI2 =D3 

2, 1, 1 

(14 

GPST 

E41 


ICOU" 

0) 


•GI3 --DS 

3, 1, 1 

(14 

GPST 

E4l 


ICOU " 

0) 


!S14 =DS 

4, 1, 1 

(14 

GPST 

E41 


ICOU" 

0) 


•SID 

5, 1, 1 

<14 

GPST 

E41 


ICOU" 

0) 


!SI6 -DS 

6, 1, 1 

<14 

GPST 

E41 


ICOU" 

0) 


rSIT ^DS 

7, i, 1 

< 14 

GPST 

E41 


ICOU" 

0) 


!3I8 ^DS 

8, li 1 

( 14 

GPST 

E41 


ICOU" 

0) 


!S1V =^03 

9, 1, 1 

(14 

GPST 

E41 


ICOU" 

0) 


•SI lO^DS 

10, 1, 1 

(14 

GPST 

E41 


"ICOU 

" 0) 


!SI1 1-DS 

11,1,1 

(14 

GPST 

E41 


"ICOU 

" 0) 


•SI 12-DS 

12, 1, 1 

(14 

GPST 

E41 


"ICOU 

•• 0) 


•EPl -D3 

1, 1, 1 

( 14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

•tP2 --^DS 

2, 1, 1 

(14 

GPEP 

E4l 


ICOU" 

"CLC" 

) 

•EPS ^DS 

3, 1, 1 

(14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

•EP4 --DS 

4, 1, 1 

( 14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

‘EPS ==03 

5, 1, 1 

(14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

!EP* -03 

6 , 1, 1 

(14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

•EP7 ^DS 

7, i, 1 

( 14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

•EPB ^DS 

8, 1, 1 

( 14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

!EP9 -OS 

9, 1, 1 

( 14 

GPEP 

E41 


ICOU" 

"CLC" 

) 

lEPlO-DS 

10, 1; 1 

(14 

• GPEP 

E41 


"ICOU 

M ..^LC" 


•£P1 l^DS 

11, 1, 1 

( 14 

GPEP 

E41 


"ICOU 

M 


?EP12=DS 

12, 1, 1 

<14 

GPEP 

£41 


"ICOU 

n hCLC" 



% 


:EP1 =SI l#EPl 
!EP2-3I2>F:P2 



!EP3=2. 0^SI3^^EP3 
!D£TE=DS 1.1.1 (14 DETE E4 1 
' EP3=EPl + EP2+EP3*D£rE 
! EP4-^3I4^frEP4 

!EP6:---2. 0^3I6*EP6 

!DETE=DS 2.1.1 (14 DETE E41 

!EP6--^EP4+£P5+EP6*DETE 

‘E-P7=SI7*EP7 

' EFB=SI8*EPn 

!EP9=2. 0-»«-SI9f^LP9 

!OETE=^OS 3.1.1 (14 DETE E4 1 

! EP9^tp7 + EPSi EP9^DETE 

! EP I O^SIlQJtEP 10 

‘EPl 1^3111-iiEPU 

' KP 12=2. 0*S I 12>EP 12 

!nETE=DS 4.1.1 (14 DETE E41 

! EP 1 2=EP 1 0-+ EP i i + EP 1 2 » DETE 


"ICOU" 0) 


•'ICQU’‘ 0> 


‘•ICOU" O) 


“ICOU" O) 


! f>SlS=EP3+EP6+EP9+-EP12 

!NU=DS l."ICOU'‘,l (14 MATP E41 0 0) 

‘3AGR=DS 8. "ICOU". 1 (14 CONN E41 O 0) 

!naTH=DS 26. "SAGR". 1 (1 SA BTAB 2 13) 

! PSIG^-PSIG/NU 
*PSl=P3I-+-PGIS 


* ICOU 
‘ PSI3 


!PSIG=FREE( ) 
!3AGR=FRt:E( ) 
! 00 VH=FREE( ) 


!NU=FREE( ) 

I 

'3I2l=SI2«SIl 
' SI54=SI5*GI4 
• 3I87=SI8'it'SI 7 
‘3I01=SI10-it3in 


!3T 1=SI1*GI1 
'GI2=SI2^^SI2 
!SI3=GI 3^^813 
! 314 = 31 4-t^GI 4 
!3I^=3I5-«GID 
‘ 316=316^316 
' GI7 = 3I7-«3I7 
'sr8=3ia*-3ra 
* 319=319*319 
‘SI 10=31 10*31 iO 
!SI n=3lll*SIli 
‘3112=3112*3112 
•I 

!SI3=GI3*3. 0+-SI i+SI2~SI21 

*Dfc'rE=DS 1.1.1 (14 DETE E41 "ICOU" 0) 

‘313=313**0. Si^DETE 

1316=316*3. 0-rSl4< BI5~SI54 

‘DFTE=DS 2.1.1 (14 DETE E41 "ICOU" O) 

! 316=316**0. 5*DETE 

1319=319*3. 0+'Sr/-*-SI8-*SI87 

ID£TE=DS 3.1.1 (14 DETE E41 "ICOU" O) 

‘319=319**0. 5*DETE 

! SI 12=3112*3 0+SI 10+SI 1 1-3 101 

‘DETE=DS 4,1,1 (14 DETE E41 "ICOU" O) 


104 


GI 12-2112^^0. 5*DETE 
Vh2T^2 1 3-^S 1 6+S 1 9+S 1 1 2 


‘AREA^-^OS t, '‘ICOU‘\ 
! VMS r-'=VhST/ AREA 
! ARr£A =^FREE( ) 

? VI-'iST 

! SI21^FREE( ) 
!SI54-FRD.E< ) 
!SI87=^^FRE£( ) 

1 SI01=^FREF( ) 

!2I1 :^FRtF() 

!SI2 ^-FRt:F<) 

!Bi3 -^frc:e<) 

!S14 ^FHEEC) 

!3I5 =FREE( ) 

!5I6 -FRFE() 

!Si7 --=fre:e<) 

*518 =^fre:e<) 

! SI9 ^FRFF( ) 
tSI 10“-FREE< ) 

!B1 1 1 = FREC:( ) 
•Sn2^FREE( ) 

!!£P1 ^FRCE. ( > 

!EP2 =^fre;e() 

!EP:3 =-FREE() 

*i;.P4 -*FRFE<) 

‘FPt? ^FREEO 
!l:P6 =fREtA) 

!EP7 =E'R£E() 

!£P8 =FRE.E() 
lEPV =FR£E() 
‘EP10=FREE( > 

!EP1 1-FR£E( ) 
»EP12^FREE( > 
•DETE>FRE£( ) 

% 

^LAI3£L 19fi 
I ICOU^^ICOUM 
*T£MP=DE41‘-1C0U'H 
^>JG2(TEMP; 200) 
!B1=FR£E( ) 

! ICOU=^FRE£( ) 
nEMP=FR£F< ) 

‘4> 

axor U1 
^SHGW 
^RETURN 
u 

'I 

ii(29 AD JO £41 0 0) 


< 14 AREA E41 00) 


ORIGINAL PAGE IS 
QE POOR QUAUTY 


ENDSC41 


ENDA041 


»XQT AUS 

' B1=JENT~1«25 

DEFINE Al-14 HIT E41 O O 

TABLE(NI=5, NJ=5):15 HITE ELEM "ADl.C" 1 

TRANSFERR(S0URCE=A1, ILIM-25. SBASE-‘'81") 

!B1^FREE( ) 



% MATRIX A 


X21 ---IjS 

3. 

‘\JEMT‘L 

1 

( 14 

ELRF 

E41 

0 

0) 

X3i =DS 


" JENT", 

1 

(14 

ELRF 

£41 

0 

0) 

X32^DS 

6. 

" JEI-4T‘L 

1 

( 14 

ELRF 

E41 

0 

0) 

X41-=DS 

7. 

JENT'L 

1 

( 14 

ELRF 

E41 

0 

0) 

X42'-=DS 

8. 

'•vJENT'L 

1 

( 14 

ELRF 

E41 

0 

0) 


‘LN22-=-l, 0-4-X:32/X21 
!EN24=^X32/X21 
!i:N:32-X3i/X21‘-i. 0 
!EN34^ -X3i/X2t 
‘ EN42^--1 . 0“;^X42/X21 
!EN44=X42/X21 
•ENb2 = X41/X21“l. 0 
1 EN54=-X41 /X21 
JtXGr U 1 


i^TI 

( 1^3 AMAl 

E41 "/^DLC" 0) 


" 1 . 0 

-1. 0 

0. 0 

-1. 0 

o 

o 

0 0 

“EN22‘ 

“EN32" 

“EN42" 

‘■EN52 

1 0 

O 

o 

0. 0 

0. 0 

0. 0 

0. 0 

■•EN24' 

' “EN34" 

'‘FN44" 

“EN54 

0. 0 

i. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

1. 0 

0. 0 

0. 0 

0 0 

0 0 

0. 0 

1. 0 

0. 0 

0 0 

0. 0 

0. 0 

0. 0 

1. 0 


' X 21 ^FREE( ) 

I X31-FREE< > 

‘ X32-FREE< ) 

I X41^FREE( > 

‘ :<42-FHEE( ) 

!t:N 22 ^FREE( ) 

‘EN24-FREE( ) 

‘EN32=FRE:E( ) 

‘EN34-FREE( ) 

!EM42=FREE( > 

!r:T444---FREl£( ) 

>EN52^FRE£( ) 

!FN54-FRt.E( ) 

«XGT AUS 

DEFU4E SAQ=:15 HITE ELEM “ADLC“ 1 
DEFINE BAW=15 AMAT E41 "ADLC*' O 

HITF E41 “ADLC“ 0=RPR0D(BAQ ,BAW) 

DEFINE BAQ-"=14 PMAT E4 1 “JENT" 0) 

DEFINE BAW^15 HITF E41 ‘*ADLC‘' O) 

!'oACR=DS 8. "JENT“/ 1 (14 CONN E41 O O) 

!QQTH=D5 26. "SAQR‘A1 (1 SA BTAB 2 13) 

Itj FHIT E41 "ADLC‘* O^RPRDD (“DOTH" BAG, BAW) 
! OaTH=FREE( ) 

•SAGR-FREE() 

Oi 

! lOOJ =1 
LABEL 400 
! I DOF 


Sli -DS 

1.1.1 

(14 

GPST 

E41 

•*JENT“ 

0) 

SI 2 =DS 

2. 1. 1 

(14 

GPST 

E41 

" JENT" 

0) 

813 -DS 

3. 1, 1 

<14 

GPST 

E41 

"JENT" 

0) 

SI 4 -=DS 

4. 1. 1 

( 14 

GPST 

E41 

"JENT" 

0) 

SI 5 -DS 

5. 1. 1 

( 14 

GPST 

E41 

"JENT" 

0) 

GIc> ~D5 

6, 1. 1 

(14 

GPST 

E41 

"JENT" 

0) 

SI7 -D3 

7. 1. 1 

( 14 

GPST 

E41 

"JENT" 

0> 


SIB ==D5 

a, 1, 1 (14 

GPST 

E41 •' 

JENT 

" 0) 


319 ^DS 

9, 1, 1 (14 

GPST 

E41 •• 

JENT 

" 0) 


61 iO=DS 

10, 1,1 (14 

CP3T 

E41 

•■JENT“ 0) 


SI 1 1==D5 

11,1,1 (14 

GPST 

EAl 

•MENT" 0) 


3112=03 

12, 1, 1 (14 

GPST 

E41 

“OEWT" 0) 


EPl =03 

1, “ IDDF“, 1 

( 15 

PHIT 

E41 

“ADLC* 

0) 

EP2 =03 

2, “ IOQF‘^ 1 

(15 

PHIT 

E41 

"ADLC” 

0) 

EP3 :=03 

3, " IOOF“, 1 

( 15 

PHIT 

E41 

“ADLC” 

0) 

EP4 =03 

4, “ IDOF“, 1 

( 15 

PHIT 

E41 

"ADLC " 

0) 

EP5 =03 

5, “ IDOF“, 1 

< 15 

PHIT 

E4l 

"ADLC" 

0) 

EP6 =0S 

6, "IDGF’', 1 

( 15 

PHIT 

E41 

"ADLC " 

0) 

E£P7 =03 

7, “IDOF“, 1 

( 15 

PHIT 

E41 

"ADLC " 

0> 

EP8 =03 

8, “IDaF'*, 1 

( 15 

PHIT 

E41 

"ADLC" 

0) 

EP9 =03 

9, "IDaF'*, 1 

( 15 

PHIT 

E41 

“ADLC " 

0) 

EP 10=03 

10. "IDOF", 

1 (15 

» PHIT E41 

"ADLC 

" 0> 

EPl 1=03 

11, "IDOF", 

1 (15 

» PHIT E41 

"ADLC 

" 0) 

EPl 2=03 

12, “IDOF", 

1 (15 

^ PHIT E41 

"ADLC 

" 0) 


i 0*SI1-SI2 

VMI=^0. 

f VM2-2. 0«SI2-SI1 
‘v;M 2 ~ 0 . 2 ?-«VM 2 
! Vh3^^*3. 0-^SI3 
• Vh4^2. 0^SI4'-SI5 
fVh4-^0. 5-»frVM4 
! VM5-^2. 0«SI5-SI4 
‘Vh5-=0. 

'VM6--3. 0-:tSI6 
I VM7-^2. 0»SI7-SI8 
•VM7=^0. 

‘ VM8-2. 0 «-SIB-SI7 
f VM8=0. 5^VM8 
! VM9^3. 0^SI9 
IVMlO-2. Oa-SnO-SI 1 1. 

!vmi^2. o-»frsin-siio 
! VMi 1~0 .‘7-aVm A 
! Vm 2 ^- 3 . 0*8112 
$ 

!SI21^SI2*SI1 
! 8154 ^^ 815*314 
!3fB7=:3IQ*SI7 

! 8101 =-sno* 3 iii 

!SU =311*311 
! 812-312*312 
! GI3==SI3*SI3 
‘314^314*314 
!SI5-=SI5*SI5 
f 316=^316*318 
‘3I7=^--SI7*GI7 
•HT8-SIQ*SI8 
!B19^SI9*SI9 
‘8110-3110*3110 
!BI 1*1=^8111*3111 
fSI12=:3I12*SI12 


’ 813^-313*3. 0-»“SI 1+SI2--SI21 
! 313' -SI 3**0. 5 
!3I6-SI6*3 0+3I4+*SI5--3r54 


n;;r'-'Tv^ man IS 

OF POOR QUAJUTM 



! 5 

!SI9-Si9a-3. 0^■^JI7■HSIa-SIa7 
!3T9^S19^i-«0. 5 

!31 13=SI 12*3. O'^SI 10+SIl 1~SI01 
*31X2=3112**0. 

$ 

: vm =VMi/si3*EPi 

• VM2-=VH2/SI3*EP2 
! VM3-=VM3/SI3*£P3 
f VM4=VH4/SI6*£P4 
! vhej=^vris/si6*EPS 
f V116=VM6/SI6*e:Pr. 

! Vh 7 = VM 7 / S 1 9*EP / 
f VM8-=VHi-3/3I9*EP8 
» v;m 9 *^VHv/SI 9*EP9 

! VPi 1 0- VM 1 0 / S n 2*EP 1 0 

• vpin =vm i/3ii2*EPii 
! V f 1 1 2 V P I i 2 / G I 1 2 > E P i 2 


‘DETE-D3 1.1,1 (14 DETE E41 ‘ 

* VM3==VH;:H VM2 + VMl*DE rE 
!DETE=^DS 2, 1, 1 (14 DETE E41 * 

* Vhs^-HVM4*DETE 

!DETE==D3 3, 1, 1 (14 DETE E41 ’ 
‘ Vh9-^VM9 I Vh8 +-VM7 « DE VE 
!DETE=DS 4,1,1 (14 DETE E41 ‘ 
! Vh 1 2 =VH 1 2'J'VM 1 1 4 VPl 1 0*DE I E 
!AREA=D5 1,"JENT",1 (14 AREA 

! VM 1 2=VM 1 2+VM9+ VH6+Vh3 
! Vril2-^Vm2/AREA 
! AREA==FRE£( ) 


JENT“ 

0) 

JENT“ 

0) 

JENT“ 

0) 

JENT" 

0) 

E41 0 

0) 


‘TEM2=IDrjP-l 
*JNZ( TEh2, 100) 
! ADl=^l. 0WM12 
*JUhP 300 
*LABEL 100 
' ]EM2=I00F-2 
n JN2 (TEH2, 110) 
! A02=l. OfiVm2 
^^ JUHP 300 


a LABEL 110 
ITEH2-=fD0F-3 

* JI'.I2(TEM2, 120) 
f AD3-=1. 0*VM12 
*JUMP 300 
*LABEL 120 

! TEM2^=IDGF-4 
*JNZ(TEH2, 130) 
! AD4==1. 0*VM12 
300 

*LADEL L30 
‘ rEN2=I0QF-5 
*xjNZ^ r£H2. 140) 

* ADS^'i. UrtVM12 
*JUhP 300 
i^LADEL 140 
ITEM2-IDQF-6 

« JNZ(TEh2, 150) 
‘ AD6=1. 0^VM12 


juhP :ioo 

150 

! T£rl2=ID0F-7 
«-JNZ(T£H2, 160) 

! AD7=1 0vVMl2 
*tJUhP 100 
»LAB£L 160 
! TEH2-r- IDGr- -8 
*JN/ ( rEh2; 170) 

‘ ADB^l 0-»VM12 
:«t.AUeL 1/0 
‘ 4 > 

^LABFL 300 
! lOOK^IDOF^l 
! TEh2 -Il>uF“ 9 
^JNZ (TEM2i 400) 
!TEM2=^FPEE( ) 

• 1D0F=^FREE( ) 
f ‘3I2l==FREE( ) 
!B154^FRE£( ) 
!SIB7=FREE( ) 
!SIOl--FREE( ) 

•SIl ^FRE£() 

1312 =FRFE() 

‘SI3 -FREE() 

I SI 4 =FREF.() 

“315 ^FR6:F<) 

•H16 -FRE£() 

!SI7 -FREF() 

•SIS -FREEO 
‘S19 -FREEO 
!SI 10=^FREE( ) 
!SI11=FRE£( ) 
•SI12=FREE() 

•EPl ^FREEO 
•EP2 =FREE() 

!fiP3 ^FREE() 

!EP4 =FREE() 

•EP5 =FREE() 

•EP6 ^FREEO 
!EP7 ^FREEO 
•EPS =FREE() 

•EP9 =^FRE£() 

I EP10^FRE£( ) 

!EP1 i -=FR££( ) 
•EP12=FREE( ) 

• Vm -FR£E( > 
•-‘M2-FREE ( ) 

! VG3^FREE( ) 

! VH4-FRFE( ) 

‘ Vh5-FREE( ) 

■ VhA=-FRE£( ) 

! Vh7=--FREE( ) 

• Vhe-FREE< ) 

• VM9-FREE( > 
“7hl0--.FREE( ) 
f Vm 1^FREE( ) 

• VM12^FR£E( ) 
IDETE=:FREE( ) 


0RlC]r^AL PAG!*: 

OF POOR QUALr 


T 
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^ NOW ASSIGrJ LOAD TO PROPER DEGREES OF FREEDOM 
! VMS! 

DEFINE ELNO £41 O O 

rABl.£(NI=^12; NJ^12):15 ROTA ELEM “ADLC" 99 

! Al-JENT-l-i^36 


TRAN3FERR ( S0URCE^=BB2, 

ILIM=3, 

SBASE=“A1“ 

, DBASE=“B1“ ) 


' A1 =jEN T- 1 tt-36+3 
! Dl-=^12 

TRANSFERR ( S0URCE^BB2, 

ILIM=3, 

SBASE=“A1” 

, DBASE='‘B1" ) 


! Al-viEN r- 1^36 + 6 
' D1^24 

TRANSFERR ( S0URCE^BB2, 

ILIM=3i 

SBASE='*A1" 

, DBASE="B1“ ) 


1 A1 ^JENT-1-^^36+9 
‘131=36+3 

TRANSFERR ( S0URCE=BB2, 

ILIM=3, 

SBASE=**A1" 

, DBASE=”B1“ ) 


‘ A1=JENT-1»36+12 
‘ D 1=46-+ 3 

TRANSFERR (S0URCE=BB2, 

ILIM=3, 

SBASE=“A1" 

, DBASE=“B1“ ) 


* A i =\JENT— 1 *36+1 5 
‘Dl=60+3 

TRANSFERR (S0URCE=BB2, 

ILIM=3, 

SBASE=‘*A1" 

, DBASE=“B1“ ) 


! A l-=OENl -1*36+18 
! B1 =72+6 

TRANSFERR ( S0URCE=BB2, 

ILIM=3, 

SBASE=*'Al“ 

, DBASE='‘B1“ ) 

• 

1 A1=JENT-1*36+21 
IB 1=84+6 

TRANSFERR ( S0URCE=BB2, 

ILIM=3, 

SBASE=“A1“ 

, DBASE="B1’' ) 


! A 1---JENT- 1*36+24 
1 B 1=96+6 

TRANSFERR ( S0URCE=BB2, 

ILIM=3, 

SBASE="A1*’ 

, DBASE="B1“ ) 


! Al=JENT-l*36+27 
*01=108+9 

TRANSFERR ( S0URCE=BB2, 

ILIM=3, 

SBASE=*‘A1" 

, DBASE=“B1“ ) 


! Al=JENT-l*36+30 
•01=120+9 

TRANSFERR (S0URCE=BB2, 

ILIM=3, 

SBASE="A1“ 

, DBASE='‘B1“ ) 

- 

•Al=vJENT-l *36+33 
1 U 1=132+9 

IRANSFERR ( S0URCE=BB2, 

ILIM=3, 

SBASE=*'A1" 

1 DBASE="B1'* ) 

• 


‘ A1=FREE< > 

! ni^FREE( ) 

TABLE(Ni.^l2. NJ^1):15 ELLO E41 "ADLC’’ 1 

J^l: ”A01“ “AD2” 0.0 “AD3" '*AD4'* 0.0 “AD5‘* “AD6“ 0.0 "AD7" “AD8" 0.0 

DEFINE A1 = 15 ROTA ELEM “ADLC" 99 
15 ROTA ELEM "ADLC" 1 '-'=RTRAN< A1 ) 

DEFINE DB2= 15 ELLO E41 “ADLC” 1 
DEFINE A1 = 15 ROTA ELEM "ADLC* 1 
14 ADLO VE41 "ADLC** 0 = RPROD ( Al, BB2) 

‘ADl^DS 1, li 1 (14 ADLQ VE41 "ADLC” 0) 

!AD2^DS 2,1,1 (14 ADLO VE41 “ADLC" 0) 

:AD3==DS 3,1,1 (14 ADLO VE41 "ADLC" O) 

•AD4^DS 4,1,1 (14 ADLO VE41 “ADLC” O) 

‘AD5 = DS 5,1,1 (14 ADLO VE41 "ADLC** O) 

*aD6--=DS 6,1,1 (14 ADLO VE41 »ADLC“ O) 

!AD7-.DS 7,1,1 (14 ADLO VE41 "ADLC* 0) 

!AD8-^DB 0,1,1 (14 ADLO VE41 “ADLC” O) 

•AD9==^^DS 9,1,1 (14 ADLO VE41 "ADLC’* O) 

'AD10=DS 10,1,1 <14 ADLO VE41 *'ADLC" 0) 
lADli^DS 11,1,1 (14 ADLO VE41 “ADLC” O) 




! AD .1 2--DS 
ip 

! -Jl-OS 
! J3-DS 
' J3=D3 
I J4-DS 

svsvec ■ 


i; 


1-1 (X4 ADLO VE41 "ADLC" 


O) 


13, "JENT-. 1 

14, "JENT", 1 
13- "JENT", 1 
16, "JEMT", 1 

APEL FORC 


I-l 


sl-"Ul " 

: "ADI" 

U-" J2“ 

; "AD4“ 

U~"U3 " 

; “AD7" 

U-"U4" 

“ADIO" 

1-2 

U=“U1" 

“AD2" 

J- " J2 " 

“ADS" 

J-" J.3- 

"ADO ■' 

'J=“ J4" 

"Aini" 


j-'-.ji- . 

“AI)3” 

U-“..>2" : 

'*AI)6*‘ 


“Al)9“ 


‘Vvf)12" 

'AD1--FREE< ) 


( 14 CONN E41 0 O) 
<14 CONN E41 0 0) 
<14 CONN E41 0 O) 
<14 CONN E41 0 0) 
“ADLC" 1 


'AD2=FRf-;E{ ) 
AD3 ^FRfE( ) 
AD4-FREE( > 
AD5^FREE( ) 
.hD6=FREE( ) 
AD7--FR£E( > 
A08-FREfc:< ) 
AD9=FREE( ) 
AD10=FREE( ) 
AD I 1 -=FREE { ) 
AD12=FREE< ) 
•Ji -=FREE( ) 
•J2--FREE( ) 

33- --FHEE( ) 

34- FREE( ) 


*XQf DCU 

erase 15 
■«return 

ih 


ORTGIMSC mns W 

OF POOR QUALITY 


% hND AD41 

-+(29 PREP E42 0 0) 

ENDPE42 

I this subroutine is for element TYP E42 

* SpORi^'LI^BjARY^fJ^ IS 16 

IxQT Eu Deleted at the end of this routine 

% 


extract. S0URCE=E42; 

CREATE 16 AREA E42 0 O 

'■V 


r.A I kaCT: 


create 16 ELRF E42 


content SPEC; 

CONTENT SPEC: 


<3E0M 1 $. area of ELEMENTS 

<5E0M 5, 12 » local ELEM, REF. FRAME 
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EXTRACT: S0URCE=E42: CONTENT SPEC: GEOM 22i 57 ^ ELEMENT-NODAL RELATION 

CREATE 16 ELNO E42 O O 

EXTRACT: S0URCE==:E42: CONTENT SPEC: MATE 1 » 2 MATERIAL PROPERTIES 

CREATE 16 MATP E42 O 0 
4 

EXTRACT: S0URCE=£42: CONTENT SPEC: INTE li 16 ^ CONNECTIVITY 

CREAIE 16 CONN E42 O 0 
% 

«XQI* DCU 

TQCC 16 CONN E42 0 O : TYPE==0 

^XGT Ui 

f ENT 1-0. 6220084 
•ENr2-l 0/6 O 
IENT3-0. 0446582 


*TI ( 16 

SHAP FUNC 

1 1) 


•'ENTl** 

■‘ENT2" 

•’ENT3‘* 

‘'ENT2" 

“EN't 2'* 

*‘ENT1‘' 

"ENT2'* 

•*ENT3'‘ 

'‘£NT3“ 

"ENT2‘' 

**ENT1‘‘ 

"ENT2" 

■•ENT2‘' 

'‘ENT3“ 

•’ENT2‘* 

"ENT1“ 


!ENT l ^-FREE( ) 

!ENT2=FREE( ) 

!ENT3-=FREE< ) 

% 

-^XGT AUS 
! ICOU^l 
! TEMP^O 
*:tLA3EL 200 
! TEMP=:FREF< ) 

!in-=icou- 

DEFINE A 1-^1 6 ELRF E42 O O 

TABLE (Nl-2. NJ=--4) : 17 ELEM REFT '‘ICOU** 99 

TRANSFER (SOURCE=:Ai, IL I M=8i SBASE=“ B 1 " ) 

• B1=FREE( ) 

DEFINE BQ1'-=17 ELEM REFT “ICOU*‘ 99 
16 ELEM REFE ‘'ICOU“ 0=RTRAN(BQ1) 

% COMPUTE DETERMINATES FOR INTEGRATION 
! 1L==1 

! GPl-^^O. 5 7735 
iU-ABEL 55 
f f>Xl-=-l (>i*OPl 
!PEr-=- 1. OaGPl 
? TEMP- IL " 2 
>JNZ<TEMP, 35) 

‘PX 1=^*1. O^GPl 
!PEV— 1. O-aOPl 
•^i-LAbEL 35 
‘ lEHP-IL-3 

TEMP, 37) 

•PXl— r 1. OOP I 
!PET— H OitGPl 
^LABEL 37 
! TEMP-lL-4 
^ JN2( TEMP, 39) 

‘PX C-— 1. 0-aGPl 
(PET-+1. O JiGPl 
^LABEL 39 

!DUM1--1. 0+PET>0. 25 
!DUM2— 1. O-iPXI->0. 25 
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!DUM3-=-»-l. 0*-PET^0. 25 
! DUM4=-1. 0-PXI itO. 25 
« r>UH5=-U. 0 KPEI -i^O. 25 
* 0Ur*l6--=-H . O+PXI J^O. 25 
.'DUn7=-l. O-PET^O. 25 
‘ DUr*1B= vl O-PXHJ-0. 25 

VAl3LF(NI=2. NJ=^4> : 17 HELP E42 “IL“ 99 


1 : 

"DUMl" 

"DUM2" 

J-2 : 

"0UM3" 

"0UM4" 

J~3 

"DUM5" 

"DUH6" 

J-.-4 

"DUH7" 

“DUMB" 

UEFI^fE 

£1^17 

HELP E42 "IL" 9‘ 


DEFINb. 002-^16 ELEM REFE *‘ICOU" O 
1/ JAC E42 ‘'ICOU'* ‘•IL"^-’^RPROD <E1,I3B2) 


DtJMl -05 

1. 1. 1 

( 17 

JAC 

E42 

"ICOU" 

"IL" ) 

0UM2=--DS 

1. 2. 1 

(17 

JAC 

E42 

"ICOU" 

"IL") 

DUM3*-=03 

2. 1. 1 

( 17 

JAC 

E42 

"ICOU" 

"IL") 

r)UM4*^QS 

2. 2. 1 

(17 

JAC 

E42 

"ICOU" 

"IL") 


!DUM5-^DUm^tOUH4 
! DUM6=DL/M2^DUH3 
‘DUHa=nuM5-ouri6 
! DOH8 

‘*IEMP=-'Il-l 
it JNZ< TEHP, 45) 

* Oin l^DUhB 
JUMP 53 

^ LABEL 45 
' I FMP^IL -2 
.i-JHZ ( TEMP. 47) 

! DFT2--DUMQ 
->JUMP 53 
-i^LAOEL 47 
! Tl- MP=--IL.-3 
^JtMZ(TEMP, 49) 

‘ UET3=DUM8 
X MJHP 53 
iH.AEEL 49 
‘ rEMP=IL~4 
i^JWZ( TEMP. 53) 

! D£T4'-^DUM8 
^LABEL 53 
' IL-ILH 

* rt:hP=IL-5 
*v)M2(TEMP, 55) 

rAULE(NI^4. NJ^l ) : 16 DETE E42 "ICOU" O 
‘OETl" "DET2" "DET3" •'DET4" 
!DET1-FREE( ) 

> DE r2--^FRE:E( ) 

IUET3-FREEO 
f DET4=FREE( ) 

* D(JM 1 =FR EE ( ) 

* DUh2-^FRE£( ) 

! DUM3=FREE< ) 

* DUM4=FREE( ) 

!DUM5=FREF( > 

‘DUM6^FR£E( ) 

•DUM7^FREE( ) 

!DUNB^FREE< ) 
f IL^FREE< ) 

•PXI=^FREE( ) 
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'per=^FRe:E( > 

!GP1-FREE( ) 

Dt:Fir4E SHAP FUNG 1 1 

DEFINE ELEM REFE "ICOU** 0 

17 OAUS POIN "ICGU" 0^= RPROD 


XGP1=DS 


1, 

1 

(17 

GAUS 

POIN 

"ICOU" 

0) 


YGPl-DS 

1/ 

2 1 

1 

(17 

GAUS 

POIN 

•*ICOU" 

0) 


XGP2^D3 

2. 

1, 

1 

( 17 

CAUS 

POIN 

"ICOU" 

0) 


YGP2^DS 

2, 

2. 

1 

( 17 

GAUS 

POIN 

"ICOU" 

0) 

- 

XGP3^DS 

3, 

1. 

1 

( 17 

GAUS 

POIN 

" ICOU" 

0) 


Y(jP3 - Db 

3. 

2, 

1 

(17 

GAUS 

POIN 

"ICOU" 

0) 


XGP4=0S 

4. 

li 

1 

(17 

GAUS 

POIN 

"ICOU" 

0) 


YGP4=Db 

4, 

2i 

1 

( 17 

GAUS 

POIN 

"ICOU" 

0) 



! XVQi=^XCPl^YGPi 
! XVG2=X0P2-i»*YGP2 
‘ XyG3=XGP3>»-YCP3 
! XYG4=XGP4« YGP4 

rABLE(NI = l2i NJ==1 1 ) : 16 PMAT E42 '‘ICOU“ O 


J=1 

1.0 0.0 

0. 0 

1 0 

0. 0 

0. 0 

1.0 0. 0 

0. 0 

1. 0 

0. 0 

0. 0 

J-^2; 

o 

o 

o 

0. 0 

0. 0 

1. 0 

0. 0 

0. 0 1.0 

0. 0 

0. 0 

1. 0 

0. 0 

J---3: 

0 0 0. 0 

i. 0 

0. 0 

0. 0 

1. 0 

0. 0 0. 0 

1. 0 

0. 0 

0. 0 

1. 0 

J=4 

"YQPl" 0.0 

0. 0 

"YGP2" 

0. 0 

0. 0 

"YGP3" 0. 0 

0. 0 

"YGP4" 

0. 0 

0. 0 

J ^-5: 

0.0 "XGPl" 

o 

o 

0. 0 “ 

XGP2" 

0. 0 

0. 0 "XGP3" 

0 0 

0. 0 " 

XGP4" 

0. 0 

J^6. 

"XGPl" 0.0 

o 

o 

"XGP2" 

0. 0 

0. 0 

"XGP3" 0.0 

0. 0 

"XGP4" 

0. 0 

0. 0 

J---7: 

0.0 "YGPl" 

0. 0 

0. 0 " 

YGP2" 

0. 0 

0. 0 "YGP3" 

0. 0 

0. 0 " 

YGP4" 

O. 0 

J-B: 

0. 0 0 0 " 

YGP 1 " 

0. 0 0 

. 0 " 

YGP2" 

0. 0 0. 0 " 

YGP3" 

0. 0 0 

. 0 " 

YGP4 


0. 0 0. 0 " 

XGPl" 

0.0 0 

. 0 " 

XGP2" 

0. 0 0. 0 " 

XGP3" 

0. 0 0 

. 0 " 

XGP4 

J = 1 0 : 

"XYGl " 0. 0 

0. 0 

"XYG2" 

0. 0 

0. 0 

"XYG3" 0.0 

0. 0 

"XYG4 

" 0. 0 

0. < 

1 1 ■ 

0 0 "XYGl 

" 0. 0 

0. 0 

"XYG2 

" 0. 0 

0.0 "XYQ3 

“ 0. 0 

0. 0 

"XYG4 

" 0. < 


‘ XGPl=FREE( ) 

> YGP1^=FREE( ) 

I XCP2----KR£E( ) 

* YGP2=FREE( ) 

! XGP3-FREE( ) 

* y'GP3=FREF( ) 

! XGP4^FRt;E( ) 

• YGP4=FRE£( ) 

1 XVG1=--^FREF< ) 

I :<YG2=PREE( ) 

! XYG3^FREE( ) 

* XVG4^FREE( ) 

‘NU^DS :L “ICOU’3 1 (16 NATP E42 O O) 

^NUP1=NU ^ 1. 0 
! NU^-1. Oirt NU 


TAI3LE(NI-12. NJ==12) : 16 SIEP E42 “ICOU“ 0 


1.0 *'NU‘' 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

J=2: "NU** 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

J=3: 0.0 0.0 "NUPl" 0.0 0.0 0.0 0.0 0.0 0.0 

J-4; 0. 0 O. 0 0. 0 1.0 "NU" O. 0 0. O O. 0 0. O 

J-5. 0. 0 0. 0 0. 0 "NU" 1.0 O. 0 O. 0 O. O 0. O 

J^6: 0. 0 0. 0 0 O O. 0 0. O “NUPl" O. 0 0.0 0. O 

J~7: O. 0 O. O 0. O 0. 0 0. 0 O. O 1.0 "NU" O. 0 

J=8: 0 O 0. 0 0. O O. 0 0. 0 O. 0 "NU" 1.0 0. O 

,j=9; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 "NUPl" 

J~=10: 0 0 0 0 0 O O. O 0. 0 0. O O. 0 0. O 0. O 


0 . 0 0 . 0 0.0 
0.0 0.0 0.0 
0. 0 0. O 0. 0 
0. O 0. O 0.0 
0. 0 0. O 0.0 
O. 0 O. 0 0. 0 
0.0 0 . 0 0.0 
0.0 0 . 0 0.0 
0.0 O. 0 0.0 
1. 0 "NU" O. O 
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J=ll: 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 "NU'* 1.0 0.0 

J^12: 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 *'NUP1** 

fNlMFREE( ) 
f MUP .l=FREE ( > 

*TEMP=FREE( ) 

?Dl = ICDU-l»i-U 

DEFINE A2=^l ESB E42 1 1 

TABl.F.(NI=:ll, NJ=^1):17 BBB E42 "ICOU" 0 
TRAraBFERR <S0URCE=^A2, ILIM=^1 1, SBASE=:“B 1 •* ) 

!B1=^FREE( ) 

DEFINE BB1«16 PMAT E42 ‘*ICOU“ 0 
DEFINE BB2=17 BBB E42 “ICOU** O 
!SAGR-DS 8. *'ICOU", 1 <16 CONN E42 0 O) 

!OOTH=DS 26> ’*SAGR", 1 (1 SA BTAB 2 13) 

•D0TM=^rj0TH«00TIH^6. 0 

IvS GPST E42 “ICOU” 0= RPROD < "OOTH" BB1>BB2) 

! OOTH-FREE ( ) 

!8A0R--FRFE< ) 

'4i 

! ICOU '-ICGUH 
' 7EMP=DE42-IC0U* 1 
M JO 2 (TEMP. 200) 

‘ 01--FRfciE< ) 
f ICLiU-=^FREE( ) 
f TEMr"^-FREE( ) 


aXQT DCU 
ERASE 17 
■»^RETURN 

* EWDPE42 

^ 

»^<29 T3ECA E42 O 0) ENDSC42 

a>:oj Aus 

' ICOU-l 
ITEMP^O 
^LADEL 200 
% 

•ELGR=DS 1. “ICCJU*'. I (1 ED42 REL 0 O ) 

! TEhP^ELGR- ICnV 
f EI..GR 
f 1C DM 

‘f:l.gr=^freec ) 

^^JiMZOEMP, 195) 

'4» 

?B1==IC0U-1k11 

DEFINE A2=-l ESB E42 "ADLC" 1 
TABLE(NI=^1 1. NJ=1):17 BBB E42 ’‘ICQU’* “CLC" 

TRANSFER (S0URCE=A2, ILIM=11, SBASE==“B1 “ ) 
*rJl=^FF<EE( ) 

DEFINE 331^^16 PMAT E42 "ICOU” 0 
DEFINE BB2=U7 BBB E42 “ICOU" “CLC" 

‘SAGR-^DS 8. "ICOU”. 1 (16 CONN E42 0 0) 
fOOTH=^DS 26. *‘SAGR'\ 1 (1 SA BTAB 2 13) 

•00TH^00TH*n0ni«-6 O 



17 QPST E£42 ’‘ICOU” "CLC*’ ^ RPROD ( ‘*OOTH’' BBliQD2) 115 

!OQTH^ FREE‘< ) 

!SAGR-^FRFE( ) 

^ ♦ 

DEFINE SIEP E42 ’'ICOU" 0 

DEFINE B12=a7 GPST E42 “ICOU“ “CLC’* 

16 GPEP E42 “ICOU” ”CLC“=^ RPROD (Bll.tn2) 


Git =DG 

1,1,1 

( 16 

GPST 

E42 

"ICQU" 0) 


GI2 ^DS 

2, 1, 1 

( 16 

GPST 

E42 

“ICOU” 0) 


SI3 ^DS 

3, 1, 1 

( 16 

GPST 

E42 

“ICOU“ 0) 


SI4 -^DS 

4,1,1 

( 16 

GPST 

E42 

“ICOU“ 0) 


S15 ^DS 

5, 1, 1 

( 16 

GPST 

E42 

“ICOU“ 0) 


GI6 =^DS 

6,1,1 

( 16 

GPST 

E42 

“ICOU" 0) 


S17 =-03 

/, 1, 1 

( 16 

GPST 

E12 

“ICOU” 0) 


sia =^03 

8, 1, i 

( 16 

GPST 

E42 

“ICOU” 0) 


SI9 =03 

9, 1, 1 

( 16 

GPST 

E42 

“ICOU” 0) 


SI 10=DS 

10, 1, 

1 ( 16 

1 GPST 

E42 

“ICOU” 0) 


31 1 1=DS 

11,1, 

1 ( 16 

► GPST 

' E42 

“ICOU” 0) 


SI 12=DS 

12, 1, 

1 < 16 

^ GPST 

E42 

“ICOU” 0) 


EPl =DS 

1,1,1 

( 16 

GPEP 

E42 

“ICOU” “CLC” 

) 

EP2 =DS 

2, 1, 1 

( 16 

GPEP 

E42 

“ICOU” "CLC” 

) 

EP3 =DS 

3, 1, 1 

< 16 

GPEP 

E42 

“ICOU” “CLC” 

) 

EP4 =DS 

4, 1, 1 

(16 

GPEP 

E42 

“ICOU” “CLC” 

) 

EPS ^DS 

5, 1, 1 

( 16 

GPEP 

E42 

“ICOU” “CLC” 

) 

EP6 =D3 

6, 1, 1 

( 16 

GPEP 

E42 

“ICOU” "CLC" 

) 

EP7 =DS 

7,1,1 

( 16 

GPEP 

E42 

**ICGU“ “CLC” 

) 

EP3 =DS 

G, 1, 1 

( 16 

GPEP 

E42 

“ICOU” “CLC” 

) 

EP9 =DS 

9, 1, 1 

( 1 6 

GPEP 

E42 

’*ICOU“ “CLC” 

) 

EP10=DS 

10, 1, 

1 ( 16 

► GPEP 

E42 

“ICOU” “CLC” 

) 

EPl 1=1)3 

11, 1, 

1 (16 GPEP 

' E42 

“ICOU” "CLC** 

) 

EP 1 2= DS 

12, 1. 

1 (16 GPEP 

' E42 

“ICOU” “CLC” 

) 


GI3=“1. 0#BI3 
EP3=-1. 0frEP3 
SI6^=-1. 0-«'SI6 
EP6^=-1. 0^EP6 
0*SI9 
EP9=-1. 0*EP9 
STl2=-i. o*«-sn2 
EP12^-“- 1. 0j^EP12 

EP l^^SI l-tEPl 
FP2*-=SI2^EP2 
LTP3=2. QJi-SI3^t£P3 

DEIE^DS 1,1>1 (16 DETE E42 “rCOU” 0) 

E P3^EP 1 4 EP2+EP3«'DETE 

FP4-SI4<iEH4 

EP5=SI5^rEPt-> 

GP6-=2. 04fSI6*lZP6 

DETE=DB 2>1>1 (16 DETE E42 **ICOU’* 0) 

EP6=^EP4+*EPf;^EP6*DETE 

EP7=SI7*FP7 

EP8=SIt3ic£PB 

EP9=^=2. 0i^SI9*h:p9 

DETE=^OS 3; 1,1 (16 DETE E42 "ICOU’* 0) 

£P9^=EP74.£P8+EP9^<•DETE 

EP10==5I 10«EP10 

EPil~G( llii-EPll 

f:;Pl2=2 0i^SI12*EP12 

DETE-'-'DS 4,1,1 (16 DETE E42 ’’ICOU" O) 


LL6 


! EP 1 2- EP 1 0 ^EP 1 i ^ EP 1 2*DETE 

*PSIS=^1. OaO. 0 

?P3I3 

' PSiri--EP3 + EP6+EP9+EP12 
!NU=^OS l,'‘ICGU“.l (16 MATP E42 0 0) 

^ PSIS- • P&IS/MU 
! PSI^PSI+PSIS 


! ICUU 

!riU=^FREE( ) 
!BAGR^PREE( ) 
?OaTH'-:PREE( ) 

$ 

fSI21-SI2^tSIl 
f SI54--SI5*S14 
•3IS7^SI3^^SI7 
! 3101=^^31 10*SI 1 1 

!Sn^SIl*3Il 
! SI2==^Sl2^i-SI2 

! SI4-^-SI4«-SI4 

' b I 6-^316*3 1 6 
* SI7^-SI7^<SI7 

•siB-sia»-sia 

!3IV^SI9»SI9 
' 31 10--^51 101*3110 
*31 11=^31 1 1 «SI 1 1 
! SI 12^SI12«*S1 12 
h 
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'SI3=^S13*3. 0+b'1 1+SI2-SI21 

'l.»ETF.=DS 1,1,1 <16 DETE E42 "ICOU" 0) 

'SI3---SI3»»0. StfDETE 

' S16-‘SI6*3. 0+SI4+SI5-SI54 

!nETE=DS 2, 1, 1 <16 DETE E42 "ICOU" 0) 

!SI6-S16«»0. 5*DEfE 

>SI9^SI9*3. 0+SI7+SI8-SI87 

'DETE---DB 3,1,1 <16 DETE E-12 "ICOU" 0) 

'SI9-8I9-t»frO. 5<^DETE 

' SI 12-SI 12»3. O+SI lOfSIl 1-SIOl 

'DE7E==DS 4,1,1 <16 DETE E42 "ICOU" O) 

‘SI 12-SI 12*«0. t>*DETE 

' Vi'IS r-SI3 t SI6+SI9+SI 12 

IAREA=DS 1, "ICOU", 1 <16 AREA E42 0 0) 

! VMST-VriST/AREA 
! Vi'IST 

'TEnP^CLC-CC0M-CDIS-CS21-CS41-l 
«-ULZ<TEMP, lf>f5) 

! TEMP^TEMP+1-CS42 
^JOZ(TE!*fP, 15f>) 

* 1 EhP-CLC -CCCM-CDI S-CS21 “CS4 1 

! JENT=^03 *‘1EMP'*, 1 .1 <1 ST42 LIST 1 1) 

f JL-.NT 

•7EhP-=lC0U JENT 

* JEN7=^FREE( > 

«JNZ(TEHP. lf?5) 
fc 

*SAOF^==DS 8, "ICOU'b 1 (16 CONN E42 0 0> 

■OOTH=-DS 26* ‘‘SAGR'b 1 (1 SA BTAB 2 13) 

‘ PS I 


tJ* tf 
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!P:r:iS 

! PS1S'-^VMS1>00TH+PSIS 
i PSI^VhSTfrODTHHPSI 
I PS I 
I PS IS 

! BAGR^FRPE< ) 
*OrjTH==FR£E< ) 
f DELT^^FREE ( ) 

‘4» 

•a.AbEL 155 
! TFMP=FRFE< ) 

>PBIS 

!PSIS--^FREE< ) 

! ) 
!SI54*-::FRFE( ) 
!3IS7=FRFE< ) 
‘SI01*^FRP:E( ) 

‘ f£m=:--FRtE< ) 

‘SI I -^rREI-<) 

‘SI 2 -PREE() 

‘SIS ~FRFE() 

*SI4 ^^FREE( ) 

!SI5 =FREE() 

*SI6 --^FRESO 
‘SI7 r.FREE() 

‘ S18 ^FRFE< ) 
fSI9 ^-FREE() 

‘ S I io=^frl-:e< ) 

‘SI 1 l^'FRffE( ) 

•Si 12=FRFE< ) 
lEPl ^^FREFK) 

!FP2 ^^FREE() 

!£P3 =FREE(> 

*EP4 =FREE<> 

TPS ^FREEO 
•FP6 ^FREE() 

•EP7 =:^FRPE() 

'FP8 ^-FRFE() 

•FP9 =^FRcE() 

*t:P10 = FREE< ) 

!EP1 1^FREE< ) 
‘EPIS^FRPE ( ) 
l OETE^^FREEC ) 

! AREA==^FRFE< ) 

■«LABEL 195 
' ICGU=--ICaUt l 
! fEr'1P=DE42-IC0U+l 
«JG7<TEMP, 200) 
'IU.--FREE( ) 

‘ ICOU^FREEC ) 

! TEiMP:=FREE< ) 

4- 

>A0f Ui 

1-^SHOU 

hpETURN 


EIMDSC42 


ENDA042 


-(29 AD JO LI42 0 O) 

•4* 

»iXGT U1 

‘OP 1=0. S773t5i^0 5 
! JL-NT 

?AO=US 2,1/1 (.16 ELFh REFE "JEIMT" 0) 
•BO=DS 4> 2, 1 <16 ELEN REFE ’'JEMT" 0) 
‘ AOO=AO*AO 
! BOO=-DO*RO 
! AO 1=2. O/AO 
‘ U01=2 0/BO 

• mBI-AO^^BO 

' AB0=2. 0/AlJl 

• AB1--*FREE( > 

•fi 

‘ H0U=0 

’ Pa 1=0. 5 -CPI 
•}'Er=0. 

K label 55 

‘A 123=12. 0«PXI"6. 0 
' A2=l . O-PE T 
‘ 0 I 1--'A123<A2/A00 
' IJ12-0. 0 

‘D13=-6. OJ^PXX i'4 0 kA2/A0 

• 014=~lil I 
' B 1 5=0 0 

:016=-6. 0*PXI+2. 0«A2/A0 
•1317= -2. 0*PXI + 1. OfrPEI O/AOO 

• BiS=0. 0 

IB19 = -3. 0*PXr +-l. 0^^PET■«■2. O/AO 
!B1 10=12. 0*HXI-6. 0<<-Fl.-T/A00 

• B1 1 1=0. O 

' B1 12=-3 O^fPXI f-2. O^^PETHv2 O/AO 
‘ RnW=ROW^ 1 

ir{I(17 AROW E42 “ROW" O) 

"-U i 1 " 

"B12" 

*bl3" 

"B14" 

"B15" 

“B16" 

' u 1 / " 

•'B19" 

•* I i 1 t O " 

"HI I i " 

"Bl 12" 

•A 123= 12. O-w-PET-6. O 

• A2=.l. 0-PXl 

‘01 l=A123»A2/B00 
•ul2-=6. O^^PEl -4. 0^iA2/U0 

• Birj- O. O 

‘014=2. 0«PET-1. 0^t-PXl4i‘6. O/BOO 
1015=3 O^^-PET-2. 0*PXIfr2. O/DO 
! i;16=0. 0 
’P17---B14 

!Bi8=3. C»^P£T-1. 0^<*PXHi2. 0/00 
ML19--0 0 

‘ rn 1 -B 1 1 
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: 0111^^6 0«-PE (--2. 0i^A2/150 
! n l 12^0. o 

.‘r?ow^vFiow-t.i 

•i<ri(|7 AHQU E42 “RlDW" 0) 

”\:<i i.“ 

'•h i2" 

13" 

•3U4" 

•*U I 7" 

"0.1 8" 

"019" 

"0110" 

":m 1 " 

" ti .1 i 2 " 

’ A 123^ 4 0-w-HLT 
« A2=^4. 0-sPX} 

•A,3^PEr-i. OiiPET'1^6. 0 
! A4=-PX I-l 0 kPX1*6. 0 
*nU-^-i 0-A3-A4*A130 

* 012=^“PET^^PEH<3 0+-A12Z-t-l. 0i*-A01 
‘ iU3~PX:i-»PXii<-3. 0 -A2+1. 0«-D01 

‘ H .1 --B J 1 

‘ulS==--Bi;? 

•016^3 0«-PXI-2 O^rPXIttOOl 
*017^-014 

‘B18-.3. Oi^PET-2. O^PETK AOl 

* LU9=-B1 A 
mil 0^13 1 4 

! 0 1 1 1 - -13 1 B 
! B 112^-013 
‘r<OW-:ROUM 

^ll(X7 AROW L-4,? "ROW" 0) 

"Bit " 

"B12" 

"013" 

•‘014" 

"Ul5" 

"0 16" 

" 0 1 / " 

".ine" 

"0 19" 

"0 t 10" 

"0111" 

"0112" 

' 7 {; NP=^R0W*-3 
^ OOZU EhP. 67) 

•Pa 1=^0. f) + C;Pl 
■PET^O. f>-CPl 
iirJUMP 77 
Ja.AOEL 67 

* TE'MP:rr<OW-o 

^ ■Jv\2 • lOrlP. 69) 

M'XI- 0. V)*K)Pi 
'PEf^'O n+CPl 
k.JUMP 77 
^.‘LAOEL 69 
! 3 EMP^RDir;-9 
TfEMP, 77) 
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'PXI~0. 5-ttf'l 
!f'K.-T=^0. 

77 

• U-MP'-RDW-l^ 
kJH2<TC-:HP. 


! r<OW-HHtE< > 

! ftiyj— FREf^ ^ ^ 

! (C'-FREE ( ) 

' (\:i-FR£E< ) 

' A4-rr<FE( ) 

‘ AO= FI<fc'E( ) 
■HCi-FRFF( ) 

! AOO-f-RL-E( ) 
'BOO=FRK F.< ) 

' A01--=FREE( ) 

! B01“FREE< > 

' ABO-FRt"E^( ) 

' OP 1 -FREE ( ) 

! P X I ----FREE ( ) 

' pet-free < ) 

> B1 1=FREE( ) 

I in 2= FREE ( > 

• ii 13 -FREE! ) 

' I(14-FREE < ) 
■n 15~FREE( ) 

' li 16-FRtE ( ) 
■Bl7-r-REE( ) 

' li 13=FREE ( ) 
nn9~FREE ( ) 

' m 1 0-- FREE < ) 

' in 1 i-=FRtE ( ) 

' tf 1 12-^ FREE! ) 
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'FMOD-DS l,”JENT".l (16 MATP E42 O 0) 

'rlUt.-D5 2, "JENT". 1 <16 MATP E42 0 0) 

!NUE2-NUEftNUE 

U4UE2-1. 0-MUE2 

'El 1--=EW0D/NUE2 

'E12=E11^^NUE 


IMUE2-1 OvNUE 
■t 33---0. 5K-EM0D/NUE2 


• ti'lUD-FREE ( ) 

'MU£ -Fm-E< ) 

! r lUE2'--FR£E ( ) 
<tTlU7 DMAT E42 


Ei 1 " 

“E12" 

0 

. 0 

El 


“E 

11“ 

0 

0 

(> 

0 

0 

1. 0 

■*E33 

ij 

0 

0 

0 

0. 

0 

i 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0, 

0 

0. 

0 

0. 

0 

0. 

0 

0 

0 

0. 

0 

0. 

o 

0. 

0 

0. 

0 

!E 1 

=FREf: ( ) 



AIJI-C" 0) 


0. 0 

0. 

0 0. 0 

0. 0 

0. 

0 0. 0 

0 0 

0. 

0 0. 0 

“EU‘* 

"E12" 0.0 

“tri2“ 

"Ell 0. 0 

0. 0 

0. 

0 “E33 

0. 0 

0. 0 

0, 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 


0. 0 


0. 0 

0. 0 

0. 0 


0. 0 

0. 0 

0. 0 


0. 0 

0. 0 

0. 

0 

0. 0 

0 0 

0. 

0 

0. 0 

0. 0 

0. 

0 

0. 0 

0. 0 

*‘E11“ 

*' 

'E12“ 

0. 0 

*'E12“ 

II 

‘Ell” 

0. 0 

0. 0 


0. 0 * 

•E33" 

0. 0 

0. 

0 0. 

0 

0. 0 

0. 

0 0. 

0 

0. 0 

0. 

0 0. 

0 


0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. o 

0. 0 

0. 0 

0. 0 

0 0 

0- 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0. 0 

0, 0 

0. 0 

0. 0 

0. 0 

0. 0 

Ell” 

..£12” 

0. 0 

E12” 

”E11” 

0. 0 

0. 0 

0. 0 

“E33” 


! Eia'^FREEl ) 
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‘K33-fri:£e:( ) 

h 

X O F AU3 

OFF INF Al = 17 ARO‘^J (142 1 0 

DEFINE A?-^0 7 AROW E42 2 0 

DEF- INF A3^*17 ARON E42 3 0 

DEFUjE A^^'-17 ARQW EA2 4 O 

1>£FINE AS^'i7 Af<OW E42 t5 O 

DEFINE A6^17 AROW E^2 6 O 

DEFINE A7^17 AROW E42 7 0 

DEFINE A8=^17 AROU E42 8 O 

DKFIhiE A9^17 ARON E^2 9 O 

OEFH.1E A10--17 AROW E42 10 0 

UEr INE /M1^17 AROW E42 11 0 

DEFINE A12=17 AROW E42 12 0 

TABl.E(NI-12> NJ^12):16 BE E42 “ADLC" 99 

TRANSFER <SOUfKE=-“Al. ILIM^12, DBASE^O) 

T R ANSFER < SnURCE=A2/ I L I N= 1 2> DBASE-^ 1 2 ) 
TRANSFER (S0URCE-^A3, ILIM=12> DBASE-24) 
rRANSFER(S0URCE=^=A4, ILIM==^12, DBASE-*36) 
TRANSFER < SOURCE*>AS. I L I M=^ 12, DB ASE=48 ) 
•|RANSFER<SaURCE^A6, ILIN:^12, DBASE=60) 
TRANSIT ER ( S0URCE=^A7, I L I M=-* 1 2, DB ASE^72 ) 
TRANSFER (SGURCE^^AS, ILIM=12, DBASE-84) 

T R AN3FEF< ( SOURCE^ A9 > I L I 1 2 , DB ASE-96 ) 
TFiAN3FF£R (SOURCE^AIO, ILIM=‘-12, DBASE^^lOB) 
TRANSFER ( SOURCE= A1 1 > IL I 1 2, DBASE^lOO ) 
i'F?ANSFER (S0URCE^A12. ILIN^12, DBAS£=^132) 
DEFINE A2-16 BB E<^2 "ADLC" 99 
IS Bi3 E42 "ADLC" 0=RTRAN(A2) 

DEFINE A 1=^1 7 DMAT E42 *'ADLC“ 0 
DEFINE A2=^16 DB E42 '*ADLC* 0 
17 DB E42 “ADLC" 0=-’RPR0D ( A1 , A2 ) 

' ID0K = 1 
kLABEL 400 

'SIl =^DS 1,1,1 (16 GPST E42 "JENT” O) 

‘EF2 =US 2. 1> 1 (16 GPST E42 “JENT” O) 

'Si3 =^DS 3>1,1 (16 GPST E42 “JENT” 0) 

•314 *DS 4,1,1 (16 GPST E42 ‘*JENT“ 0) 

!3ID --^D3 5,1,1 (16 GPST E42 “JENT” O) 

‘SI6 -DS 6, 1,1 <16 GPST E42 “JENT‘* O) 

■S17 =^^DS 7,1,1 (16 GPST E42 “JENT** O) 

‘318 ^DS 3, 1,1 (16 GPST E42 “JENT” 0) 

■S19 -DS 9, 1,1 (16 GPST E42 "JENT“ 0) 

•SI10=DS 10,1,1 (16 GPST E42 “JENT" O) 


SI 1 l=DG 

11, 

1,1 (16 

GPST 

E4 

2 "JENT" 0) 


SI 12=DS 

12, 

1,1 (16 

GPST 

E4 

2 "JENT" 0) 


EPl =-'03 

1, •' 

IDOF", 1 

( 17 

DB 

E42 

"ADLC" 

0) 

EP2 

2, “ 

IDOF", 1 

( 17 

DB 

£42 

"ADLC" 

0) 

IIP 3 ^DS 

3, •• 

IDOF“, 1 

(17 

DB 

E42 

"ADLC" 

0) 

EFT =^DS 

4, “ 

IDOF", 1 

( 17 

DB 

E42 

"ADLC" 

0) 

EPS =^DS 

5, '• 

I DOF", 1 

(17 

DB 

E42 

"ADLC" 

0) 

lP6 =DS 

6, “ 

IDQF", 1 

( 17 

DB 

£42 

"ADLC" 

0) 

r:P7 

7, " I OOF", 1 

( 17 

DB 

E42 

"ADLC" 

0) 

EPS =^DS 

8, " 

I DOF", 1 

( 17 

DB 

E42 

"ADLC" 

0) 

EP9 =DS 

9, 

I DOF", 1 

( 17 

DB 

E42 

"ADLC" 

0) 

t_P lO^DS 

10, 

"I DOF", 

1 ( 17 

DB 

E42 

"ADLC" 

0 

FPl l=DS 

11, 

"IDOF", 

1 ( 17 

DB 

E42 

"ADLC " 

0 

EP 12= OS 

12, 

"I DOF" 

, 1 ( 17 DB E42 "ADLC 

1 


122 


OifSl3 
•Sr6=~l. 0^316 
!BI9--^-l. 0*S19 
‘31 12= -1. OfrGI 12 
•00S2^-1. 0/2. 0 
‘ VMl-2. OJtSI 1-SI2 
' vm --^oos2»vni 
!VM2-*2 0^t-BI2-Sn 
’ VH 2 =^- 0 QB 2 tt-VH 2 
!Vrl3= 6 0^:-S I 3^^-0082 
! VM4^2. 0*BI4-Slt? 
f Vi14^0(JB2*VM4 
.‘VMS>^2 0>SIb-SI4 
! VHt 3 ^DOS 2 »^Vi' 1 S 
!Vi'l6- 6. 0^816*0082 
: VM 7 ^ 2 . 0 *SI 7 -Sia 

‘ VMS=^=2 0*SIB~2I7 
! VI IB^^0US2*VHB 
!VI‘1V^6 o*siv*ons2 
IVHtv:)-2. 0«^SI10- 8I11 

: Vmo=GQS2*VMLO 
•VHl J.-=2. 0*SI11~SI10 
! Vhl I ^OOB2*Vm 1 
'VN 12=^6. 0*3112*0082 
•0DS2=FREE< ) 
fc 

iSl 21 =^S 12 *SI 1 
I BI!>4=^SI5*SI4 
'3I07---SIQ*SI7 
‘3I01=SI10*SI11 

4 » 

!Si 1^^SI1*SI1 

'312=^812*512 

•3f3--*SI3*3I3 

•SI4-SI4*SI4 

'31^-SI5«SIt3 

:3r6-S16*SI6 

'BI7-SI7«SI7 

•3l8=^SIBfrSI8 

'SI9-S19*SI9 

'SI 10-=GU0*SI10 

"313 1-3111*5111 

"Ui:>-GI 12 tfSI 12 

H* 

•3i:i*--S13*3. 0+SI1+SI2-S121 
'513-^513**0. 5 
'SI6*-SI6*3. 0+SI4+SI5-SI54 
! S I 6--S T 6**0 5 
‘SIV--SI9*3. 0+SI7+SI8-SI87 
' 519- 519**0. 5 

!Sn2=SI12*3. O+SIlO+SIll-SIOl 
! 31 12^51 12**0. 5 
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!VMi = vm/si 3 
! VM 2 ~ VM 2 /aKJ 
! VH3=Vtl3/SI3 
) VM4=VM4/!51 A 
' VM5=VM5/Sr.(. 



! VM6 -- Vh6/lSi6 
I VM*/- Vri7/S l V 
» VMB-VilO/SIV 
' VM9^-VM9/3IV 
' Vm0-=^VH10/SU2 
‘ vm 1 =-VH 1 l/S L 12 
! vm?'-^vm2/si J2 

I vMi-^vm M-f-Hi 
? VhZ'^VnP-rrfc'P^ 

♦ VM3-Vh3^EP:i 
» Vn4==VH4>EP4 
f Vht>=VMtj*EPt> 


‘ VM6^-=Vn6*E‘P6 
! '/M7-Vli7^^iVF7 
! VMB-vMIS+i-EPB 
‘ VMV=^Vri9iiEP9 
! vmo=^vmotti:pio 
' vhi i = vmu^EPn 

* VH12==Vf112fr£PI2 
\DE‘lh=OS 1, 1. 1 (16 DETE E4i 
? VM3=^ Vti3+VH2 ^-VM 1 4^DETE 
‘DF.TE^^Db 2. 1, I (16 DETE E42 *’JENT" O) 
! V IM 6 ‘ M 6 V hi + VM 4 * D E T E 
IDETF--DS 3.. i> 1 (16 DETE E4; 

! Vh9=^'VM9+VH8 +Vn7^^DFTE 
•PETE^DB 4>1.1 (16 DETE E42 *‘JENT' 
f VN 1 2=^ VM 1 2 + Vmi +VH 1 0 ^DETE 
‘AREA=^U3 1,“JENT">1 (16 AREA E42 0 0) 

•SA0R=^DS3 8>"JENT'\1 (16 CCJl 

'00TM=D3 26. "SAGR’'> 1 (1 SA BTAB 

' VM 1 2--^ VH 1 2 ^ VH9 4-Vl‘16+VN3 

* VM12=VH12/AREA 


JEMT“ 

0) 

JENT" 

0) 

JENT” 

0) 

JENT" 

0) 

E42 0 

0) 

E42 0 

0) 


13) 


!Vm2- (3. t>»VH12 
12/DOTH 
! AREA^EREE( ) 

. SAGR=^FREE( > 
!OUTH^^FREE( ) 




• rEh12=^U>0F-l 
^MZ ( TEH2. 100) 

• AD1--1. OuVrilS 
X- JUMP JOG 

^ LABEL 100 
! rEM2==^lI)0r -2 
n‘EM 2 . 110 ) 

• Ai'/L'-^l. 0«VM12 
»vK.JMP 300 

LABEL 110 
‘ r£M2-.ir)or-3 
*-J»MZ ( TEH2, 120) 
' A03==l. 0>UM12 
frJUMP 300 
*«L.ABEL 120 
■ TEM2>-^1 DOF-4 

TEM2. 130) 
‘ AD4^1. 0-K VM12 
{^JUHP 300 
^LABEL 130 
• TEM2-IDDF-5 
^^JN2(TEM2> 140) 


'Ans-1. 

>v^.iUhP 300 
‘•I ABEL 140 
! l Eh 2 -I 00 r-O 
*^JNZ (TEH2> lt:50) 
•AD6^*1. Or VH.12 
*JOMP 300 
ABEL. ibO 

‘ rEK2^^aD0F 7 

160 ) 

‘ A 0 7 ^ J o«vm 2 
i^JUhP 300 
LABEL 

• rEM 2 ^ID 0 K -S 
AjrJ2CIEH2. 170) 

! A 0 S -:1 0 KVHI 2 
;n^UMP 300 

i. LABEL 170 
‘ ] EM2--ID0E -9 
JrJMZ<TEH2> 100 ) 
•A09=l 0KVM12 
J^JUriP 300 
'H.ABEL 130 

• rCM 2 =^ai)C)!“ -10 
< vJi4Z ( TEHLL 190) 

• A010^=-1. OJtvmS 
rvJUHP 300 

f LABEL IVO 
‘ ( En2 =inof--i 1 
( 1E32, LOO) 

! ADI L- 1. 0« VM12 
i^rJUMP 300 
^il.ADEL 200 
‘ lEM2^=IDOP-12 
^ JM7(TEri2. 210) 

! AD I 2=1. 0*fVM12 
fi-lABEL 210 
li 

*LABEL 300 
f IDOr^IDGEM 
' TEM2=^ ID0F -*13 
^JNZ (TLM2, 400) 
f TEM2--FREE ( ) 

I IDUP=FREE( ) 

' 312 1- FREE ( ) 

! S t54=FREE( ) 

! SI87-FREE< ) 

^ 310l=FRr:E( ) 

!311 -4"REE() 

•SI2 =FRFE() 

!3I3 ---FREEO 
•314 =FREE() 

!S15 =FRLE( ) 

•3 16 =FRFE() 

*SI7 -FRbE() 

!B1B -FREE() 
fSI9 - r-REEC) 

!SI 10=FRFE< ) 

!SJ1 l-FREE( ) 

*S1 12=f-RBE( ) 


original 
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‘L-;pi ^*FRir:E() 

't'P2 -FRFE() 

'EP3 --=Fr<f.E<) 

!EP4 ^^PREE() 

'EPS ^FRtP() 

!P.P6 :=FREF() 

‘FF7 :=FRFFv) 

'EPS =^rRFF( ) 

IFF 9 ==FRFF() 

'EP10=^FREF( ) 

! E.Pl t =-FHEF( ) 

!E:P12^^FRFE( > 

' VMi^*FRbF( > 

• VM2-^FR£E< ) 
f V)43^rRL"F< ) 

! VM4=^PHCF< ) 

!VMS=FHEF(> 

! VHft^FRFE < ) 

! VN?^PRFF( ) 

I VNS=.4-R£F( ) 

‘ VN9-^FrtEE< ) 

! Vh I \0=-- FREE ( ) 

• VHl l-FREt ( ) 

• VH1.^^FREE( ) 

‘ i>FlF = FRFE( ) 

p 

3 iiOW ASSIGN LOAD TO PROPER DEGREES OF FREEDOH 

• VMSl 

DL"F N4C I3B2=16 ELNO CA2 O 0 

rABLE(NI=24,. NJ^^24) :17 ROTA EL EH "ADLC" 99 
! Al^- JEMT" liP36 
! li i ^0 


TRAN3FERR (S0URCE=DB2, 

ILIM==3. 

SBASE=-"A1 “i 

DBASE="B1 

! Al=^ v>ENT -1 5P36 V3 
' U.t -24 

f N ANSFERR ( B0URCE^BB2, 

ILIM^‘3, 

SBASE^”A1*‘, 

DBASF>"B1 

! A 1 — 'JEN 1”* 1 *36+6 
! B J -“43 

PRANSFERR < St3URCE=DB2. 

ILIM=3, 

SDASE=’‘A1 

DBASE=“B1 

‘ Al=JENT-l*36+9 
' Bi-72-*3 

1 RAN3FERR ( S0URCE“-^BB2i 

ILiri^^3> 

SBASE=”A1 "i 

D0ASE=“B1 

' A1=JENT -1*36M2 
131-96+3 

fRANSFERR ( S0URCE^BB2> 

ILIH=3> 

SBASE="A1 'S 

DBASE=*'B1 

! Ai==jeNT-l*36+lS 
!Bl==).20+3 

TRAN3FERR (S0URCEr.=BB2> 

ILiri='3> 

SBASE=”Al“i 

. DBASE="B1 

' Ai-^s3EN r-l*36+18 
»R1.^1444,6 

TRANSFERR < S0URCE^BB2i 

ILIH==3> 

SBASF>'*A1 

> DBASE==‘*B1 

! A1==. JENT-* 1*36*21 
' 01 = 160+6 

TRAfJSFERR (S0URCE = DB2> 

ILIH=3, 

SBASE= “Ai '9 

> DBASE="B1 

! A1=.JENT -1*36+24 
*01=192+6 

■fRANSFERR (S0URCE=BB2. 

ILIH-=3, 

SeASE='**Al“; 

• DBASE=”B1 

! Ai -^ JFNT - 1*36+27 
‘01=216 1-9 

TliANSf-ERR ( B0URCE=BB2. 

IL,IH=3> 

SBASE="A1 

. DBASE='*B1 
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‘ A1 ♦^JEN r-l >36-S“30 
I 

Ti^ANL'FERR (S0URCE=BB2, 
! Ai=-- JEN r-‘1^^36>3:3 
‘Bi^264-^•9 

THANGFERR <SOUHCE--^DI32, 
' 01^288+12 

TRANSFERR ( S0URCE=BB2> 
! A1-- JEN r-1 1^36 *-3 
' B I = 3 1 2 1 2 

i’HANSFFRR ( B0URCE'*^BB2, 
! A1=^J£NI ’-I ’-36+6 
! B i-=^336+12 

TR ANSKFHR ( BDURCE==DB2. 
?Al ~x>FN I -‘1-5^36+9 
‘31-360+15 

TRANGFERR ( B0URCE-DB2; 
! A1=^ JENT-1 «-36+12 
‘ B1=^384H5 

PRANGFERR (30URCE=BI32i 
f Al=^JffWT-l^>36+15 
•31=40S^ 15 

TUANGFEHR (S0URCE=^BB2, 
! A1=-< JFNT-1 }^36+1S 
: U 1=^432+18 

TP ANGFERR < bOURCE-DB2. 
’ Al^JENT-1 »-36+2l 
‘ B 1==456+1G 

n-^ANGFERR (S0URCE-BB2, 

• A 1 =■* JENT- lit 36 +24 
! B 1 H B0-»' 1 8 

TRANSFERR < SCIURCE-BB2, 

• A i JFNT- 1 it36+27 
in 1=^504-^21 

I RANSFERR ( S0URCE=--'BB2, 
' Al=^JENT-lit-36+30 

• B1==52B+21 

TRAN3FERR ( S0URCE==BD2. 

I Al=^JENT-lit36+33 
‘Bi-‘552+:>l 

TRANGFERR (SOURCE- DB2, 
' AI^TREE ( ) 


• Bl=hREE( ) 

TABLE (NI^l. NJ=24 


J= 1 

0. 0 

J= 2 

O 0 

J-^3 

'*AJ)1 

^J--4 

"AO? 


"AU3 

vl^6 

0. 0 

J=.--7 

0. 0 

J^S 

0 0 

J.:.9 

"A 04 

J 0 

"AOS 

1 1 

"A06 

J- J 2 

0 0 

J^i3 

O. 0 

J i4 

0 0 

J-J5 

“AOV 


17 


ILIM=3, 

SBASE^“A1 •*> 

DBASE="B1 


SBASE=“A1 ‘S 

DBASE="B1 

ILIM==3, 

SDASE=^^"A1 ", 

DBASE=^"B1 

ILIM^=3, 

SBASE-"A1 ", 

DBASE="LU 

ILIM=3> 

SBASE=--"A1 ", 

DBASE="B1 

ILin==3, 

SBASE="A1 ", 

DBASE="B1 

ILIM=^3, 

SBASE="A1 ", 

DBASE="B1 

ILiri==3, 

SBASE=^"A1", 

DBASE="B1 

ILIM==3, 

SBASt>"Al*', 

DBASE=--"B1 

ILIM==3. 

SBASE="A1 ", 

DBASE="B1 

ILIM=3i 

SBASE~"A1 ", 

DBASE="B1 

ILIM=*3, 

SBASE-"A1 ", 

DBASE^"B1 

ILIM=3. 

SBASE=^"A1 ", 

DBASE^“B1 

ILIM=-3, 

SBASE^"A1 ", 

DBASE="B1 


EL.LO E42 "ADLC" 99 
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J --16 ; 

"ADO" 

J=-17 : 

"AD9" 

J"-1B : 

0. 0 

J^\9 : 

0. 0 

J^-20 : 

0. 0 

J=21 ■ 

“ADlO 

J--22 : 

“ADI 1 

J-23 : 

"AD 12 

J-24 ; 

0. 0 


DEFINE A I 17 EZLLO L4? **ADLC“ 99 
17 ELLO E^2 '‘ADL.C" 1 =^RTRAN < A1 ) 

DEFINE A1 ^ 17 ROTA EL £h '‘ADLC’ 99 
17 RQ-fA ELEM •'ADLC'' 1 =r<7RAN< A1 ) 

DEFINE 17 ELL.O EA2 "ADl.C*' 1 

DEFINE Al ^ 17 ROTA ELEN “ADLC" 1 
16 ADlO VE42 “ADLC“ O RPROD ( Al> BB2 ) 


‘ADl^DS 1,1, 

1 


(16 ADLO VE42 

"ADLC" 

0) 

' AD2^DS 2.1, 

1 


(16 ADLO VE42 

“ADLC" 

0) 

' AD3==0S 3, 1, 

1 


(16 ADLO VE42 

“ADLC" 

0) 

IAD4=:0S 7, 1, 

1 


(16 ADLO VE42 

"ADLC" 

0) 

‘ADb:=D5 a, 1, 

1 


(16 ADLO VE42 

“ADLC" 

0) 

•AD6^=DS 9, 1, 

1 


(16 ADLO VE^2 

"ADLC" 

0) 

•AD7=DS 13.1,1 


< 16 ADLO VE42 

“ADLC" 

0) 

•ADvO^-DS 1M,1 

. 1 


<16 ADLO VE42 

“ADLC" 

0) 

!AD9=-DS 15. 1 

. 1 


( 16 ADLO VE42 

"ADLC" 

0) 

•ADIO^-DS 19. 

1. 

1 

(16 ADLO VE42 “ADLC 

“ 0) 

•ADll-DS 20. 

1, 

1 

(16 ADLO VE4 

2 “ADLC 

“ 0) 

!AD12:=DS 21. 

1. 

1 

(16 ADLO VE4 

2 “ADLC 

“ 0) 


% 

•J1=DS 13> “JE-:NT“, 1 (16 CONN E42 0 0) 

! J2=D3 14. “JENT". 1 <16 CONN E42 0 0) 

!J3=^0S 15, “JENT", 1 (16 CONN E42 0 0) 

!J4=-DS 16, “JENT", 1 (16 CONN E42 0 O) 

3VSVEC: APPL FORC "ADLC" 1 

I^-l 


J^“ J1 “ 

“ADI " 


“AU4 “ 

J3“ 

“AD7" 

'• J4» 

“ADlO" 

1 "--2 


J^.“ J1 “ 

“AU2" 


"aD5'* 

J-" J3" 

“ADQ" 

J4" 

"ADI 1" 

I- 3 


J=u*' jj •• 

“A) >3" 

J2" 

“AD6" 

J3" 

“AD7" 

J^“ J4“ 

“AD 12" 


AD1=DS 

4, 1. 1 

<16 

ADLO 

VE42 

"ADLC" 

0) 

AD2=DS 

5. 1. 1 

(16 

ADLO 

VE42 

"ADLC" 

0) 

AD3=DB 

6. 1. 1 

(16 

ADLO 

VE42 

"ADLC" 

0) 


A04^^DS 

10, 1, 1 

< 16 

ADLO 

VE42 

"ADLC " 

0) 

Ai>5^^i:>3 

11.1.1 

( 16 

ADLO 

VE42 

"ADLC" 

0) 

7iD6~DG 

12. 1. 1 

( 16 

ADLO 

VE42 

"ADLC" 

0) 

AD7=DS 

16. 1. 1 

< 16 

ADLO 

VE42 

“ADLC" 

0) 

ADB-^DS 

17. 1. 1 

( 16 

ADLO 

VE42 

“ADLC" 

0) 

AD9=DS 

18. 1. 1 

( 16 

ADLO 

VE42 

“ADLC" 

0) 


•ADIO^DS 22.1.1 <16 ADLO VE42 “ADLC" 0) 
lADll^DS 23.1,1 (16 ADLO VE42 “ADLC" 0) 





1. 1 (16 ADLO VE42 "ADLC“ 0) 

1-4 


~ ** J 1 '* 

“ADI 

“ j:>" 

•‘AD4" 


••AD7‘* 


••ADKV 



ji “ 

“A02** 

J-“ J2*’ 

•*AD‘3'* 

J - Jj“ 

“ADO" 

•J- “ J4“ 

“ADI t" 

1-6 



“Al)3“ 


“AD6" 

•' J3“ 

“AD9“ 

J-" J4 “ 

“ADI 2" 


‘ AD l-^rKb’E( ) 

! AD2^^KRL-ri::( ) 

‘ AD3--FREE( ) 

I AD4-'-^rRL-E( ) 

’ AD!5^-FREE( ) 

! AD6=^FREE( ) 

I A07--bRL: lf ( ) 

' Ar>8- h Rlr E< ) 

» Ai:>9-^FHte:< ) 

! A01O^FRF.E( ) 

‘ ADI 1=^!'REE( ) 

* ADi;?---^Frvl-E( ) 

' JJ ^-KNFE( ) 

!v;2-bRhE() 

‘ J3~FREE( ) 

IJ4=^FRIt'.E< ) 

^XQT DCU 
ERASE 17 
-f(ETURM 

EN):jAD42 
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^ 

4 This t^li 0 a subroutine in EAL Pinal starts running the 
V proyraf/j Everything before is stored as subroutines in the 
'*> EAL ii at abase. 
h 

-OCAl.L. (29 DR IV GLOB O 0) 

% 

^XQl EXIT 
f r ME 



Addendun to Technical Report No. 86-2 

RESULTS OF THE DESIGN SENSITIVITY ANALYSIS 

The results of the design sensitivity analysis program in EAL are all 
stored in EAL-llbrary file L012. 

If a stress constraint design sensitivity of the appropriate element type 
is specified in the input control parameter, the following data sets are 
created to store the appropriate stress constraint values: 

DVAL E21 1 0 The values of the stress constraint 

functionals are given for constraints 
listed in the input data set ST 21 LIST 

DVPO E21 1 0 A set of pointers that indicate the location 

of points where maximum stress occurs in 
elements are given for constraints listed in 
input data set ST21 LIST. 

DVAL E41 1 0 The values of the stress constraint 

functionals are given for constraint listed 
in input data set ST41 LIST 

DVAL E42 1 0 The values of the stress constraint 

functionals are given for constraints listed 


in input data set ST42 LIST 



The names of the data sets for design sensitivity vectors are given in the 
EAL-library ' file L012 that have the following basic form: 


where 

DSVE 

"ETYP" 


"CONT" : 


« 

DSVE "ETYP" "CONT" "CIND" 


Design sensitivity vector 

- E21, E41, or E42: "ETYP" indicates the element type for which 

the design sensitivity vector is stored. 

For element types E41 and E42, the data set 
contains one sensitivity vector, indicating 
design sensitivity with respect to thickness 
design parameter, t. For element type E21, 
the data set contains two sensitivity 
vectors, indicating design sensitivity with 
respect to width, b, and height, h, of the 
beam. 

In the input file for each element, a 
control parameter has to be specified as the 
following arrays: 


ED21 

REL 

0 

0, 

ED41 REL 0 

0 and/or 

ED42 

REL 

0 

0. 

The parameter 

"CONT" 


i 




indicates the control parameter for which 
a design sensitivity vector is computed. 


CIND": 


CIND" is the constraint indicator 


10,000 


between 

21,000 

and 

22,000 

between 

22,000 

and 

23,000 

between 

23,000 

and 

24,000 

betweem 

30,000 

and 

40,000 

between 

40,000 

and 

50,000 

between 

50,000 

and 

60,000 


Compliance Constraint 

Displacement constraint in the x-direction 

Displacement constraint in the y-direction 

Displacement constraint in the z-direction 

Stress constraint in beam elements 
Stress constraint in membrane elements 
Stress constraint in plate elements 


For displacement and stress constaints, the last three digits of the 
constraint indicator "CIND" give the node number of the displacement 
constraint and the element number of the stress constraint, respectively. 


